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This paper summarizes recent developments in the theory of Bogomol'nyi-Prasad- 
Sommerfield (BPS) state counting and the wall crossing phenomena, emphasizing in 
particular the role of the statistical mechanical model of crystal melting. 

This paper is divided into two parts, which are closely related to each other. In the 
first part, we discuss the statistical mechanical model of crystal melting counting BPS 
states. Each of the BPS state contributing to the BPS index is in one-to-one correspon- 
dence with a configuration of a molten crystal, and the statistical partition function of 
the melting crystal gives the BPS partition function. We also show that smooth geome- 
try of the Calabi-Yau manifold emerges in the thermodynamic limit of the crystal. This 
suggests a remarkable interpretation that an atom in the crystal is a discretization of 
the classical geometry, giving an important clue as to the geometry at the Planck scale. 

In the second part we discuss the wall crossing phenomena. Wall crossing phenomena 
states that the BPS index depends on the value of the moduli of the Calabi-Yau mani- 
fold, and jumps along real codimension one subspaces in the moduli space. We show that 
by using type IIA/M-theory duality, we can provide a simple and an intuitive derivation 
of the wall crossing phenomena, furthermore clarifying the connection with the topo- 
logical string theory. This derivation is consistent with another derivation from the wall 
crossing formula, motivated by multi-centered BPS extremal black holes. We also explain 
the representation of the wall crossing phenomena in terms of crystal melting, and the 
generalization of the counting problem and the wall crossing to the open BPS invariants. 
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1. Introduction and Overview 

One of the most fundamental problems in theoretical physics in the 21st century 
is to construct a theory of quantum gravity. General relativity and quantum me- 
chanics, which are the cornerstones of the 20th century physics, are mutually in- 
consistent. For example, general relativity is not renormalizable, and conventional 
perturbative techniques in quantum field theory does not apply. Not only are these 
two theories inconsistent, but there are indications that theoretical consistency re- 
quires unification of these two theories. In general relativity, there is a celebrated 
singularity theorem which predicts the existence of singularities where general 
relativity breaks down. In this sense general relativity is not a consistent theory in 
itself, and singularity occurs precisely when quantum effects are expected to play 
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crucial roles. The existence of singularities is classically not a problem if cosmic 
censorship conjecture ^ holds and singularity is surrou nded by the horizon. How- 
ever, quantum mechanically black holes emit particles ^El. This then leads us to 
the paradox of information loss ^ (see Refs. HI [7] for review), since we only have a 
thermal radiation coming out of the event horizon regardless of what goes inside the 
horizon. One should always keep in mind, however, that these arguments are based 
on semiclassical quantization of general relativity, and proper treatment of these 
problems require a theory of quantum gravity which unifies the general relativity 
and quantum mechanics in a consistent framework. 

Over decades string theory has been the most promising candidate for quantum 
gravity. Just as in general relativity, string theory is severely constrained by various 
consistency conditions, and it is often impossible to tune the parameters by hand. 
This means that if we can extract quantitative predictions from string theory, we can 
perform pass/fail test of string theory as a theory of quantum gravity. Computation 
of black hole entropy is an excellent example for such a stringent test. One of 
the most successful predictions of string theory, as shown by Strominger and Vafa 
in 1996^, is that string theory correctly reproduces the semiclassical Bekenstein- 
Hawking entropy I 9 | 10 | H | 12 | 13 | Q j a c j ags f supersymmetric extremal black holes 
(see e.g. Refs. HU [15] for review). In statistical mechanics, entropy is given by the 
logarithm of total degeneracy of states, and Strominger and Vafa showed that string 
theory reproduces the correct number of states in the limit of large charges. Their 
analysis has subsequently been generalized to many other black holes, including 
four-dimensional black holes and non-extremal black holes. Most remarkably, string 
theory now reproduces not only the semi-classical Bekenstein-Hawking entropy of 
general relativity, but also the subleading entropy c ontributions from the higher 
curvature corrections to the Einstein-Hil bert actio n | 16 | 17 | 18 LL9| ( see j^ e f ^ for 
review), as dictated by the Wald's formula [ 21 | 22 | 23 | This gives a rather remarkable 
quantitative check of string theory as a theory of quantum gravity. 

However, there are many issues that remain to be solved. For example, the 
entropy is typically determined only by the asymptotic growth of the microstate 
degeneracies (given by the Cardy formula ^21) ; in the limit of large charges. How- 
ever, we hope that string theory gives a more complete and detailed theory of the 
microstates, not just the asymptotic growth of their numbers. This will hopefully 
lead to rich and yet unknown aspects of quantum gravity. 

This is closely related to the question: what is the geometry at the Planck scale? 
One of the key ideas in general relativity is the "geometrization of physics" , where 
the physics notion (e.g. mass) are translated and reformulated in terms of geometry 
(e.g. curvature). If we follow a similar path, the central question in quantum gravity 
should be to identify the "quantum geometry" , geometry at the Planck scale where 
spacetime itself fluctuates quantum mechanically. 

In this paper, we will make small steps towards these ambitious goals. Un- 
fortunately, solving string theory in gravity backgrounds is a notoriously difficult 
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problem. The strategy we take is to simplify the problem — to replace the problem 
of gravity with a problem of gauge theory. In string theory compactifications, this 
corresponds to taking the Calabi-Yau manifold to be non-compact. Of course, the 
notion of a black hole is subtle in this limit since the gravity decouples and the New- 
ton constant becomes zero. However, part of the important data in gravity theory 
still remain. For example, we can still discuss entropy of black holes since we can 
take a scaling limit where the mass of the black hole goes to infinity, while the en- 
tropy is kept finite. Moreover, it is generally believed that microstate degeneracies of 
black holes are captured by the BPS index, which is deformation invariant and stays 
constant in the decoupling process. The counting of black hole microstates is now 
turned into a counting problem of BPS states in supersymmetric gauge theories. 

The counting problem of BPS states in string theories and in supersymmet- 
ric gauge theories is an important problem in itself, even if we forget about the 
motivation from black hole physics. For example, they provide primary tools to 
test various dualities in sup ersymmetric gauge theori es a nd in string theory (e.g. 
non-perturbative checks of Montonen-Olive duality for TV = 4 supersymmet- 
ric Yang-Mills theory). BPS solitons in supersymmetric gauge theories have a rich 
structure, and provides a classic example of frui tful co llaboration bet ween physics 
and mathematics, as exemplified in the ADHM I 27 | 28 | anc j Nahm ^3 construction. 
Furthermore, as we will see in later sections BPS state counting has an intimate 
connection with another counting problem in string theory, the topological string 
theory (mathematically Gromov-Witten theory). 

In the first part of this paper, we show that when X is a toric Calabi-Yau 
manifold (this roughly means that X has an action of the three-dimensional complex 
torus; see Refs. l30ll3T1 for review), we can give explicit answers to the BPS counting 
problem. More precisely, each of the BPS states contributing to the BPS index 
(defined in Section ^. 1[) is in one-to-one correspondence with a configuration of a 
molten crystal, and the BPS partition function Zbps (defined in Scction B.ip is the 
same as the statistical partition function of a crystal melting model ^* 

ZbPS = ^crystal- (1.1) 

Scction[3]is devoted to the explanation of this result. Remarkably, the derivation 
of the formula above depends on the new ly dev eloped mathematical theory, the non- 
commutative Donaldson-Thomas theory^EH This theory gives a new invariant for 
Calabi-Yau manifolds, which exactly coincides with the BPS index we are interested 
in. This means that BPS counting problem is important not only to physicists but 
also to the mathematicians alike. 

In the next section (Section|4|), we discuss the implication of these results to 
quantum gravity. We show that the thermodynamic limit of the crystal gives a 
projection of the shape of the mirror Calabi-Yau manifold ^1 This in particular 
suggests that if we start from a classical smooth geometry and goes to smaller and 
smaller scales, the geometry gets discretized into a set of atoms when we finally 
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arrive at the Planck scale. In this sense an atom in the crystal melting model is 
"an atom/quark of internal geometry", a discrctized version of the geometry at the 
Planck scale. We therefore see that the two problems posed earlier are now related. 
Each of the microstate, which is an atom of the crystal, is the discretized version 
of the geometry; thus the problem of identifying black hole microstates is solved by 
the quantum structure of geometry! 

In the second half of the paper, we move on to the discussion of the wall crossing 
phenomena (see Section [2~2l for an introduction) . Wall crossing phenomena states 
that the BPS degeneracy jumps as we change the value of the moduli of the Calabi- 
Yau manifold. Wall crossing phenomena, first discussed in Ref. [36] in the context 
of supersymmetric TV = (2,2) theories in two dimensions, have a long history of 
nearly two decades. They also play import ant roles in the Seiberg-Witten theory 
| 37 | 38 | 39 | 40 | and its stringy realization SDH ( see also Refs.HHIHas wel1 as Ref. US]). 
In these old days, it was possible to derive the jump of BPS states in simple cases, 
but generalization seemed to be difficult. 

The re cent breakthrough was partly triggered by the paper of Kontscvich and 
Soibelman who proposed a rather general formula for the jumps of BPS de- 
generacies (or mathematically "generalized Donaldson-Thomas invariants"), gen- 
eralizing t he " semi-primitive" wall crossing formula previously proposed by Dencf 
and Moore S3 Physical interpretations of the formulas were subsequently given in 
Refs.HS, 49, 50]|nil[511 We discuss these formulas in Section[H]and in |Appcndix A. 3 

The wall crossing formulas arc rather general, and can be applied to our setup, 
namely compactification on the toric Calabi-Yau manifold. In particular, the exam- 
ple of the resolved conifold is analyzed in Ref.[53]and independently in Ref. [54] There 
it was shown that the non-commutative Donaldson-Thomas invariants (which are 
just generalized Donaldson- Thomas invariants in a certain chamber) discussed in 
Scction[3] is related by wall crossing to the commutative (i.e. ordinary) Donaldson- 
Thomas invariants. In physics language, this means that the crystal melting parti- 
tion function is related by wall crossings to the topological string partition function. 
This clarifies the connection of our crystal melting model and another crystal melt- 
ing model proposed in Ref. ESI which describes the topological vertex 

This is not the end of the story, however. First, it was observed in the litera- 
ture that the BPS partition function computed by the non-commutative Donaldson- 
Thomas theory takes a beautiful infinite product form, and there should be an intu- 
itive explanation of this result. From the viewpoint of non-commutative Donaldson- 
Thomas theory, this seems miraculous: we first compute the BPS indices separately 
by intricate mathematics, and only after summing up all of them and going through 
combinatorics can we see that the partition function takes such a simple form. Sec- 
ond, the partition function obtained in Section[3] is very similar to the topological 
string partition function, and it is necessary to understand the precise relation and 
to explain why the topological string partition function enters into the story. 

In Section[5] we give a simple, yet another derivation of the wall crossing phc- 
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nomena from the viewpoint of M-theory This answers the questions raised in 
the previous paragraph. By lifting t ype H A branc configurations to M-theory and 
by using the 4d/5d correspondence l d° l ' 3 ^ the problem of counting BPS states is 
mapped (under certain conditions explained in Scction[3]) to a counting problem of 
free M2-brane particles in five dimensions, which span the free particle Fock space. 
This naturally explains the infinite product form of the BPS partition function. Also, 
the counting proble m of M 2-brane particles is a generalization of the Gopakumar- 
Vafa (GV) argument EQEI] which explains the appearance of the topological string 
partition function. More precisely, we prove a formula which schematically takes a 
form 



Zbps - Z top 



i b ' (L2) 

chamber 



This derivation is consistent with another derivation from the wall crossing formula. 
As a bonus, we have new mathematical predictions for non-toric examples, which 
can be tested by future mathematical developments. 

There is a generalization of the above-mentioned story to the case of open BPS 
invariants. Closed BPS invariants discussed up to this point are defined by count- 
ing D2-branes wrapping 2-cycles of the Calabi-Yau manifold. Open BPS invariants 
are defined by counting D2-branes wrapping on disks ending on another D-branes 
(D4-branes). Open BPS invariants are natural generalizations of closed invariants. 
Moreover, they give a useful computational tool to study closed in var iants for com- 
plicated geometries, as shown in the topological vertex formalism ESI. 

In Section[7] we first give a definition of the "non-commutative topological ver- 
tex" , which gives a basic building block for computing open BPS invariants. The 
definition uses the crystal melting model, and we perform several consistency checks 
of the proposal. We also discuss the wall crossing phenomena for the open BPS in- 
variants both with respect to the open and closed string moduli, by again using the 
viewpoint of M-theory. 

In the final section (Section[8]), we close this paper by pointing out several inter- 
esting problems which suggest directions for future research. We also collect slightly 
technical results in the appendices. 

This paper is essentially the author's Ph.D thesis (submitted to University of 
Tokyo on November 2009), which in turn is mostly based on his papers Refs.l32 j l35 | 
[57l 162"! l63l We tried our best to make the presentation consistent and to present a 
unified perspective on the problem. However, BPS state counting and wall crossing 
phenomena are very active areas of research with long history, and is still growing 
very rapidly as of this writing. This inevitably means that there are many important 
omissions. For example, in this review there is almost no discussion of the wall 
crossing phenomena from 4d M = 2 gauge theory viewpoint (see e.g. Refs HH 
HUGO] in thi s direction), and our discussion is mostly on the 6d Calabi-Yau side. 
Moreover, in order to keep this paper to a reasonable length we could not include 
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enough introductory material in this review. The complete introduction will require 
another review, but the good starting point is the lecture note Ref. [64j 



The interdependence of the sections are shown in Fig. [I] 




§A.2. Multi-Centered Black Holes 



§A.3. Kontsevich-Soibelman Formula 



Fig. 1. The interdependence of the sections of this paper. 



For the convenience of the reader, each section begins with a short summary 
of the section. The connections between different sections are emphasized through- 
out this review in order to organize the whole material into a unified framework. 
However, some sections are mostly independent, and do not require prior reading of 
earlier sections. For examples, most of section 5 does not require the prior reading 
of part I. 



2. Basic Notions 

Before venturing into the study of crystal melting and wall crossing phenomena 
beginning Section|31 we here give a brief summary of the BPS states, BPS indices 
and their wall crossing phenomena. The material in this section is standard, and 
can safely be skipped by readers familiar with these topics. 



2.1. BPS States and BPS Index 

In theories with extended (J\f > 1) supcrsymmetry in four dimensions, the super- 
symmetry algebra has a central extension Zmn(M, N = 1, . . . , AT), which appears 
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in the commutation relation of two supercharges Q^ 1 (in the two-component nota- 
tion) : 

{Qa \Qp} — CapZMN- (2-1) 

From the above equation it follows that Zmn is an antisymmetric matrix. In the 
case of A/" = 2 supcrsymmctry treated in this paper, Zmn has just one component 
(a complex number) Z = Z\2. We can then show, purely from symmetry arguments, 
that the mass M of a particle always satisfies the BPS inequality @ 

M>\Z\. (2.2) 



States which saturates this inequality are called BPS states. 

BPS states have played special roles in the study of supersymmetric gauge the- 
ories and string theories. This is because BPS states are protected from quantum 
corrections. In other words, quantum corrections may renormalize the value of Z, 
but the BPS inequality holds for all values of the coupling constant. This is ba- 
sically due to the fact that BPS states are annihilated by certain combination of 
supercharges and thus their representation have smaller dimensions (i.e. they belong 
to short representations) as compared with non-BPS states (which belong to long 
representations). Therefore BPS states provide an indispensable tool to study the 
strong coupling behavior of supersymmetric gauge theories. They are also impor- 
tant for understanding black hole microstates, and also are related to mathematical 
invariants, as discussed in the introduction. 

Let us try to be more specific. In this paper, we arc going to study the counting 
problem of BPS states arising from type IIA string theory compactificd on a Calabi- 
Yau manifold X . Depending on whether X is compact or non-compact, we have a 
4d J\f = 2 supcrgravity or 4d J\f = 2 supersymmetric gauge theory. 

We want to study BPS states which preserve half of the Af = 2 supcrsymmetry 
in the theory. These states span the BPS Hilbert spaced 

%BPS7 

and we want to study the structure of this Hilbert space. 

The BPS states have charges 7, which takes value in the (Z- valued) cohomology 
of X 

7 G T = H 6 {X; Z) © H\X; Z) © H 2 (X; Z) © H°(X; Z), (2.3) 
and this means that the BPS Hilbert states can be decomposed into sectors with 



a This inequality was originally discovered by Bogomol'nyi, Prasad and Sommerfield IS5ISSI in the 
context of topological solitons, and was later sho wn t o arise from the supersymmetry algebra 
b This space is equipped with a product structure EE] 
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specific charge 7: El 

Hbps = 0Hbps, 7 . (2.4) 

7 

We define the BPS index fi(7) as an important quantity to understand the 
structure of Kbps- This index is defined in such a way that the index does not 
jump when two short multiplcts combine into a long multiplet. In four-dimensional 
TV = 2 supersymmetry the answer is unique, and is given by0 

^(7) = -^Tr WBPS , 7 (-I)^(2J 3 ) 2 , (2.6) 

where J3 is the 3rd component of the 5*0(3) spatial rotation symmetry. Also, the 
factor — I is just a convention, which will turn out to be useful later. 

We can calculate the index explicitly for an M = 2 SUSY multiplet. For our 
purposes we only need a massive multiplet, and a massive multiplet of J\f = 2 
theory is either 

(1) a massive long representation, 

Lj : [j]®([1]©4[1/2]©5[0]), (2.7) 

or (2) a short BPS massive multiplet 

S 3 -:[7]®([l/2]©2[0]), (2.8) 

where 5*0 is a hypermultiplet and S1/2 is a vectormultiplet. Their contributions to 
the BPS index are 

fi(5 i ) = (-l) 2j '(2j + l), n{L j )=Q. (2.9) 

This explicitly verifies that the long multiplets do not contribute, and that hypcr- 
multiplets and vectormultiplets contribute 1 and —2 to the index. 

Every BPS particle in an M = 2 theory has a universal half-hypermultiplct 
obtained from quantizing the fermionic degrees of freedom associated to its center 



c More precisely, the RR-charges of D-branes are classified by K-theory (Refs. 1691 1701) - For our 
purposes, however, cohomology is sufficient. 

d The natural candidate for index is the n-th helicity supertrace defined by 

Q„( T ) = (_l)"/ 2 l T r WBpS7 (-l) aJ a(2J s ) n . (2.5) 

However, there are several restrictions one has to worry about. First, n should be an even integer 
because we have CPT invariance and CPT send J3 to — J3. Second, for n = the helicity supertrace 
vanishes, and if n is large then not only the short multiplcts, but also the long multiplcts, begin 
to contribute. In general, if we have four-dimensional M supersymmetry, the only possibility is an 
even integer n satisfying M < n < 2N ', and for M = 2, n = 2 is the only possible option. The 
quantity ^(7) in the main text is ^2(7) in this notation. The situation is different, for example, 
for = 4. In this case, index is not unique and is given either by ^4(7) or fig (7), or their linear 
combination. See appendix E of Rcf. 1711 for a nice summary. 
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of mass (sec (|2~51)). We can therefore write "H B ps, 7 = ([1/2] ©2[0]) ®%bps,7- Writing 
J3 for % B ps 7 we have 

n( 7 ) = -2 Tt <ps7 (-l) 2J 3 (^2(J 3 ) 2 - (J 3 - i) 2 - (J 3 + \f 

The final result is simply that we have the Witten index Tr(— 1)^, the number 
of boson minus the number of fermions, where bosons/fermions refer to the spin j 
in (|2.7|) . This quantity is important since in many examples it is this index which 
gives the correct microstate degeneracies of black holes 

Let us define the generating function of the BPS index. As we will see later in 
examples, it is often more useful to study the generating function rather than to 
study each of the BPS index separately. Decompose the charge lattice T into electric 
and magnetic charges lattices: 

r = r e ©r m . 

In type IIA compactification of a Calabi-Yau manifold, we have 

r e = H 6 (X;Z) ffii7 4 (A;Z), Y m = H 2 (X; Z) © H°(X; Z), (2.11) 
and when we write 

7=(q ,qi,p I ,P°) G H 6 (X;Z) © H 4 (X;Z) © 77 2 (A;Z) © 77°(X;Z), (2.12) 

the electric charges are q\ = (qo, qi) and the magnetic onesp A = (p 1 ,p°). Here 7 and 
A run over 7 = 1,..., dimi7 4 (A; Z) and A = 0, . . . , dim i7 4 (A; Z), respectively. 
We define the generating function Zbps by a 

Z BPS (cb A ,P A ) = E n (7 = (<ia,p a )) e~*V (2-13) 

This means that we are considering microcanonical ensemble of magnetic charges p A 
and a canonical ensemble for electric charges q\ with chemical potentials </> A . Such 
a mixed microcanonical/canonical ensemble appears in the Ooguri-Strominger-Vafa 
conjectured^. In particular, it was shown there that 

T(4> K ,P K ) ■= log (Z BPS (0 A ,P A )) (2.14) 
is the Legendre transform of the entropy <Sbh(<ZAiP A ) of the black hole 

Sbh(<ZA,P A ) = H4> A ,P K ) - A ^^(0 A ,P A ). (2.15) 

For this reason, Zbps is sometimes denoted Zbh m the literature. 



There are some discussions in the literature as to whether it is the index or the total degeneracy 
which gives the correct black hole entropy. See Ref . 1721 for a recent discussion. 
For compact examples, we necessarily have dim H 4 (X;1,) = dim_ff 2 (X;Z) because of Poincare 
duality. These can be different in non-compact examples we discuss later. 
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2.2. Wall Crossing Phenomena 

In the previous section we defined the BPS index ^(7) and its generating function 
^bps(0 A ,P A )- From the viewpoint of counting BPS states all we need to do is to 
determine these numbers. However, there is an interesting twist, the wall crossing 
phenomena. The wall crossing phenomena states that the BPS index ^(7) depends 
on the value of moduli of the Calabi-Yau manifold at infinity @. More concretely, 
consider the Kahlcr moduli of the Calabi-Yau manifold. Then wall crossing phe- 
nomena states that the BPS index jumps along the real codimension 1 subspace of 
the Kahler moduli space. The codimension 1 subspace is often called the walls of 
marginal stability (or sometimes simply walls), and the walls of marginal stability 
divides the Kahler moduli space into several regions, called chambers. The BPS 
index is piecewise constant inside each of the chambers, but jumps as we cross the 
walls of marginal stability. We sometimes use the notation ^(7; t) to explicitly show 
the dependence on Kahlcr moduli. 

What is the physics behind such a jump? Naively, the jump seems to be in- 
consistent with the common wisdom that index stays invariant under continuous 
deformations. The resolution to this puzzle is given by the fact that Hbps-j used in 
the definition of the BPS index (|2.6|) is in fact a one-particle Hilbert space. 

Consider the following situation: there is a stable particle with charge 7 when 
the Kahlcr moduli is t = t±. When we change the Kahler moduli to t = £2, the 
particle becomes unstable and decays into two particles with charges 71 and 72 (see 
Fig. [2] Of course, the charge conservation dictates that 7 = 71 + 72.). This means 
that the one-particle Hilbert space at t = t\ is larger than that at t = t%. This 
explains the jump of the BPS index: although the index is deformation invariant, 
the Hilbert space over which we are taking the trace jumps as we move around the 
Kahler moduli space 0. Of course, the reverse process is also possible: the stable 
particles can combine to make a new bound state. 



wal I of marginal stabi I ity 




Fig. 2. Decay and formation of BPS bound states under the change of moduli. The values of 
moduli arc different for left and right. 



s We added "at infinity" here in order to emph asize t hat the value of the moduli is not the horizon 
value determined by the attractor mechanism | 73 | 7 4] 

h There is a neat supergravity picture for this, which will be explained in | Appendix A. 2] 
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Now the next question is how to determine the positions of walls of marginal 
stability. In fact, this is easy. Consider the decay 7 — >• 71 + 72- Then the particle 
with charge 7 is stable (unstable) if its mass M( 7 ) is larger (smaller) than the mass 
or particles 71 and 72. In other words, on the walls of marginal stability we have 

M(7) = M( 7l ) + M( 72 ). 

At the same time, since the particles are BPS the mass is given by the central charge 

M( 7 ) = |Z( 7 )|. 

Since the central charge is a linear function with respect to charges, we have 

z{i) = + 2(72). 

The only way to make the above three equations consistent is to have 

Arg(Z(7i))-Arg(Z( 72 )). (2.16) 

This is the equation determining the position of walls of marginal stability. From 
(|2.16[) it follows that the position of walls of marginal stability for the decay process 
7 — >• 7i + 72 is the same as that for the process 7 — > Niji + ^272 ■ 

We can also consider more complicated decay patterns, for example 7 — > Niji + 
A272 + -^373- However, such a decay occurs in a higher codimension subspace of 
the Kahler moduli, which in this case is determined by 

Arg (Z( 7l )) = Arg (Z( 72 )) = Arg (Z(j 3 )) . 

Since these decay patterns are non-generic, we are not going to consider these com- 
plicated decay patterns in the rest of this paper. 

In this section we have formulated the problem: to determine the BPS index 
or its generating function in each chamber of the moduli space. In the following 
sections we will see that this problem is beautifully solved using the crystal melting 
model. 



3. Crystal Melting 

In this section we discuss the crystal melting model proposed in Ref. 1321 The work 
E2l i s motivated by and is based on earlier results in mathematics I 33 | 34 | 



3.1. Overview: Crystal Melting and Topological Strings 

3.1.1. Formulation of the Problem 



We already discussed BPS state counting in Scction l2TTl but here we would like to 
be more specific and would like to prepare some notations used throughout the rest 
of this paper. 
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There arc different versions of BPS counting problems depending on which string 
theory and what kind of branes we consider. As in the previous section, in this paper 
we concentrate on the BPS counting problem of Type IIA string theory compactificd 
on a Calabi-Yau manifold X Q. We emphasize here that the problem in this paper 
is the counting in the physical string theory, not in the topological string theory. 
Topological string theory comes later into the story, somewhat surprisingly (see 
Section[5] for clarification on this point). 

In this section and the following, we assume in addition that X is a toric Calabi- 
Yau manifold. The advantage of this choice is that thanks to the torus action on X 
we can often perform explicit computations, therefore making quantitative state- 
ments possible. Since X is non-compact, the four-dimensional theory is an A/" = 2 
supcrsymmetric gauge theory. 

We consider 1/2 BPS states of this N — 2 theory. In string theory, BPS states are 
realized as supersymmetric BPS bound states of D-branes wrapping holomorphic 
cycles. We concentrate on the configuration of D-branes which arc particles in 4d: in 
other words, we consider Dp-brancs wrapping p-cyclcs. D-brancs which are vortex 
strings and domain walls in 4d are also of potential interest, and we will see such 
examples in the discussion of open BPS invariants in Scction[7] 

In type IIA string theory, we have Dp-branes with p even. We therefore con- 
sider D0/D2-branes wrapping holomorphic compact 0/2-cycles of the Calabi-Yau 
manifold. Due to the slightly technical reasons explained later, we do not include 
D4-branes wrapping 4-cycles. We also consider a single D6-brane filling the en- 
tire Calabi-Yau manifold. The D6-brane is different from other D-branes in that it 
wraps the whole X, which is non-compact (recall toric Calabi-Yau manifold is non- 
compact). We will see later that the D6-brane plays crucial roles in the discussion 
of the wall crossing phenomena. 

Let us denote the D-brane charges of the particles by© 



7 = (n, P, 0, 1) e H 6 (X; Z) © H 4 (X; Z) © H 2 (X; Z) © H°(X; Z), (3.1) 



where /? = {/?/} collectively denotes an element of H (X ; Z). As in the introduction, 
the question is how to compute the BPS index 



'By string duality this is sometimes related to the BPS state counting in different string theories. 
For example, type IIA string theory on K3 X T 2 is dual to Heterotic string theory on T 6 . 
J We are changing the notation from the previous section. The notation here will be used throughout 
the rest of this paper. 



n( 7 ) 



or its generating function 




(3.2) 



where Q = {Q/} and Q 13 := UiQi- 
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3.1.2. The Connection with Topological String Theory 

Before discussing the crystal melting model, let us summarize the known connec- 
tion between the BPS partition function and the topological string theory. In the 
literature, the partition function Zbps we define here is related to topological string 
theory partition function Z top in two different contexts. 

First, as explained in Section l2~T1 the quantity Zbps is the same the black hole 
partition f unc tion Zbh which appears in the context of Ooguri-Strominger-Vafa 
conjecture The conjecture states that when the D-brane charges are chosen 
such that bound states become large black holes with smooth event horizons, the 
generating function Zbh of a suitable index for black hole microstates is equal to 
the absolute value squared of the topological string partition function Z top , 

Zbh = \Zto P \ 2 , (3.3) 

to all orders in the string coupling expansion. 

Second, when there is a single D6 brane with DO and D2 branes bound on it, it 
has been proposed that the bound states arc counted by the Donaldson-Thomas 
invariants I 76 | 7? l f the moduli space of ideal sheaves on the D6 brane. For a non- 
compact toric Calabi-Yau manifold, the Donaldson-Thomas invariants are related 
to the topological string partition function | 55 | 75 | 78 | us j n g topological vertex 
constructional. Recently the connection between the topological string theory and 
the Donaldson-Thomas theory for toric Calabi-Yau manifolds was proven mathe- 
matically in Ref.[7ni Given the conjectural relation between the counting of D-brane 
bound states and the Donaldson- Thomas theory, it is natural to expect the relation, 

Zbh = Z top - (3.4) 

Note that these two conjectures are supposed to hold in different regimes of validity, 
as we will discuss in the final section (Section[5]) . 

From the viewpoint of counting BPS states in type IIA string theory, the param- 
eters q and Q are just formal variables for the generating function. We will see later 
that in the context of topological string theory they are related to the topological 
string coupling constant and the Kahlcr moduli (see (|3.34l) ). 

k An ideal sheaf I is mathematically defined to be a torsion-free sheaf of rank 1 with trivial 
determinant. Given such an X, we can determine the subscheme Y C X (roughly speaking, the 
subspace spanned by the D-brane) by 

Donaldson-Thomas invariants N n p counts the "number" of such ideal sheaves with 

X (Oy)=n, [Y]=fS€H 2 (X,Z), 

where n and /3 are DO and D2-brane charges as defined previously. In this paper we do not use this 
explicit definition of (commutative) Donaldson-Thomas invariants, and assume no prior familiarity 
with them on the side of the reader. 
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3.1.3. What is Discussed in This Section 

The purpose of this section is to understand the case p.4[) better. We start with the 
left-hand side of the relation, namely the counting of BPS states. Recently, the non- 
commutative version of the Donaldson-Thomas t heor y is formulated by Szendroi 
1331 f or the conifold and by Mozgovoy and Reineke 1331 for general toric Calabi-Yau 
manifolds @. In this section, we will establish a direct connection between the non- 
commutative Donaldson-Thomas theory and the counting of BPS bound states of 
DO and D2 branes on a single D6 brane. Using this correspondence, we will find a 
statistical model of crystal melting which counts the BPS states. The formula we 
prove will be written asF°l 

ZbH = ^crystal- (3.5) 

Some readers will be more familiar with the crystal melting description of the 
topological string theory on the right-hand side of (|3.4j) . It was shown in Refs.^SEll 
that the topological string partition function on C 3 , the simplest toric Calabi-Yau 
manifold, and the topological vertex can be expressed as sums of three-dimensional 
Young diagrams, which can be regarded as complements of molten crystals with 
the cubic lattice structure 0. Since the topological vertex can be used to compute 
the topological string partition function for a general non-compact toric Calabi-Yau 
manifold, it is natural to expect that a crystal melting description exists for any 
such manifold. To our knowledge, however, this idea has not been made explicit. 
The crystal melting model defined in this section is different from the one suggested 
by the topological vertex construction, and the precise relation will be explained in 
the discussion of wall crossing phenomena in part II. 

We should next explain how the crystal arises. The short answer is that it arises 
from quiver quantum mechanics B. 

Instead of studying the full dynamics of four-dimensional gauge theory, we can 
use the low energy effective theory on the D-brane worldvolume to compute the 
index. The low energy effective theory of DO and D2 branes bound on a single D6 
brane is a one-dimensional supersymmetric gauge theory, which is a dimensional 
reduction of an TV = 1 gauge theory in four dimensions. Since the index is deforma- 
tion invariant and thus invariant under the process of dimensional reduction, the 
BPS index reduces to the Witten index of this one-dimensional theory, i.e. quantum 
mechanics. The Witten index of quantum mechanics is in turn given by the Eulcr 

'See Refs. |8U|81|54|82|83|84|3^|35|ti2|85|53|86|87|88|8iJ|9U] for further developments . 

111 As will be discussed in detail in Section[5] this is a specialization of the formula 11.21 . 

n See Refs.[9l][92l[l9]|93]|9l|95l[96]for further developments. 

°Onc should keep in mind that the appearance of quantum mechanics (and the computation of the 
BPS index by it) is more universal than in our toric setup, and should apply to any Calabi-Yau 
manifold. The advantage of our setup is that the precise form of the quantum mechanics is known, 
thanks to the brane tiling techniques. The torus action also helps to compute the index of the 
quantum mechanics, since we compute the Euler character of the moduli space by localization 
with respect to the torus action originating from the toric condition. 
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character of the vacuum moduli space This statement will be made more precise 
in Scction [3~2l 

The field content of the quantum mechanics is encoded in a quiver diagram and 
the superpotential can be found by the so-called brane tiling techniquesl 97 | 98 | 99 | 100 | 
0. From these gauge theory data, we define a crystalline structure in three dimen- 
sions. The crystal is composed of atoms of different colors, each of which corresponds 
to a node of the quiver diagram and carries a particular combination of DO and D2 
charges. The chemical bond is dictated by the arrows of the quiver diagram. There is 
a special crystal configuration, whose exterior shape lines up with the toric diagram 
of the Calabi-Yau manifold. Such a crystal corresponds to a single D6 brane with 
no DO and D2 charges. We define a rule to remove atoms from the crystal, which 
basically says that the cry stal m elts from its peak. By using the non-commutative 
Donaldson- Thomas theoryEHEl] we ghow that there is a one-to-one correspondence 
between molten crystal configurations and BPS bound states carrying non-zero DO 
and D2 charges. The statistical model of crystal melting computes the index of 
D-branc bound states. 

The number of BPS states depends on the choice of the stability condi- 
tion, and the BPS countings for different stability conditions are related to each 
other by the wall crossing formulas. In this section, we find that, under a cer- 
tain stability condition, BPS bound states of D-branes are counted by the non- 
commutative Donaldson-Thomas theory. As we will see in Section[5j we can re- 
late non-commutative Donaldson-Thomas theory to the commutative Donaldson- 
Thomas theory by a series of wall crossings. Since the topological string theory is 
equivalent to t he commutative Donaldson-Thomas theory for a general toric Calabi- 
Yau manifold the relation (|3.4p is indeed true for some choice of the stability 
condition, as expected in Ref. [75] In general, the topological string partition func- 
tion and the partition function of the crystal melting model are not identical, but 
their relation involves the wall crossing, 

^crystal melting ~ Z top (modulo wall crossings) . (3.6) 

This does not contradict with the result in Refs. l55l l75l since there is no wall cross- 
ing phenomenon involved between commutative and non-commutative Donaldson- 
Thomas invariants for C 3 . This shows that wall crossing phenomena (to be explained 
in part II) is crucial for understanding the precise relation between crystal melting 
and topological string theory. 

The rest of this section is organized as follows (see Fig. [3]). In Scction l3~2l we will 

p More precisely, the correct statement is 

r2(7) = (-i) dimTAW x(.Mvac), 

where .Mvac is the vacuum moduli space and T.M vac is the tangent space of the moduli space El 
The sign has the effect of shifting the sign of chemical potentials of the crystal melting model, as 
will be explained in more detail in (1 3 . 22 ft in Scction|3]5] 

q See Refs. [ToTI 11021 for reviews of the quiver gauge theory and the brane tiling method. 
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summarize the computation of D-branc bound states from the gauge theory perspec- 
tive. In Section r3.31 we will discuss how this is related to the recent mathematical 
results on the non-commutative Donaldson-Thomas invariants. In Scction[3'"4j we 
will formulate the statistical model of crystal melting for a general toric Calabi-Yau 
manifold, whose partition function is computed in Section[3?5] The final section is 
devoted to the summary of our result and discussion on the wall crossing phenom- 
ena. The equivalence of a configuration of molten crystal with a perfect matching 
of the bipartite graph is explained in appendix | Appendix A.l 



wall crossing 



§3.1. Counting of BPS states 
in type IIA on X 



§3.2. Quiver Quantum Mechanics 



§3.3. Counting of Modules 
(Non-commutative DT theory) 



§5. Topological String on X 



mirror symmetry 



Topological String on the Mirror X 




Thermodynamic Limit 



§3.4. Crystal Melting 



Fig. 3. The organization of this section and the relation to Section|4] and Section[5] 



3.2. Quiver Quantum Mechanics 

Let us begin with the description of quiver quantum mechanics. The first problem 
is to give a precise identification of the quantum mechanics on the worldvolume of 
D-branes. 

In the classic paper by Douglas and Moore^^l, it was shown that the low energy 
effective theories of D-branes on a class of orbifolds are described by gauge theories 
associated to quiver diagrams. Subsequently, this result has been generalized to an 
arbitrary non-compact toric Calabi-Yau threefold. A toric Calabi-Yau threefold X& 
is a T 2 x K-fibration over R 3 , where the fibers are special Lagrangian submanifolds. 
The toric diagram A tells us where and how the fiber degenerates. For a given 
and a set of DO and D2 branes on Xa, the following procedure determines the field 
content and superpotential of the gauge theory on the branes. We will add a single 
D6 brane to the system later in this section. 

3.2.1. Quiver Diagram and Field Content 

The low energy gauge theory is a one-dimensional theory given by dimensional 
reduction of an TV* = 1 supcrsymmctric gauge theory in four dimensions. The field 
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content of the theory is encoded in a quiver diagram, which is determined from the 
toric data and the set of D-branes, as described in the following. A quiver diagram 
Q = (QoiQi) consists of a set Qq of nodes, with a rank Ni > associated to 
each node i £ Qoi and a set Qi of arrows connecting the nodes. The corresponding 
gauge theory has a vector multiplct of gauge group U(Ni) at each node i. There is 
also a chiral multiplet in the bifundamental representation associated to each arrow 
connecting a pair of nodes. 

In the following, we will explain how to identify the quiver diagram. Readers 
are encouraged to consult Fig.|H which describes the procedure for the Suspended 
Pinched Point (SPP) singularity, which is a Calabi-Yau manifold defined by the 
toric diagram in Fig.[D-(a) or equivalently by the equation, 

xy = zw 2 , (3-7) 




Fig. 4. (a) The toric diagram for the Suspended Pinched Point singularity, (b) The configuration 
of D2 and NS5 branes after the T-duality on T 2 . The green exterior lines are periodically identified. 
The red lines representing NS5 branes separate the fundamental domain into several domains. The 
T-dual of DO branes wrap the entire fundamental domain, and fractional D2 branes are suspended 
between the red lines. The white domains contain D2 branes only. In each shaded domain, there is 
an additional NS5 brane. There are two types of shades depending of the NS5 brane orientation. 
The white domains are connected by arrows through the vertices, and the directions of the arrows 
are determined by the orientation of the NS5 branes. (c) The quiver diagram obtained by replacing 
the white domains of (b) by the nodes. 



To identify the quiv er diagr am, we take a T-dual of the toric Calabi-Yau mani- 
fold along the T 2 fibers ESSQH The fibers degenerate at loci specified b y th e toric 
diagram A, and the T-duality replaces the singular fibers by NS5 branes USSl Some 
of these NS5 br anes divide T 2 into domains as shown in the red lines in Fig.|4]-(b) 
| 106 | 107 | 108 [l02 Thc DQ branes become D2 branes wrapping the whole T 2 . The 
original D2 are still D2 branes after the T-duality, but each of them is in a par- 
ticular domain of T 2 suspended between NS5 branes. In addition, there arc some 
domains that contain NS5 branes stretched two-dimensionally in parallel with D2 
branes 0. Let us denote the domains without NS5 branes by i <E Qo an( l the domains 



r The NS5 branes are also filling the four-dimensional spacetime R 1,3 while the D2 branes are 
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with NS5 branes by a G /. In Fig.0]-(b), the Qo _ type domains are shown in white, 
and the /-type domains arc shown with shade. There are two types of shades, cor- 
responding to two different orientations of NS5 branes. This distinction will become 
relevant when we discuss the superpotential. 

The Qo-tyP c domains are identified with nodes of the quiver diagram since open 
strings ending on them can contain massless excitations. The rank of the node 
% € Qo is the number of D2 branes in the corresponding domain. On the other 
hand, /-type domains give rise to the superpotential constraints as we shall see 
below. Though two domains i,j S Qo never share an edge, they can touch each 
other at a vertex. In that case, open strings going between i and j contain massless 
modes. We draw an arrow from i — )■ j or i <— j depending on the orientation of the 
massless open string modes, which is determined by the orientation of NS5 branes. 
Note that the quiver gauge theory we consider in this section are in general chiral. 
This completes the specification of the quiver diagram. 

As another example, the quiver diagram for the conifold geometry has two nodes 
connected by two sets of arrows in both directions. The ranks of the gauge groups 
arc no and no + 712, where no and are the numbers of DO and D2 bra nes H . The 
gauge theory is a dimensional reduction of the Klebanov-Witten theory UQQi when 
?i2 = and the Klcbanov-Strassler theory when 772 > 0. 



3.2.2. Superpotential and Brane Tiling 

Each domain a € I containing an NS5 brane is surrounded by domains i\, 12, i n G 
Qo without NS5 branes, as in Fig.|4]-(b). By studying the geometry T-dual to Ja 
in more detail, one finds that the domain is contractiblc. Since open strings can 
end on the domains ii, 12, i n , the domain a can give rise to worldsheet instanton 
corrections to the superpotential. 

This fact, combined with the requirement that the moduli space of the gauge 
theory agrees with the geometric expectation for D-branes on determines the 
superpotential. Depending on the NS5 brane orientation, the /-type domains are 
further classified into two types, /+ and 7_, and thus the regions of torus is divided 
into three types Qo, I + and i_. Such a brane configuration, or a classification of 
regions of T 2 , is called the brane tiling §. In Fig.rj}(b), the brane tiling is shown by 
the two different shades. The superpotential W is then given by 

w=£Tr( n ^)-E Ti ( n a X ( 3 - 8 ) 

aGl+ \e(Ea, clockwise / aG/- \ e£a, counterclockwise / 



localized along a timclike path in four dimensions. 

s In the context of 4d Af = 1 quiver gauge theories, these branes are called D3 and D5-branes, the 
latter often called fractional branes. 

'In the literature the word brane tiling refers to the bipartite graph explained below. Here the 
word brane tiling refers to a brane configuration as shown Fig.|4}(b). Such a graph is called the 
fivebranc diagram in Ref. II 111 
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where the relation e £ a means that an edge e is one of the boundary edges of the 
face a £ I±, and the product is over all such edges in a (countcr)clockwise manner. 
This formula is tested in many examples. In particular, it has been shown that the 
formula reproduces the toric Calabi-Yau manifold as the moduli space of the 
quiver gauge theory | H2 | 113 | 

In the literature of brane tiling, bipartite graphs arc often used in place of brane 
configurations as in Fig.|U-(b). A bipartite graph is a graph consisting of vertices 
colored either black or white and edges connecting black and white vertices. Since 
bipartite graphs will also play roles in the following sections, it would be useful to 
explain how it is related to our story so far. For a given brane configuration, we 
can draw a bipartite graph on T 2 as follows. In each domain in I + (/_), place a 
white (black) vertex. Draw a line connecting a white vertex in a domain i £ 1+ and 
a black vertex in a neighboring domain j £ i_. The resulting graph T is bipartite. 
See Fig. [5] for the comparison of the brane configuration and the bipartite graph in 
the case of the Suspended Pinched Point singularity. We can turn this into a form 
that is more commonly found in the literature, for example in Ref. I98[ by choosing 
a different fundamental region as in Fig-El 




brane configuration 



bipartite graph 



Fig. 5. The correspondence between the brane configuration on T and the bipartite graph. The 
white (black) vertex of the bipartite graph corresponds to the region /+ {!-) in light (dark) 
shade. The edge of the bipartite graph corresponds to an intersection of I— and in 1+ . From this 
construction, it automatically follows that the graph so obtained is bipartite. 




Fig. 6. By choosing a different fundamental 
commonly found in the literature. 




of T 2 , we find a bipartite graph which is more 
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3.2.3. D-term Constraints and the Moduli Space 

The F-term constraints are given by derivatives of the superpotential, which can 
be determined as in the above. The moduli space of solutions to the D-term con- 
straints is then described by a set of gauge invariant obs eryab les divided out by the 
complexified gauge group Gc l ^^ i The theorem by King UJM states that an orbit of 
Gc contains a solution to the D-term conditions if and only if we start with a point 
that satisfies the (9-stability, a condition defined in the next section. Thus, we can 
think of the moduli space as a set of solutions to the F-term constraints obeying 
the (9-stability condition, modulo the action of Gc- 

3.2.4. Adding a Single D6 Brane 

To make contact with the Donaldson-Thomas theory, we need to include one D6 
brane. Since the D6 brane fills the entire Calabi-Yau manifold, which is non- 
compact, it behaves as a flavor brane. In the low energy limit, the open string 
between the D6 brane and another D-brane gives rise to one chiral multiplct in the 
fundamental representation for the D-brane on the other end. The D6 brane then 
enlarges the quiver diagram by one node and one arrow from the new node. To 
understand why we only get one arrow from the D6 brane, let us take T-duality 
along the T 2 fiber again. The D6 brane is mapped into a D4 brane which is a point 
in some region in T 2 . This means that we only have one new arrow from the new 
node corresponding to the D6 brane to the node corresponding to the D2 branes in 
the region Lj. See Refs. 194 } II 161 for related discussion in the literature. 

3.3. N on- commutative Donaldson- Thomas Theory 

In the previous section, we discussed how to construct the moduli space of solutions 
to the F-term and D-term constraints in the quiver gauge theory corresponding to 
a toric Calabi-Yau manifold with a set of D0/D2 branes and a single D6 brane. 
In this section, we will review and interpret the mathematical formulation of the 
non- commutative Donaldson-Thomas invariant in Refs. [33j [34] for We find that 
it is identical to the Euler number of the gauge theory moduli space. 

3.3.1. Path Algebra and its Module 

For the purpose of this section, modules are the same as representations. Consider a 
set of all open paths on the quiver diagram Q = (Qo, Qi)- By introducing a product 
as an operation to join a head of a path to a tail of another (the product is supposed 
to vanish if the head and the tail do not match on the same node) and by allowing 

"This intuitive picture holds only in the special chamber where the BPS indices coincide with the 
non-commutative Donaldson-Thomas invariants. This is because the basis of fractional D2-branes 
corresponding to the nodes of the quiver diagram arc different in different chambers of the moduli 
space. See Scction EHO 



Crystal Melting and Wall Crossing Phenomena 21 



formal sums of paths over C 3 , the set of open oriented paths can be made into an 
algebra CQ called the path algebra. Note that we are including both gauge variant 
as well as gauge invariant operators. In physics we usually do not consider such a 
huge algebra, but it has turned out to be a useful language for our purposes. 

We would like to point out that there is a one-to-one correspondence between 
a representation of the path algebra and a classical configuration of bifundamcntal 
fields of the quiver gauge theory. Suppose there is a representation M of the path 
algebra. For each node i £ Qq, there is a trivial path e,; of zero length that begins 
and ends at i. It is a projection, (ej) 2 = e^, and satisfies e^a,; = cti for any arrow 
starting from vertex i. Since every path starts at some node i and ends at some node 
j, the sum acts as the identity on the path algebra. Therefore, M = (Bi£Q Mi, 

where Mi = aM. Let us write Ni = dim M, . For each path from i to j, one 
can assign a map from Mi to Mj. In particular, there is an Ni x Nj matrix for 
each arrow i — } j £ Q± of the quiver diagram. By identifying this matrix as the 
bifundamcntal field associated to the arrow i — > j, wc obtain a classical configuration 
of bifundamcntal fields with the gauge group U(Ni) at the node i. By reversing the 
process, we can construct a representation of the path algebra for each configuration 
of the bifundamcntal fields. 

3.3.2. F-term Constraints and Factor Algebra A 

Let us turn to the F-term constraints. Since the bifundamcntal fields of the quiver 
gauge theory is a representation of the path algebra, the F-term equations give 
relations among generators of the path algebra. It is natural to consider the ideal T 
generated by the F-term equations and define the factor algebra A = CQ/J 7 . The 
bifundamcntal fields obeying the F-term constraints then generate a representation 
of this factor algebra. Namely, classical configurations of the quiver gauge theory 
obeying the F-term constraints are in one-to-one correspondence with A-modules. 

As an example, the algebra A for the conifold geometry contains an idempotent 
ring C[ei, e 2 ] generated by two elements and is given by the following four generators 
and relations: □ 




A = C[ei,e 2 ](ai, 02,61,62)/ {aiha% = 026*01, 61^62 = b 2 a i b 1 ) t=12 , (3.10) 

Each A-module for this algebra corresponds to a choice of ranks of the gauge groups 
and a configuration of the bifundamental fields ai,bi satisfying the F-term con- 



F-term constraints have a nice geometric interpretation on the quiver diagram, 
which wc will find useful in the next section. We observe that each bifundamental 
field appears exactly twice in the superpotential with different signs of coefficients 




straints. 



v The center Z(A) of this algebra A is generated by Xij = aibj + bjai(i,j = 1, 2), and is given by 
Z(A) = C[a;ii,a;i2,iz;2i,^22]/(a;iiX22 - £122:21), (3.9) 
which is the ring of functions of the conifold singularity. 
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in the supcrpotcntial shown in (|3.8[) . By taking a derivative of the supcrpotential 
with respect to a bifundamental held corresponding to a given arrow, the resulting 
F-term constraint states that the product of bifundamental fields around a face 
of the quiver on T 2 on one side of the arrow is equal to that around the face on 
the other side. See Fig. [7] for an example. Therefore, when we have a product of 
bifundamental fields along a path, any loop on the path can be moved along the 
path and the resulting product is F-term equivalent to the original one. In Ref. 1341 
it is shown that for any point i, j £ Qo, we can find a shortest path Vi_j from i 
to j such that any other path a from i to j is F-term equivalent to Vi.jU> n with 
non-negative integer n, where a; is a loop around one face of the quiver diagram. 
It does not matter where the loop ui is inserted along the path Vi j since different 
insertions are all F-term equivalent. This means that any path is characterized by 
the integer n and the shortest path Vij. 




Fig. 7. Representation of F-term constraints on the quiver diagram on T 2 . In this example, 
if we write by Xj^g the bifundamental corresponding to an arrow starting from vertex A and 
ending at B etc., then the superpotential 113.81 1 contains a term W = — tr(X^gXgcXcj^) + 
tr(Xj^gXgpXuE^EA)j an d the F-term condition for Xj\g (multiplied by X \g) says that the 
product of bifundamcntals fields along the triangle ABC and that along the square ABDE is the 



In the next subsection, we will impose the D-term constraints on the space of 
finitely generated left A-modules, mod A. Before doing this, however, it is instruc- 
tive to discuss topological aspects of mod A by considering its bounded derived 
category D b (modA). In mathe matic s, the algebra A gives the so-called "non- 
commutative crepant resolution" F or singular Calabi-Yau manifolds such as 
X\ 1 the crepant resolution means a resolution that preserves the Calabi-Yau con- 
dition 0- For a crepant resolution Ya of we have the following equivalence of 
categories @: 

D b (coh(Y A )) = £> b (mod A), (3.11) 



w See Ref, H17l for an introductory explanation of derived categories in the context of string theory. 
x The origin of the name "non-commutative" stems from the fact that the algebra A is non- 
commutative. 

y Mathematically, wc mean a resolution / : — > Xa such that uiy = f*uix> where uix and 
loy are canonical bundles of X and Y, respectively. For the class of toric Calabi-Yau threcfolds, 
the existence of crepant resolut ion is known and different crepant resolutions related by flops are 
equivalent in derived categories ESI. 

z This is well-known in the case of the conifold (cf. Ref. ll20ll . For general toric Calabi-Yau threcfolds, 
this is proven in Ref. [ml See also Refs. [1221 [1231 
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where L> b (coh(Y A )) is a bounded derived category of coherent sheaves of crepant 
resolution Ya, and D b (modA) is the bounded derived categories of finitely gener- 
ated left A-modules. The equation (|3.11[) is also interesting from the physics view- 
point. Since _D 6 (coh(Y"A)) gives a topological classification of A brancs on the 
resolved space Ya, the equivalence means that D b (modA) also classifies D-branes, 
which is consistent with our interpretation above that A-modules are in one-to-one 
correspondence with a configuration of bifundamental fields obeying the F-term 
constraints. 

We should note that the paper I24|, which computes the non-commutative 
Donaldson-Thomas invariants for general toric Calabi-Yau manifolds, requires a 
set of conditions on brane tilings, namely on the supcrpotcntial. It is easy to show 
that the conditions specified in lemma 3.5 and conditions 4.12 of Ref. [34] are au- 
tomatically satisfied for any quiver gauge theories for D-branes on general toric 
Calabi-Yau manifold. To prove the condition 5.3 is more difficult (see Ref. 11241 for 
a proof EB)- 

3.3.3. D-term Constraints and 9-Stability 

We saw that the derived category D b (modA) of A- modules gives the topological 
classification of D-branes in the toric Calabi-Yau manifold To understand the 
moduli space of D-brancs, however, we also need to understand implications of the 
D-tcrm constraints. This is where the 0-stability comes in 0. Let 6 <E be a 
vector whose components are real numbers. Consider an A-modulc M, and recall 
that this M is decomposed as M = (Bi£Q Mi with Mj = aM. The module M is 
called 0-stable if 

]T 0;(dime;M') > 0. (3.12) 
ieQo 

for every submodule M' of A/0- When > is replaced by >, the module M is called 
#-semistablc. 

In the language of gauge theory, the stability condition (|3.12|) is required by the 
D-term conditions. Some readers might wonder why the D-term conditions, which 
are equality relations, can be replaced by an inequality as (|3.12p . In fact, the simi- 
lar story goes for the Hermitian Yang-Mills equations. There instead of solving the 
Donaldson-Uhlenbeck-Yau equations, we can consider holomorphic vector bundles 
with a suitable stability condition, the so-called //-stability or Mumford-Takemoto 

aa Some readers may wonder whether this should be B branes, since it is often stated in the 
literature that coherent sheaves classify B branes. However, in those references we consider D- 
branes filling the whole R 3,1 , whereas in our context D-branes spread only in the time direction in 
R 3,1 . This means coherent sheaves also describe particle-like D-branes in type II A string theory. 

See Refs. IT2T1 |T25 126 for more on consistency conditions in the mathematics literature. 
ac The 0-stability is a special limit of II-stability as discussed in Refs. 127 128 
ad In some literature, an additional condition 5ZigQ &i(dimMj) = is imposed for a choice of 9. 
This is trivially satisfied for the choice 9 = (0, 0, . . . , 0) we choose below. 
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stability H2SH3DI As we mentioned at the end of section 2, it is known that a config- 
uration of bifundamcntal fields is mapped to a solution to the D-term equations by 
a complexified gauge transformation Gc if and only if the configuration is 0-stablc. 
Since each A-module M gives a representation of Gc = iLeQo GL(Ni,C), where 
GL(Ni, C) is represented by Mj = e^M at each node, each A-module specifies a par- 
ticular Gc orbit. Thus, finding a ^-stable module is the same as solving the D-term 
conditions. 

Up to this point we have not specified the value of 9. Physically, 0's correspond 
to the Fa yet-I liopoulos (FI) parameters, which are needed to write down D-term 
equations Although the Euler number of the space of 0-stable A-modules does 
not change under infinit esimal deformation of 9, it does change along the walls of 
marginal stability ^^l The noncommutative Donaldson- Thomas invariant defined 
by Ref. [33] is in a particular chamber in the space of 0's. Following Ref. |34j we 
hereafter take 6 = (0. 0. . . . . 0) R. We will comment more on this issue in the final 
section. 



3.3.4. D6 Brant and Compactification of the Moduli Space 

We have found that solutions to the F-term and D-term conditions in the quiver 
gauge theory are identified with 0-stable A-modules. We want to understand the 
moduli space of such modules and compute its Euler number. 

Since D-brane charges correspond to the ranks of the gauge groups, we consider 
moduli space of (9-stable modules with dimension dim Mj = Ni (i G Qo), which we 
denote by A4 (A). To compute its Euler number, we need to address the fact that 
the moduli space of stable A-modulcs is not always compact. In mathematics liter- 
ature, the necessary compactification is performed by enlarging the quiver diagram 
by adding one more node in the following way. 

Let us fix an arbitrary vertex io, and define a new quiver Q = (Qo, Q\) by 

Qo = QoU{*}, Qi = Qi U {a» : * -> i }. (3.13) 

Namely, we have added one new vertex * and one arrow * — > i$ to obtain the 
extended quiver diagram Q. As in the previous case for Q, we can define the path 
algebra CQ, the ideal JF generated in CQ by J 7 , and the factor algebra A = CQ/J 7 . 
Define 9 G by 9 = (9, 1) and define ^-stable and semistable A-modules using 

stability parameter 9. It is shown in lemma 2.3 of Ref. |34] that (9-semistable A- 
modules are always ^-stable 0, and the moduli space M.^ o (A) of ^-stable modules 
with specified dimension vector N G N^ 0+1 is compact (the boundary of moduli 
space of 0-stable modules are given by ^-semistable modules which are not ^-stable) . 

ac This sounds like a rather degenerate choice, making all module to be 0-semistable but not 0- 
stable according to |I3.12|I . However, as we will see in the next subsection we add an extra node 
corresponding to the D6-brane and use the stability parameter 6 = (8, 1). 

af It is easy to see this fact from the definition l|3.12|l . Suppose that a module M is 0-semistable 
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Adding the extra node allows us to compactify the moduli space. In the lan- 
guage of D-branes, this corresponds to adding a single D6 brane filling the entire 
Calabi-Yau manifold, which is necessary to interpret the whole system as a six- 
dimensional U(l) gauge theory related to the Donaldson-Thomas theory. As we 
mentioned in Section r3.2.4[ the D6 brane serves as a flavor brane and adds an extra 
node exactly in the way described in the above paragraph. Note that, in the above 
paragraph, the ideal JF is generated by the original ideal T . In the quiver gauge 
theory, this corresponds to the fact that the flavor brane does not introduce a new 
gauge invariant operator to modify the superpotential. In this way, we arrive at 
the definition of non-commutative Donaldson-Thomas invariant as the Eulcr char- 
acteristic x(A4^(A)) of cohomologies El- With our identification of Mf Q (A) with 
the moduli space of solutions to the F-term and D-term conditions, x(-^i^(^)) 
computes the Witten index of bound states of DO and D2 branes bound on a single 
D6 brane ignoring the trivial degrees of freedom corresponding the center of mass 
of D-branes in R 1 ' 3 (see Section [2~lT) . 



We have chosen a specific vertex ig to define the non-commutative Donaldson- 
Thomas invariant. The %q dependence drops out in simple cases such as C 3 and 
conifold, but in general x{-M^(A)) depends on the choice of the We note that 
the quiver gauge theory discussed in Section l3~2l also has an apparent dependence 
on ig. 

The most systematic method for computing the number and direction of arrows 
in the quiver diagram is to use fractional branes, or more mathematically exceptional 
collections EB- The resolved conifold X is a normal bundle Oc( — 1) © Og?(~l) over 
a curve C ~ P 1 , and the fractional branes are given by Oc and Oc{— 1)[1] 
where [I] means the shift of the complex by onelfj. These sheaves correspond to the 
nodes of the quiver diagram, and we have 4 (=2+2) arrows since 

Ext^Oc, O c (-l)[l]) Ext 1 (O c (-I)[I],O c ) ~ C 2 . (3.15) 

In this language, the D6-brane is given by the structure sheaf Ox [1] of the whole 
X, with a suitable shift (in this case, I). By computing the Ext groups, we can 
determine the quiver diagram including the D6-branc node and the arrows starting 



but not 0-stable. This means that M has a submodule M' such that 

^2 §i(dimeiM') = 0. (3.14) 
i£Qo 

This is impossible, however, since the left hand side is dime*M' = 1. 

ag More precisely, this invariant should be defined from the Euler class of the obstruction bundle 
over the moduli space. The exis tence of obstruction theory guarantees the equivalence of these two 
quantities up to a possible sign 1321121] 

ah Exceptional collection appear in the work of Ref. 11321 and in Ref. 11331 See Refs. IT341 11351 

al Recall that fractional branes are in the derived categories of coherent sheaves, which consists of 

complexes of coherent sheaves (localized with respect to the quasi- isomorphism). 
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from it EH 

As a digression, this understanding gives a bonus. As is well-known, the choice 
of the fractional branes is far from unique. For example, instead of Ox[l]? Cc an d 
Oc(— we can choose Ox [1], Oc(— n) and Oc{—n — 1)[1], where n is a positive 
integer. Computing the Ext group of the quiver diagram, we have a different quiver, 
shown in Fig.[8j 




Fig. 8. The quiver with a different choice of fractional brancs, for n = 2. In general, we have 
n + 1 (n) arrows starting from (ending at) the D6-brane node, corresponding to the chamber C n 
discussed later in Section r5.4.1l 

Physically, the quivers are related by Seiberg dualities (called mutations in math- 
ematics 0)- Note that in the construction given below, different quivers (and dif- 
ferent superpotential) give a different path algebra, and thus a different crystal. 
This means that the BPS index changes, which is reminiscent of the wall crossing 
phenomena. Indeed, this is not an analogy; as clarified in Refs.l86 l l54 ( in the quiver 
language wall crossing phenomena can be understood as Seiberg dualities on the 
quiver. We will return to the issue of wall crossing in part II. But before that, let 
us describe the structure of the crystal. 



3.4. Crystal Melting 

In this section, we define a statistical mechanical model of crystal melting and show 
that the model reproduces the counting of BPS bound state of D-branes. Using 
the quiver diagram and the superpotential of the gauge theory, we define a natural 
crystalline structure in three dimensions. We specify a rule to remove atoms from 
the crystal and show that each molten crystal corresponds to a particular BPS 
bound state of D-branes. We use the result of Ref. [34] to show that all the relevant 
BPS states are counted in this way. 



aJ Computation of the superpotential by this method is in general difficult, although possible in 
some examples. See Ref. 11361 

ak T here is a conjecture that Seiberg duality is the same as the mutation of exceptional collections 
[1341 
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3.4.1. Crystalline Structure 

Mathematically, the three-dimensional crystal we define here is equivalent to a set 
of basis for Aei , where A is the factor algebra A — CQ/J- of the path algebra CQ 
divided by the ideal J- generated by the F-term constraints and ei is the path of 
zero length at the reference node i , which is also the projection operator to the 
space of paths starting at io- Colloquially, the crystal is a set of paths starting at io 
modulo the F-term relations. As we shall later, it corresponds to a single D6 brane 
with no DO and D2 charges. We interpret Aei in terms of a three-dimensional 
crystal as follows. 

The crystal is composed of atoms piled up on nodes in the universal covering 
Q on R 2 . By using the projection, tt : Q — > Q, each atom is assigned with a color 
corresponding to the node in the original quiver diagram Q. The arrows of the 
quiver diagram determines the chemical bond between atoms. We start by putting 
one atom on the top of the reference node io- Next attach an atom at an adjacent 
node j € Qo that is connected to io by an arrow going from io to j. The atoms at 
such nodes are placed lower than the atom at io- In the next step, start with the 
atoms we just placed and follow arrows emanating from them to attach more atoms 
at the heads of the arrows. 

As we repeat this procedure, we may return back to a node where an atom 
is already placed. In such a case, we use the following rule. As we explained in 
Section r3.3.21 modulo F-term constraints, any oriented path a starting at io and 
ending at j can be expressed as Vi 0> jW n , where ui is the loop around a face in the 
quiver diagram and Vi j is one of the shortest paths from io to j (Fig. [9] for an 
example). This defines an integer h{a) = n for each path a. The rule of placing 
atoms is that, if a path a takes io to j and if h(a) = n, we place an atom at the 
n-th place under the first atom on the node j fFig.lTO)). If there is already an atom 
at the n-th place, we do not place a new atom. 




Fig. 9. This example shows that the green path is F-term equivalent to the shortest path (blue 
path) plus a loop. 

By repeating this procedure, we continue to attach atoms and construct a pyra- 
mid consisting of infinitely many atoms. Since atoms are placed following paths 
from io modulo the F-term relations, it is clear that atoms in the crystal are in 
one-to-one correspondence with basis elements of Aei . Note that by construction 
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Fig. 10. Structure of the crystal, representing an atom a at position j and depth n. Different 
paths from the top atom to the atom a (for example, the green path and the green path) are all 
F-term equivalent. 

the crystal has a single peak at the reference node iq. 

This defines a crystalline structure for an arbitrary toric Calabi-Yau manifold. In 
particular, it reproduces the crystal for C 3 discussed in Refs. [551 1751 and the one for 
conifold in Rcf.[33J See Fig.imfFig. [T2|) for the crystalline structure corresponding to 
the resolved conifold (Suspended Pinched Point singularity). In these two examples, 
the ridge of the crystal (shown as blue lines in Fig. [T2|) coincides with the (p, g)-web 
of the toric geometry. As we will discuss later, this is a general property of our 
crystal. 



Fig. 11. The crystal for the resolved conifold. The pyramid consists of infinite layers of atoms, 
each layer colored altcrnatingly by white and black. The four ridges of the pyramid are shown by 
four black lines extending to infinity. 



3.4.2. BPS State and Molten Crystal 

In the forthcoming discussions, the crystal defined above will be identified with a 
single D6 brane with no DO and D2 charges. Bound states with non-zero DO and 
D2 charges are obtained by removing atoms following the rule specified below. 
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Fig. 12. Starting from the universal cover Q of quiver Q shown on the left, we can construct 
a crystal on the right. Each atom carries a color corresponding to a node in Q, and they are 
connected by arrows in Q\. The green arrows represent arrows on the surfaces of the crystal, 
whereas the red ones are not. In the case of the Suspended Pinched Point singularity, the atoms 
come with 3 colors (white, black and gray), corresponding to the 3 nodes of the original quiver 
diagram Q on T 2 shown in Fig. [4] 

In Ref. EH [34j the Donaldson-Thomas invariants x(-M^(A)) are computed by 
using the U(l) 2 symmetry of the moduli space M.f Q corresponding to the transla- 
tional invariance of T 2 . By the standard localization techniques, the Euler number 
can be evaluated at the fixed point set of the moduli space under the symmetry 
Correspondingly, in the gauge theory side, BPS states counted by the index are 
those that are invariant under the global U(l) 2 symmetry acting on bifundamcntal 
fields preserving the F-tcrm constraints since those do not have extra zero modes 
and do not contribute to the index. We are interested in counting such BPS states. 

In order for a molten crystal to correspond to t/(l) 2 invariant 0-stable A- 
modules, we need to impose the following rule to remove atoms from the crystal. 
Let n be a finite set of atoms to be removed from the crystal. 

The Melting Rule: 

If aa e £1 for some a E A, then a S fi. (3.16) 

Since atoms of the crystal correspond to elements of Aei Q , we used the natural 
action of A on Aei Q to define aa in the above. This means that crystal melting 
starts at the peak at io and takes place following paths in Aei . An example of a 
molten crystal satisfying this condition is shown in Fig.[T3l 

The melting rule means that a complement X of the vector space spanned by £1 
in Aei gives an ideal of A. To see this, we just need to take the contraposition of 
the melting rule. It states: For any (3 £ I and for any a 6 A, af3 is also in I. 

Generally speaking, an ideal of an algebra defines a module. To see this, consider 
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Fig. 13. Example of a molten crystal and its complement f2. In this example Q contains 12 atoms, 
one hidden behind an atom on the reference point represented by a blue point. It is easy to check 
that f! satisfies the melting rule mentioned in the text. 

a vector \I) which is annihilated by all elements of the ideal I. From \I) , we can 
generate a finite dimensional representation of the algebra A by acting with elements 
of A on it. However, the converse is not always true. Fortunately, when modules 
are ^-stable and invariant under the U(l) 2 symmetry, it was shown in Refs. [33l [34l 
that there is a one-to-one correspondence between ideals and modules. It follows 
that our molten crystal configurations are also in one-to-one correspondence with 
A-modules and therefore with relevant BPS bound states of D-branes. This proves 
that the statistical model of crystal melting computes the index of D-brane bound 
states. 

It would be instructive to understand explicitly how each molten crystal config- 
uration corresponds to a BPS bound state. Starting from a molten crystal specified 
by £7, prepare a one-dimensional vector space V a with basis vector e a for each atom 
a G 0. For each arrow a of Q, define the action of a on V a by a(e a ) = ep when 
the arrow a begins from a and ends an another atom /3 G fi. Otherwise a(e Q ) is 
defined to be zero. Since an arbitrary path is generated by concatenation of arrows, 
we have defined an action of a G A on each V a . By linearly extending the action of 
a onto the total space M = (B a enV a , we obtain a A-module M. 

There are several special properties concerning this module M. First, the F-term 
relations are automatically satisfied. This is because when there exist two different 
paths a, b G A starting at a and ending at /3, a(e Q ) and b(e a ) are both defined 
to be ep. Second, by construction M is generated by the action of the algebra A 
on a single element ei a G Vi Q . In such a case M is called a cyclic ^-module, and 
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by lemma 2.3 of Ref. [Ml is also ^-stable. Third, by the cyclicity of the module it 
follows that M is U{1) 2 invariant up to gauge transformations. Therefore, M is a 
C/(l) 2 invariant (9-stablc module. It follows from the result of Ref. [M] discussed at 
the beginning of this section that M indeed corresponds a bound state of D-branes 
contributing to the Witten index. 

At the beginning of this subsection, we stated without explanation that the 
original crystal corresponds to a single D6 brane with no DO and D2 charges, and 
removing atoms correspond to adding the D-brane charges. To understand this 
statement, let us recall that, in Scction l3~2l we started with a configuration of DO 
and D2 branes on the toric Calabi-Yau manifold and took a T-duality along the fiber 
to arrive at the brane configuration. Thus, the number of D2 branes at each node 
j of the quiver diagram Q is a combination of DO and D2 charges before T-duality. 
It is this number that is equal to the rank of the gauge group at j. 

By using the projection it : Q Q — > Qo, the A-module M is decomposed as 
M = @j£Q Mj as we saw in Section [3~3l where 

Mj = V a . (3.17) 

CK£fi,7r(a)=,7 

In particular, the formula (|3.17|) means that the rank of the gauge group Nj = 
dim Mj at the node j is equal to the number of atoms with the color corresponding 
to the node j that have been removed from the crystal. Thus, removing an atom at 
the node j is equivalent to adding DO and D2 charges carried by the node j. It is 
interesting to note that each atom in the crystal does not correspond to a single DO 
brane or a single D2 brane, but each of them carries a specific combination of DO 
and D2 charges. In the crystal melting picture, fundamental constituents are not DO 
and D 2 branes but the atoms. This reminds us of the quark model of Gcll-Mann and 
Zwcig 1 1 where the fundamental constituents carry combinations of quantum 
numbers of hadrons, as opposed to the Sakata model! I, where existing elementary 
particles such as the proton, neutron and A particle are chosen as fundamental 
constituents. 

3.4.3. Observations on the Crystal Melting Model 

We would like to make a few observations on the statistical model of crystal melting 
that counts the number of BPS bound states of D-brancs. 

We have studied several examples of toric Calabi-Yau manifolds and found that 
the crystal structure in each case matches with the toric diagram. In particular, 
the ridges of the crystal, when projected onto the R 2 plane, line up with the (p, q) 
web of the maximally degenerate toric diagram. This phenomenon is discussed in 
|Appcndix A.l| There, we also explain the correspondence between molten crys- 
tal configurations and perfect matching of the bipartite graph introduced in Sec- 
tion l3X2l 

In the last subsection, we found it useful to describe BPS bound states using 
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ideals of the algebra A. In the case when the toric Calabi-Yau manifold is C 3 , ideals 
are closely related to the quantization of the toric structure as discussed in Ref. [73] 
The gauge theory for C 3 is the dimensional reduction of the Af = 4 supersymmetric 
Yang-Mills theory in four dimensions down to one dimension, and the bifundamcntal 
fields are three adjoint fields. The F-term and D-term conditions require that they 
all commute with each other. Thus chiral ring is generated by three elements x, y, z 
which commute with each other without any further relation. In this case, any ideal 
is characterized by the three-dimensional Young diagram ■k. Locate each box in 
the 3d Young diagram 7r by the Cartesian coordinates (i, j, k) (i, j,k = 1, 2, 3, ...) of 
the corner of the box most distant from the origin, and define J7 W to be a set of the 
3d Cartesian coordinates k) for boxes in ir. We can then define the ideal X v of 
the chiral ring by, 

X, = {x i - 1 y>- 1 w k - 1 \(i,j, k) £ fij. (3.18) 

In Ref. [751 this description was obtained by quantizing the toric geometry by using 
its canonical Kahler form and by identifying x t ~ 1 y : '~ 1 'w k ~ 1 as states in the Hilbcrt 
space. 

This can be generalized to an arbitrary toric Calabi-Yau manifold as follows. 
One starts with the quiver diagram corresponding to and use the brane tiling 
to identify the F-term equations. This gives the chiral ring generated by bifunda- 
mental fields obeying the F-term and D-term relations. As we saw in Section [3.4.31 
each BPS bound state is related to an ideal of the chiral algebra. We expect that 
such ideals arise from quantization of the toric structure. BPS bound states of D- 
branes emerging from the quantization of background geometry is reminiscent of the 
bub bling AdS space of Ref. 11401 and Mathur's conjecture on black hole microstatcs 

am 

3.5. The BPS Partition Function 

We have seen that a molten crystal is in one-to-one correspondence with a BPS 
state contributing to the index. It automatically follows that the BPS partition 
function, which is defined in Scction ["2.1l bv summing over BPS states contributing 
to the index, is the same as the statistical partition function of the crystal melting 
model. 

From the explanation around (|3.17p . the partition function expressed in terms 
of the crystal is given by 

Zbps(<Zo, • • • , On) = Yl ±C ( %P (Q) • ■ • lT m , (3-19) 
a 

where iu,(0), which is the same as dimMj, is the number of atoms in 57 with the 
i-th color. The parameters qi, counting the rank of the dimension of the vector space 
at vertex i, is a combination of D0/D2-brane chemical potentials. For example, in 
the conifold example dimension vectors are related to the D-brane charges rio,ri2 
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by (see the discussion at the end of Section l3.2.1l) 

dimMi = no, dim M 2 = n + n 2 . (3.20) 

This means that the relation between parameters are given by 

= ?o?i, Q = qi- (3.21) 

In the discussion so far we have not discussed the sign in (|3.19j) . which should 
be properly taken into account when we discuss the BPS index. Namely the correct 
statement is that the BPS index is the same as the number of configuration of 
molten crystals times a sign, which mathematically comes from the dimension of 
the tangent space of the vacuum moduli of the quantum mechanics, as mentioned 
in footnote[p]in Section l3.1.3l Recall that we are actually computing the BPS index, 
which is defined as the number of boson and that of fcrmions, and therefore can 
take negative as well as positive values. In general, for a configuration of molten 
crystal f2, Theorem 7.1 of Ref. |M]says that the sign is given by 

(_l) w io( n )+(«'( n )' w ( n )) j (3.22) 
where (a,/3) for a — (ai)i 6 g ,/3 = (/3i)i g Q is the Ringel form of the quiver, given 

by 

(a lJ 8>= ( 3 - 23 > 

ieQo aeQi:i-yj 

In general we have to insert this sign to each term of the crystal melting partition 
function separately. When the toric Calabi-Yau manifold X does not have a compact 
4-cycle, however, these signs are simply taken care of by the change of the signs of 
the chemical potentials in the partition function. To show this, we begin with the 
fact that in these examples the quiver diagram takes the following form 

Suppose that the quiver has TV nodes, together with one node corresponding to 
the D6-brane. Then we can number the N nodes by i = 1 . . . , N such that (a) there 
is always exactly one arrow from node i to node i + 1, (b) there is always exactly 
one arrow from node i to node i — 1, and (c) there may or may not be an arrow from 
node i to node i. In this situation contributions from arrows of type (a), (b) cancel 
out in the second term of (|3.22|) . and the only remaining contribution is from (c): 

= w io (n) +J2 n i - J2 Qi ( mod2 )- 

i i-3a£Qi ,a:i^ri 

Thus the total sign can also be written as 

(-l)£«es n < (3.24) 

al This statement is well-known in the brane tiling literature. See for example Fig. 6 of Ref. 1 1081 for 
the basic building blocks of the bipartite graph (fivebranc diagram), from which we can construct 
the bipartite graph. The dual graph gives a quiver diagram. 
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for a set S = {io} U {i\ 3 a G Qi, a : i — > i}. This means that the sign is conveniently 
taken care of by the replacement |ffj (we use pi for the parameters in the generating 
functions of the crystal without signs) 

q, -*Pi = -Qi (if i€S), qi^Pi = qi (if i f S). (3.25) 

In other words, the precise relation between the BPS partition function (|3.19[) and 
the crystal partition function (defined by (|3.19p without signs) is given by 

Zbps(?o, . . . , gjv) — ^crystal (po,...,pn) (3.26) 

under the identification in (|3.25|) . 

As an example, for the resolved conifold there is no arrow of type (c), and we 
have {S} = {1} and 

Po = qo, Pi = ~qi- (3.27) 
This concludes the discussion of signs in the BPS partition function Q3.19p . 

Since we already know the explicit form of the crystal and the sign, we can 
compute Zbps from lowest order terms. For example, for the conifold example we 
have (see Fig. [T4|) 

^BPs(?0,gi) = 1 + <7o - 2q <7i - 4<7o<ji + q Q qf + . . . , (3.28) 

which becomes, under the parameter identification (|3 . 2 1 1) . 

^bps(<Z, Q) = 1 + qQ- 1 -2q- Aq 2 Q^ + qQ + . . . , (3.29) 

We can perform similar computations for any crystal, however complicated it is. 

In some examples we can write down a closed expression for the BPS partition 
function, inste ad o f computing each term separately. For example, for the resolved 
conifold Young EH has proven combinatorially that 

^crystal (PO, Pi) = M ( PoPl ) 2 TJ(1 + (p Pl)>l)" + ipoPlT P^T ■ (3.30) 

n n 

This means, using (f3"T2T|) and (f3~2"Tj) . 

^bps(<7, Q) = M{-qf - {-q) n QT - (-qTQ^Ti (3.3i) 

n n 

where M (q) is the MacMahon function 

M(q)= TJ(1- ST" (3-32) 

n>0 

This is very similar to the Gopakumar- Vafa^^^ expansion of the topological string 
partition function 

Zto P (q, Q) = M(-q) JJ(1 - {-q) n Q) n , (3.33) 

n 

am As emphasized previously, this is not the case in general, e.g. for the canonical bundle of F 1 xP 1 . 
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Fig. 14. The examples of the states of the conifold crystal, starting with the ground state. 

where the parameters q and Q in the above expression are related to the topological 
string coupling constant g s and the Kahler moduli t bvF"l 

q=-e' 9 % Q = e~ t . (3.34) 

In fact, we can write 

Zbps(<Z, Q) = Z top {q, Q)Z top {q, Q^ 1 ). (3.35) 

Similar expressions hold in all the toric Calabi-Yau 3-folds without compact 4-cycles. 
To explain the appearance of the topological string theory, as well as the relation 

an For notational simplicity, in the following sections we will be sloppy and often absorb the minus 
sign in 113,340 into the definition of q to write q = e~ 9s . 
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(I3.35[) . will be the topic of Scction[SJ 
3.6. Summary 

In this section, we established the connection between the counting of BPS bound 
states of DO and D2 branes on a single D6 brane to the non-commutative Donaldson- 
Thomas theory. We studied the moduli space of solutions to the F-term and D-term 
constraints of the quiver gauge theory which arises as the low energy limit of the 
brane configuration. We found the direct correspondence between the gauge theory 
moduli space and the space of modules of the factor algebra of the path algebra 
for the quiver diagram quotiented by its ideal related to the F-term constraints, 
subject to a stability condition to enforce the D-term constraints. Using this cor- 
respondence, we found a new description of BPS bound states of the D-branes in 
terms of the statistical model of crystal melting. The crystalline structure is de- 
termined by the quiver diagram and the brane tiling which characterize the low 
energy effective theory of D-branes. The crystal is composed of atoms of different 
colors, each of which corresponds to a node of the quiver diagram, and the chemical 
bond is dictated by the arrows of the quiver diagram. BPS states are constructed 
by removing atoms from the crystal. 

We naturally encounter the wall crossing phenomena, the theme of the sec- 
ond half of this paper, by studying the relation between the commutative and 
non-commutative Donaldson-Thomas invariants. The degeneracy of D-brane bound 
states changes when the value of 6, used to define the stability condition, which is 
the wall crossing phenomena discussed in Section l2T2lF°l. T he jump in the degen- 
eracy can be computed by the wall crossing formulaSZES], and if we start from a 
particular chamber and apply the wall crossing formula, we can obtain the value of 
xi-Mf {A)) in any chamber we want. We will see this explicitly in the conifold exam- 
ple in Section[6] As will be discussed there, wall crossing relates non- commutative 
Donaldson-Thomas invariants to commutative Donaldso n- Th omas invariants and 
to new invariants defined by Pandharipande and Thomas This story is further 
generalized by Ref. [52] when A has no internal lattice point. 

It has been proven recently that the topological string theory is equivalent to 
the co mmutative Donaldson-Thomas theory for a general toric Calabi-Yau manifold 
| 75 | 79 | g mce t ne commutative Donaldson- Thomas theory count BPS states for some 
choice of stability condition, 

■Zbh = Zto-p, (3.36) 

is indeed true in some chamber. On the other hand, our result shows that the 
relation, 

■^BH -^crystal melting; 
ao Thc FI parameters is a function of the Kahler moduli. 



(3.37) 
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holds in another chamber, where Zg H is the BPS state counting for another choice 
of the stability condition. Combining these two results, we find that the topological 
string theory and the statistical model of crystal melting are related by the wall 
crossing, and we have 

^crystal melting ~ ^to P (modulo wall crossings) . (3.38) 

Since there is no wall crossing phenomenon between the commutative and non- 
commutative Donaldson-Thomas theory on C 3 , this result does not contradict with 
Ref. [7SJ where a direct identification of the topological string theory and the crystal 
melting is made for C 3 . In general, we expect that a proper understanding of the 
relation between the topological string theory and the crystal melting requires that 
we take the wall crossing phenomena into account. This motivates the study of wall 
crossing phenomena in part II. 



4. Thermodynamic Limit 

In this section we study the thermodynamic limit of the crystal melting model 
defined in the previous section. This section is based on the work^S. 



4.1. Overview: Why Thermodynamic Limit? 

In the previous section we have seen that the BPS counting problem is beautifully 
solved by the crystal melting model, making an explicit computation possible. How- 
ever, the question still remains whether the crystal expansion is just a technical tool 
to obtain the answer, or whether it has much deeper significance. 

There is one notable feature of the crystal melting model: it is a strong coupling 
expansion with respect to the (topological) string coupling constant. In the usual 
string theory, we consider a perturbative, weak coupling (g s <C 1) expansion 

Z top =exp [f^Q^-^gj , (4.1) 

whereas the crystal expansion gives 

^crystal = ^ P)^^ , (4.2) 

n,f> 

where we used the parameter identification (|3.34p , which holds for any toric Calabi- 
Yau manifold. Due to the exponential factor, this expansion is good when g s ^> 1. 

This in particular means that crystal melting provides a tool to study the physics 
at the Planck scale When g s is small, as is usually the case, the string scale l s is 

ap The following explanation is a bit imprecise in that g s here is the topological string coupling, 
and is different from the physical string coupling constant which really determines the relation 
between the string scale and the Planck scale. However, the topolog ical string theory is relevant 
for counting of microstates of black holes in the superstring theory ED and we expect that our 
result sheds some light on quantum nature of spacctime in the superstring theory also. 
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larger than the Planck scale l p . This means that we see stringy corrections first, and 
quantum corrections are obscured. When g s is large, however, l s is smaller than the 
Planck scale l p kicks in first. This means that we can explore the 'quantum' aspects 
of geometry. 

The natural question is to understand the relation between the two descriptions. 
In this section, following the result of Rcf. [35] we will show how the smooth Calabi- 
Yau geometry emerges from the discrete structure of the crystal melting model in 
the limit g s — ¥ 0. Since g s is the chemical potential for the total number of atoms 
removed, g s — > means that we are in the thermodynamic limit (high temperature 
limit), where an infinite number of atoms are removed. It is reasonable to expect 
that a classical geometric picture emerges in the limit of large DO and D2 charges 
since it can represent a large black hole in the superstring theory. The corresponding 
thermodynamic limit in the dimer model was studied in Ref. 11431 

There are two main results proved in this section. First, the limit shape (see 
Fig. ll5[) of the crystal in the thermodynamic limit coincides with a projection 
(amoeba) of the mirror Calabi-Yau manifold. 




Fig. 15. Limit shape of the crystal, in the C 3 example. In the thermodynamic limit, the boundary 
of the molten region (where the atoms are taken away) becomes a smooth curve, which is called 
the limit shape. The claim is that this curve is the same as the amoeba of the mirror curve, which 
is shown in Fig. 1171 (b). 



This is proven as follows. A toric Calabi-Yau 3-fold X is a Kahler quotient of 
C£+f_by U(1) F , and its mirror manifold X is defined by the polynomial equation 



iL 1 "by 
114411321 



uv + P(z, w) = 0, ii.ueC, z,wGC x . 



(4.3) 
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Here P(z, w) is a Newton polynomial of the form, 



F+3 



P(z, w) = Ci(t)z 



w 



(4.4) 



and Cj(t)'s are functions of the Kahler moduli t of the original toric 3-fold X. 
The exponents (nj,m,) € Z 2 correspond to lattice points of the toric diagram. For 
example, for the mirror of C 3 , P(z, w) is given by 



and for the mirror of the resolved conifold (0{— 1) © O(-l) bundle over P 1 ) by 



In this section, we will show that the Newton polynomial P(z, w) for the mirror 
of a Calabi-Yau manifold is identical to the characteristic polynomial of the corre- 
sponding dirtier model, which is the partition function of the model on a torus. The 
relation between P(z, w) and the characteristic polynomial had been discussed ear- 
lier in Ref. 11041 Here, we will prove their precise equality including the dependence 
on the moduli t. As discussed in |Appcndix A.l| the dimcr model is equivalent to 
the crystal melting model discussed in the previous section. Since we can show from 
the result of Ref. 11431 that the amoeba of the characteristic polynomial gives the 
limit shape of the crystal, this proves our result. 

As the second important result in this section, we will show that the partition 
function of the dimer model evaluated in the thermodynamic limit is equal to the 
genus-0 limit of the partition function of the topological string theory on A. As 
discussed in Scction ["3.6[ the dimcr model has b een formulated to describe the non- 

IOOIO/1 1 

commutative Donaldson-Thomas theory 1 1 while the topological string theory 
for a general to ric C alabi-Yau manifold is equivalent to the commutative Donaldson- 
Thomas theory as shown by Ref. [79] Despite the non-trivial wall crossing factors 
between the two theories discussed in Scction l3~6l the thermodynamic limit of the 
crystal partition function gives the topological string partition function. 

Our results suggests that the atoms in the crystal are a discretization of the 
Calabi-Yau geometry at the Planck scale, and the thermodynamic limit is a process 
where a smooth geometry emerges from a collection of such discretized piece of the 
geometry. The emergence of Calabi-Yau geometry from the thermodynamic limit 
has been observed in Ref. [55] in the case of C 3 , and is interpreted in Ref. [75] from 
the viewpoint of Kahler gravity Ref. 145. In this section, we make the connection 
sharper and more explicit by showing the direct connection between the partition 
functions of the crystal melting model and the topological theory for a general toric 
Calabi-Yau 3-fold. 



P(z, w) = 1 + z + w, 



(4.5) 



P(z, w) = 1 + z + w + e t zw. 



(4.6) 
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4.2. Thermodynamic Limit of the Crystal Melting Model 

The main object of study in this section is the partition function, 

Z=J2 V{n, p^e- 3 ' 71 -*' 131 , (4.7) 

n,Pi 

where 0(n, 0) is the Witten index for bound states of n DO branes and f3 = (/?/) 
D2 branes on the i-th 2-cycle (/ = 1, dimi?2(A/)) with a single D6 brane on a 
toric Calabi-Yau manifold. According to the dictionary in the previous section, the 
Witten index is equal (up to a sign) to the number of molten crystal configurations 
where n is the total number of atoms removed, whereas the relative numbers of 
different types of atoms removed from the crystal are specified by /3/'s. The variables 
g s and t 1 arc chemical potentials for DO and D2 charges, and as we have seen in 
p. 21 [) El we w iU identify g s as the topological string coupling constant and t 1 as 
the Kahler moduli of the toric Calabi-Yau manifold X . 

The behavior of Z as g s — > can be evaluated by the result of Ref. 11431 Consider 
a finite covering the original quiver diagram, N times in one direction and N times 
in another direction on T 2 . N is introduced as an infrared regulator, and we will 
take N — > oo at the end of the computation so that we have the dimer model on the 
plane R 2 . The surface of the crystal is determined by the height function h over the 
plane. To define h, we start with the canonical perfect matching mo of the dimer 
model corresponding to the initial crystal configuration with no DO or D2 charges 
(in this section, we heavily use the language of dimer models, whose equivalence 
with crystal melting is explained in | Appendix A.l ). For any other perfect matching 



m, the superposition of mo and m gives a set of closed loops on the dimer graph. 
If m corresponds to a bound state with finite DO and D2 charges, m and mo differ 
only in a finite region on the graph. The value of the height function h at i-th node 
of the N x N cover is defined so that it is far away from the region where m 
and mo differ, and it increases by 1 every time we cross a closed loop as we move 
inside of the region (sec Fig.|37]for an example). The corresponding molten crystal 
configuration is obtained by removing h(i) atoms of the initial crystal over the node 
i. In particular, 



(4.8) 



where n is (as before) the DO-brane charge, or equivalently the total number of 
atoms removed. 

To take the thermodynamic limit, it is useful to introduce the Cartesian coordi- 
nates (x, y), < x,y < 1, on the N x N covering of T 2 . In the limit where g s <C 1 
and 1 C JV, the height h(i) becomes a smooth function h(x,y). We rescale the 



aq For simplicity we here neglect the sign in 13.211 . 
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height function by the factor of 1 /N so that 

f 1 n 
N 2 dxdy h(x,y) = -, (4.9) 

to take into account the large N scaling of the partition function discussed in 
Ref. 11461 and quoted as Theorem 2.1 in Ref. H43l H. The statistical weight in the 
thermodynamic limit is given by an integral of a surface tension a(dh), which is a 
function of the gradient of h, asQ 



exp 



»i 

r2 



N max/! / dxdy [~a(dh) — g s N h(x, y)] 



n 



(4.10) 



The integral of g s N h(x,y) in the exponent comes from the weight factor e~ 9s " 
in (|4.2[) . and we used (|4.9[) . Note that interestingly the combination g s N appears, 
which is the 't Hooft parameter. In the thermodynamic limit, we look for a height 
function h(x,y) which maximize the exponent. 

In order to derive the macroscopic surface tension from the microscopic crystal 
melting model, we need a few definitions. We first define a characteristic polynomial 
P(z, w) as the partition function per fundamental domain, i.e. the s um o f perfect 
matchings with weights assigned to edges of the dimcr model on T 2 11421. The pa- 
rameters z and w specify the height difference in the language of dimcr models. 

Let us explain the meaning of parameters z and w in more detail, using the 
Suspended Pinched Point shown in Fig.Ql)] as an example. In this example, the 
bipartite graph has 7 edges and 6 perfect matchings. When we denote the weight 
for the i-th edge by e Wi , the characteristic polynomial is given by 

This is a function of seven variables, but there is a redundancy in the description. 
We call this redundancy "gauge symmetry" . For example, we can simultaneously 
multiply by a constant c all the weights for edges incident to the vertex a (i.e. 
wi, u>4, 1U5, wq), and the partition function stays the same up to the overall multi- 
plication; this is because every perfect matching has exactly one edge incident to 
vertex a. This means that we can use this gauge degrees of freedom to set W5 = 0. If 
we repeat the similar procedure for the other vertices, we can also set W3 = wj = 0, 
and we have fixed all gauge degrees of freedom. The characteristic polynomial now 
simplifies to 

1 + e W2 + e Wi + e W2+Wi + e W6 + e Wl . (4.12) 

Instead to introducing redundant variables with gauge symmetry as above, we 
can be more economical and express this result in terms of gauge invariant variables. 

ar This statement is consistent with our intuition that the height diverges in the same order as the 
size of the enlarged torus. 

as Here we only show the g s dependence explicitly and the dependence on the Kahler moduli t 1 is 
in a. 



42 Masahito Yamazaki 



A 5 













3 



5 

bipartite graph 



/ 












5, 


7 


/ 

















=S<. 




-A 





2,5 




2,4 



3, 6 

perfect match ings 











-A 


ff 


4 




/ 










PI 


i, 


3 



Fig. 16. The 6 perfect matchings for Suspended Pinched Point bipartite graph (Fig.[5]l. The 
numbers assigned to each perfect matching shows which edge is included in the matching. 



The gauge invariant obscrvablcs are given by the "mag nctic flux" EH. To define 
this, let us fix an oriented loop 7 = {wi,bi, ■ ■ ■ , Wk-i, bk-x, Wk}, = Wq, where 
Wi (hi) are white (black) vertices. We then define the magnetic flux along 7 to be 



1) 



(4.13) 



where W( Wi ^ is the weight assigned to an edge (wi,bi) connecting a white vertex 
u>i and a black vertex hi. This is clearly gauge invariant. In general, when we denote 
the number of faces of the bipartite graph by F (this is the same as the number 
of nodes of the quiver diagram), we have F + 1 independent gauge invariants: 2 of 
them comes from a and /3-cycles of the torus, and F — 1 comes from the flux around 
each of the F faces, with one condition that the sum of the fluxes of all faces sum 
up to zero. 

In the Suspended Pinched Point example, we have F + 1 = 4 independent 
parameters, and correspondingly we define 



z = e " , w = e 

— tl _ -W 3 ,+W l —W 5 +W 3 -W 2 +W7 



where the parameters z and w correspond to the fluxes along the a and /3-cycles, 
and t\ and to those along cycles around a face. In these variables the previous 
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equation (|4.12[) becomes (after setting = w$ = W7 = as before) 

P(z, w) = l + z + er il z + e^z 2 + e~ tx ~ H zw + e~ tl z 2 w. (4.14) 
Given P(z,w), we next define its Ronkin function 033 R(x 7 y) by 



2tt 



ded4> 



R(x,y) = lnP(e^,e^)^=. (4.15) 

The integrand in this definition diverges when P(e x+lB , e v+l ^) = 0, however the 
function itself is well-defined. To see this, let us compute the its derivative: 



d x R{x,y) 



2jv 







d z p{e x + ie , e y+^) x+i0 d6d(t> 

p( e x+ie >e y+i<f>) C (2tt) 2 



dz f( Z,W h 1 dzdw (4 16) 

\z\=e*,\w\=ev P{z,w) (2tt) 2 Z W ' 

This can be evaluated as a residue integral. For example, suppose that (x, y) is 
outside amoeba which is a subset of M 2 defined by l au l 

Amoeba = {(x, y) £ R 2 : P(e x+l9 , e y+l ' p ) = for some (0, 0)}. (4.18) 

Then the only pole contributing to (|4.16[) is located at infinity. When the lead- 
ing contribution of P(z,w) at infinity is given by z m w n , we can easily obtain 
{d x R{x 1 y),d v R(x,y)) = (m, n) . In other words, Ronkin function is linear in the 
complement of amoeba | 147 | 152 [ 

Let us here give the simplest example: C 3 , whose toric diagram is given in Fig. 1171 
(a). The characteristic polynomial is given by 

P(z, w) = z + w + 1, (4.19) 

and its amoeba is shown in Fig. [17] (b). The derivatives of the corresponding Ronkin 
function are computed to be 



for x < log \e y - 1 

/ e 2x_ £ 2y_ l \ 

d x R(x,y) = I 1 - - cos 1 ( giv J f ° r bg ~ ^ " X ~ bg ^ + ^ 

for x > log \ e y + 1 



at See Ref. 1 1481 for a survey of aspects of amoebae, and Rcfs. 149, 150 for their uses in other aspects 
of supersymmetric gauge theories. 

au There is a similar notion, called alga or coamoeba, which is defined by 

Alga= {(0,4>) S T 2 : P{e x+iB ,e y+i '*') = for some (x,y)}. (4.17) 

It is known that coamoebas are equivalent to fivebranc diagrams, and therefore yield the brane 
configurations and the bipartite graph (see Ref. 11021 for deta iled discussions). They also have 
interesting connections with mirror symmetry | 104 | 122 | 123 | 1 51] 



44 Masahito Yamazaki 



(O.li 



(0,0) (1,0) 




(a) toric diagram (b) amoeba 

Fig. 17. The toric diagram (a), and amoeba (b) for C 3 , adapted from Ref. 11501 



d y R(x, y) 



1 cos 

7T 



2e x 



for y < log le^ — 1 1 
for log \e x - 1| < y 

for y > log \e x + 1| 



2y 2x i ' 

1 for log \e x - 1| < y < log \e x + 1| 



and their plots are given in Fig. [TH1 Note that the derivatives of the Ronkin function 
takes a constant value at each complement of amoeba, as expected. 





(a) Ronkin function (b) gradient of Ronkin function 

Fig. 18. (a) Ronkin function and (b) the gradient flow of the Ronkin function, from Ref. 11501 



We now come back to the thermodynamic limit of the partition function (|4.10[) . 
The reason we introduced the Ronkin function of the characteristic polynomial 
is that the surface tension a(d x h, d y h) is the Legendre transform of the Ronkin 
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function with respect to (s, t) = (d x h, d y h) as (Theorem 3.6 of Rcf. 114=3)) 

R(x, y) = max s , t [-cr(s, t) + xs + yt] . (4.22) 

The first step in relating the dimer model to the topological string theory is to 
show that the characteristic polynomial P(z, w) of the dimer model is equal to the 
Newton polynomial P(z,w) (|4.4[) for the mirror Calabi-Yau manifold (|4.3[) . 

P{z,w) = P(z,w). (4.23) 

According to Ref. I112( there is a one-to-one correspondence between perfect 
matchings of the dimer model on T 2 and bi-fundamental fields of gauged linear 
sigma model appearing in the Kahlcr quotient construction of the toric Calabi-Yau 
manifold X. They are then related, by change of variables described in Ref. 11121 to 
lattice points of the toric diagram and to terms z ni w mi in the Newton polynomial 
(|4.4j) . This shows that there is a one-to-one correspondence between terms in P(z, w) 
and P(z, w), as pointed out by Ref. 11041 Furthermore we can show that their coef- 
ficients agree. The Kahlcr moduli of X are the Fayet-Iliopoulos (FI) parameters of 
the quiver quantum mechanics. According to the dictionary between quiver gauge 
theories and dimer models (Scction f3.2|) . FI parameters are associated with nodes of 
the quiver diagram, or equivalcntly to faces of the dimer model. Thus, we can iden- 
tify the FI parameters with the magnetic fluxes through faces of the dimer model 
(|4.13[) , which parametrize the energy of each perfect matching of the dimer model 
(recall the example of the Suspended Pinched Point discussed previously; there the 
magnetic fluxes along faces are denoted by t\ and £2)- Each perfect matching ap- 
pears in P{z, w) with the weight given by an exponential of the fluxes. On the other 
hand, the Newton polynomial is a sum of z ni w mi , each of which corresponds to a 
lattice p oint of the toric diagram and is weighted by an exponential of the Kahlcr 
moduli f E31 We have verified that the weight for perfect matchings and the weight 
for lattice points of the toric diagram agree, and this proves the identity (|4.23[) . 

Combining (|4.10[) and (|4. 22[) and discarding a term in total derivative in (x,y), 
which is justified by the subtraction of the linear piece in R{x, y) discussed in the 



av Let us given a rough outline of why the Ronkin function arises here. I|4,l0j l shows that in the 
thermodynamic limit the surface tension is related to 

lim -^logZjv, (4.20) 

where Zjv is defined for dimer model with fundamental domain enlarged by N X N. Now since the 
char acter istic polynomial can be written as a determinant of the Kasteleyn matrix EEH it follows 

that ma 

Z N {z,w)= Yl EI ^O^o). (4.21) 



Plugging 114.211 1 into 114.201 1. the product is turned into a Riemann sum, which is nothing but the 
integral in the definition of the Ronkin function in \ 1.1511 . 
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next paragraph, we find 
Z ~ 



exp 



By reseating (x, y) by the factor of g s N/2, this becomes 

—n [ dxdyR(x,y) 
.9s J 



TV / dxdyR 



exp 



g s N g s N 
-rrx,—-y 



(4.24) 



Note that the TV dependence has disappeared except that the range of the (x, y) 
integral has been rescaled by the factor of g s N/2. For TV — > oo with small but fixed 
g s , we have an integral over the whole (x, y) plane. 

The integral (|4.24|) in the large TV limit is divergent. To identify and subtract the 
divergent part, we use the amoeba of the characteristic polynomial defined in f|4. 18|) . 
In the ther modynamic limit, the amoeba corresponds to the liquid phase of the 
crystal ESQ If there are no interior points in the toric diagram, the complement 
of the amoeba is the solid phase, where the crystal retains its original shape. There, 
the Ronkin function R(x, y) is linear. If there are interior lattice points in the toric 
diagram, the amoeba acquires holes, inside of which are in the gas phase, where the 
Ronkin function is again linear but the slope of the crystal surface is different from 
the original one. The integral (|4.24j) becomes finite if we subtract the linear piece of 
the Ronkin function in the solid phase so that the partition function is normalized 
to be 1 for the initial crystal configuration. 



4.3. Topological String at Genus 

Our next task is to compute the genus-0 topological string partition function To 
of the toric Calabi-Yau manifold X and compare it with the thermodynamic limit 
of the partition function of the crystal melting model (|4.24j) . For this purpose, it is 
convenient to use the mirror Calabi-Yau manifold X defined by the equation (|4.3j) 
since J-"o can be evaluated by the classical period integral as, 

T Q = f n, (4.25) 

where f2 is the holomorphic 3-form on the mirror, and Cq is the Lagrangian 3-cycle 
which is the mirror of the 6-cyclc filling the entire toric Calabi-Yau manifold. 

A ccor ding to the microscopic derivation of the mirror symmetry by Hori and 
VafaESl the sigma-model on the toric Calabi-Yau manifold is equivalent to the 

"The paper Rcf. 143 gives three different phases of the dimer model: solid, liquid and gas. In the 
solid phase (frozen phase), the height difference between points arbitrary far away are deterministic. 
Non-solid phases are classified into either a gas phase (rough phase) or a liquid phase (smooth 
phase) depending on whether or not the height difference between the two point is bounded 
independently of the choice of the points. It was shown in Rcf. 143 that an unbounded exterior 
region, an interior region and a bounded exterior region of the amoeba correspond to a solid, liquid 
and gas region, respectively. 
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Landau-Ginzburg model with the supcrpotcntial 

W{u, x, y) = e u P(e x , e y ). (4.26) 

It was shown in Ref. 11321 that an integral of e in a 3-dimensional subspace of the 
(it, x, y) plane can be transformed into a period integral of the holomorphic 3-form 
f2 on the mirror Calabi-Yau manifold. Thus, we should be able to evaluate (|4.25[) 
as an integral of e w . To do so, we need to identify the contour of the integral. 

Since g s and t s in (|4.2p are taken to be real in the dimer model, the Newton 
polynomial P{x, y) in our case is with real coefficients. The mirror manifold (|4.3[) 
has the complex conjugation involution, and thus the fixed point set is a natural 
candidate for Cq. In fact, following the mirror symmetry transformation as described 
in Ref. 11441 we find that the 6-cycle of the original toric Calabi-Yau manifold X 
corresponds to the real section in the (u, x, y) space in the mirror X. Thus, we find 

poo /-OO /• 

Jo = / du dxdye eUp{e *< eV) = - / dxdy In P(e x ,e v ). (4.27) 

JQ J-oo J 

The divergent part of the integral in the (x, y)-plane can be removed by subtracting 
a linear term in \nP(e x , e y ) for x, y — > oo as we did for (|4.24p . The integral (|4.27|) is 
almost equal to the exponent of the partition function (|4.24j) of the crystal melting 
model, except that we do not have the averaging over the phases (9, (f>) to define 
the Ronkin function as in (|4.15|) . It turns out that the integral over (x,y) in (|4.27|) 
removes the dependence on (#, 4>), and thus the averaging process is not necessary. 

To see this, let us define a generalization of J-q for an integral of (x, y) with 
arbitrary phases as 

.F(M) = J dxdylnP(e x+lf> ,e v+zf> ). (4.28) 

Taking derivatives of T with respect to (9, 4>), we find the integrand becomes a total 
derivative in (x, y) as in 

4e + ^)^ = J dxdy 1 {4x + 4) lnP(eX+ " ^ 

If we choose (a, /?) so that it is not in the directions of the tentacles of the amoeba 
(|4.18[) . the boundary term is removed by the regularization and we find (adg + 
l3drf,)J- = 0. Since a, /? is arbitrary except in the directions of tentacles of amoeba, 
T is independent of (9, <fi) and agrees with its average. Namely, 

J"o = - J dxdy In P(x, y) = - J dxdy R(x, y). (4.29) 

Thus, we found that the thermodynamic limit of the partition function of the 
crystal melting model given by (|4.24j) is equal to exp(— -yJ-o), which is the genus-0 
partition function of the topological string theory. This is what we wanted to show. 

The fact that that crystal melting gives the topological string theory in the 
thermodynamic limit, despite the non-trivial wall crossing phenomena between the 
two, poses a puzzle. We will see later in Section[5] how this puzzle is solved in the 
analysis of the wall crossing phenomena. 
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5. Wall Crossing and M-theory 

In section [3] we succeeded in computing the BPS partition function as a statistical 
partition function of crystal melting. However, the crystal melting gives only an 
expansion of the BPS partition function, and whether the partition function can be 
expressed in a closed expression is a different issue. In Scction ["3.5[ we have seen that 
in the example of the resolved conifold the BPS partition function takes a simple 
expression, which is essentially a square of the topological string partition function. 
The question we would like to address in this section is to understand this curious 
fact. Interestingly, this also gives a nice intuitive derivation of the wall crossing 
phenomena. This section is based on the results of Ref. [57] 

5.1. Overview: Topological Strings and Wall Crossings 

Since one of the goals of this section is to relate the BPS partition function with the 
topological string partition function, let us begin with an overview of topological 
string theory. 

The topological string theory gives solutions to a variety of counting problems 
in string theory and M-theory (see Refs. 11541 11551 [T5d1 for review). From the world- 
sheet perspective, the A-model topological string partition function Z top generates 
the Gromov-Witten invariants, which count holomorphic curves in a Calabi-Yau 
3- fold X . On the other hand, from the target space perspective, Z top computes the 
Gopakumar-Vafa (GV) invariants, which count BPS states of spinning black holes 
in 5 dimensions constructed from M2 brancs in M-theory on X ^Hl. Moreover, the 
absolute-value-squared |-Zt op | 2 has been related to the partition function of BPS 
extremal black holes in 4 dimensions, which are bound states of D-branes in type 
II string theory on X 033. 

The topological string partition function Z t0 p also counts the numbers of DO and 
D2 branc bound states on a sing le D6 brane on X, namely the Donaldson-Thomas 
(DT) invariants defined in Refs. I 76 | 77 | The relation between the Gopakumar-Vafa 
invariants and DT invariants was suggested and formulated in Refs. |75j l78l a nd its 
physical explanation was given in Ref. [59] using the 4D/5D connection ^1. More 
recently, a mathematical proof of the GV/DT correspondence was given in Ref. [791 
when X is a toric Calabi-Yau 3-fold. 

However, the number of BPS states has background dependence. As we vary 
moduli of the background geometry and cross a wall of marginal stability, the num- 
ber can jump (this is the wall crossing phenomena discussed in I2.2[) . This means 
that the Kahler moduli is implicitly assumed in most of the above computations, 
and we need to understand at which value of the Kaher moduli the statement holds 
and what happens if we change the value of the moduli. 

In this section we will generalize the results of Ref. [59] to include the background 
dependence of the M-theory computation. We show BPS bound states are organized 
into a free field Fock space, whose generators correspond to BPS states of spinning 
M2 branes in M-theory compactificd down to 5 dimensions by a Calabi-Yau 3-fold 



Crystal Melting and Wall Crossing Phenomena 49 



X. This enables us to write the generating function Zbps of BPS bound states of 
D-branes as a reduction of the square of the topological string partition function, 

Zbps(<7, Q) = Z t o P {q,Q)Z t0 p(q,Q~ ) chambcri 

(5.1) 

in an appropriate sense described in the following, in all chambers of the Kahlcr 
moduli space. Our results apply to the BPS counting for an arbitrary Calabi-Yau 
(whether toric or non-toric) without compact 4-cycles. 

For the conifold, the change of the numbers of BPS states across a wall of 
marginal stability h as been studied by physicists in Refs. [531 [Ml (see also Ref. I157[) 
and mathematicians I 80 | 54 [ The case of generalized conifold geometries was studied 
in Ref. [55] The formula (|5.1[) derived from the perspective of M-theory reproduces 
these results. Our results also provide a sim ple derivation of the "semi-primitive" 
wall crossing formula of Denef and Moore S3 (discussed in Section|H]) , in the present 
context. 

The rest of this section is organized as follows. In Section [S~2l we will explain 
the basic idea: to use M-theory to count bound states of a single D6 brane with DO 
and D2 branes on a Calabi-Yau 3-fold. In Section [531 we will describe the counting 
procedure in more detail and derive the generating function for the numbers of 
BPS bound states using a free field Fock space in any chamber of the background 
Kahler moduli space. In Scction[5]4j we will compare the Fock space picture with the 
known results for the the resolved conifold and its generalizations. We also discuss 
in Section [5~5l the representation of the wall crossing phenomena in the language of 
the crystal melting model. 



5.2. The Basic Idea 

In this section we will explain the basic idea. We will apply this idea, in the fol- 
lowing sections, to find a concrete expression for BPS state degeneracies in various 
chambers for Calabi-Yau 3-folds with no compact 4-cycles. As we will see, this tech- 
nical assumption about compact 4-cycles will ensure that M2 brane particles are 
mutually BPS, and thus the partition function takes an infinite product form F*l. 

We are interested in counting the BPS partition function of one D6 brane bound 
to an arbitrary number of D2 and DO branes. The idea is the following: In M-theory, 
the D6 brane lifts to the Taub-NUT space with the unit charge. D2 branes are 
M2 branes transverse to the S 1 , and DO branes are gravitons with Kaluza-Klein 
momenta along the S 1 . The Taub-NUT space is an S fibration over R , and S 
shrinks at the position of D6. Thus the problem of finding bound states to the 
D6 brane becomes simply the problem of finding BPS states in the Taub-NUT 
geometry. Suppose we have BPS states for flat R 4 ' 1 background. Then for each such 
BPS state we can consider the corresponding possible BPS states in the Taub-NUT 
geometry. For each single particle BPS state we can consider its normalized wave 



See Ref. 158 for cases with compact 4-cycles. 
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functions in this geometry. Such states would constitute BPS states which in the 
type IIA reduction correspond to BPS particles bound to the D6 brane. However, 
this would only constitute single particle BPS states bound to the D6 brane. 

Now consider multiple such particles in the Taub-NUT background. This prob- 
lem may sound formidable, because now we will have to consider the interaction of 
such particles with each other and even their potentially forming new bound states. 
We will now make the following two assumptions: 

Assumption 1: We can choose the background moduli of Calabi-Yau as well as the 
chemical potential so that a maximal set of BPS states have parallel central charge 
and thus exert no force on one another. Therefore, at far away separation, the bound 
states correspond to single particle wave function in the Taub-NUT geometry 
Assumption 2: The only BPS states in 5D are particles. In other words there are 
no compact 4-cyclcs in the Calabi-Yau and thus we can ignore BPS string states 
obtained by wrapping 5 branes around 4-cvcles F^l 

Assumption 1 can be satisfied as follows: Consider the Euclidean geometry of 
M-theory in the form of Taub-NUT times S 1 , where we have compactified the 
Euclidean time on the circle. The BPS central charge for M2 branes wrapping 2- 
cycles of Calabi-Yau, but with no excitation along the Taub-NUT, is given by 



Z{M2) = iA(M2) - C(M2), 



where A(M2) denotes the area of the M2 brane and C(M2) corresponds to the 
coupling of the M2 brane to the 3-form potential turned on along the Calabi-Yau 
2-cycles as well as the S 1 of the Taub-NUT. However we need to include excitations 
along the Taub-NUT. As discussed in Ref. [59] these are given by the momenta along 
the Taub-NUT circle. Let us denote the total momentum along the circle by n (as 
we will review in Section (|5.3[) , this can arise both due to internal spin as well as the 
orbital spin in the SU(2)l C <S0(4) rot ation)- Let us denote the radius of the Taub- 
NUT circle by R. In this case the central charges of the BPS M2 brane becomes 



ay From the viewpoint of geometric engineering | 159 | 160 | y^g means that the corresponding gauge 
theory is Abelian. 
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s 



Z(M2,n) = iA(M2) - C(M2) - n/R. 



(5.3) 



To satisfy assumption 1, wc need to make sure that differently wrapped M2 
brancs all have the same phase for Z. This in particular means that we need to 
choose the Kahler classes so that the 2-cycles of Calabi-Yau have all shrunk to 
zero size, i.e. A(M2) = for all the classes. Even though this may sound singular 
and it could lead to many massless states, by turning on the C(M2) we can avoid 
generating massless states in the limit. The condition that different states have the 
same central charge is simply that 



Note that, in going to type IIA, this condition is simply the statement that the 
B fields are turned on along 2-cycles of Calabi-Yau and the M2 branes wrapping 
them will have B(D2) branes induced. Moreover n, being the momentum along 
the Taub-NUT translates to DO brane charge and as long as the net number of 
branes is positive, they correspond to BPS states of the same type, i.e. preserving 
the same supcrsymmetry. 

Now we are ready to put together all these mutually BPS states as a gas of 
particles in the Taub-NUT geometry. By the fact that they are mutually BPS, they 
will exert no force on one another. Moreover, as long as they are far away, we can 
simply consider the product of the individual wave functions. One may worry what 
happens if they come close together. Indeed they can form bound states, but that 
is already accounted for by including all single particle bound states of M2 branc. 
Here is where the assumption 2 becomes important: If we in addition had 4-cycles, 
then wrapped M5 branes along 4-cycles, which also wrap the S 1 of Taub-NUT can 
now form new bound state with the gas of M2 brane particles on the Taub-NUT 0. 
But in the absence of 4-cycles of Calabi-Yau, we can simply take the single particle 
wave functions (taking their statistics into account) and write the total degeneracy 
of such BPS states, by taking suitable bosonic/fermionic creation operators, one 

az Note that C(M2) is periodic with period l/R. To see this, note that we can view it as a 
holonomy of the gauge field obtained by reducing the 3-form on the 2-cycle of Calabi-Yau, around 
the Taub-NUT circle. The holonomy of a gauge field on a circle of radius R is periodic, with period 
l/R. In terms of the IIA quantities, we have 



where B(D2) is the NS-NS B-field through the 2-cycle in IIA on the Calabi-Yau wrapped by the 
corresponding D2 brane (which has periodicity B — > B + 1). We are denoting by l/R the central 
charge of the DO brane. As such, it does not have to be positive, and in fact it does not have 
to be real either. A better way to think about this is that the quantity wc denote by l/R is 
proportional, up to a complex constant, to the inverse radius of the M-theory circle. For simplicity 
of the notation, we will identify the two, but this fact has to be kept in mind. Note also that relative 
to Rcf. 53 wc arc keeping the D6 branc charge fixed (for example, to Z(D6)=1), and varying the 
DO and D2 brane central charges. 
ba From gravity viewpoint, we have black strings. 



C(M2)+n/R > 0. 



(5.4) 



C(M2) = B(D2)/R, 



(5.2) 
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for each state satisfying C(M2) + n/R > 0. Finally, while the assumption 1 is 
satisfied only for special backgrounds where A(M2) vanishes, the degeneracies are 
guaranteed to be the same everywhere within a given chamber, and are independent 
of this choice. This is all we need to compute all the degeneracies of BPS states in 
various chambers as we will show in the following sections. 

5.3. BPS State Counting and Wall Crossing 

We will use this section to spell out, in a little more detail, how to use M-theory 
to compute the degeneracies of one D6 brane on X bound to D2 branes wrapping 
2-cycles in X and DO branes. The D6-D2-D0 partition function is the Wittcn index 

E 

Tr[(-l) F e- eH ] 

of the theory on 

X x R 3 x Si, 

where we have compactificd the Euclidean time on a circle of radius e. The type 
IIA geometry with one D6 brane lifts to to M-theory on 

X x Taub-NUT x S 1 , 

where the asymptotic radius of the Taub-NUT circle R is related to IIA string 
coupling. Since the D6 brane is geometrized, the computation of the BPS bound 
states of D2 branes and DO branes with D6 brane lifts to a question of computing 
the degeneracies of M2 branes with momentum around the Taub-NUT circle. 
Suppose we know the degeneracies of M-theory in the 

X x R 4 x S 1 . 

This corresponds to taking the R — > oo limit, where the Taub-NUT just becomes 
R 4 (Fig.dSJ). As is clear from the previous section, at fixed B the degeneracies are 
unchanged by varying R since no states decay in the process — all the central 
charges simply get re-scaled. Thus, the knowledge of these allows us to compute the 
degeneracies on Ax Taub-NUT xS 1 background as well. 

The Kaluza-Klein momentum around the Taub-NUT circle gets identified, in 
terms of the theory in the R —> oo limit, with the total spin of the M2 brane. This can 
be understood by comparing the isometries of the finite and the infinite R theory, as 
explained in Ref. 1591 We can view taking R to infinity as zooming in to the origin of 
the Taub-NUT. The isometry group is the rotation group SO (A) = SU(2) L x SU(2) R 
about the origin of K 4 . The SU{2) R is identified with the SO (3) that rotates the 

bb Here we are ignoring the fermionic zero modes in the 4 non-compact directions. Otherwise, 
additional factors need to be inserted to absorb these. See the definition of BPS index in Section fe.ll 
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Fig. 19. The pictorial explanation of 4d/5d correspondence. When we send R — > oo, the Taub- 
NUT space, which is asymptotically R 3 X S 1 as shown on the left, becomes flat IR 4 as shown on 
the right. The M2-brane particles, which are originally bound to the tip of the Taub-NUT, moves 
freely in R 4 after this process. This figure is the reproduction of Fig. 1 in Ref. 1591 

sphere at infinity of the R 3 base of the Taub-NUT. Moreover, the rotations around 
the S 1 of the Taub-NUT, end up identified with the 

U(l) C SU(2) L . 

Thus, the Kaluza-Klcin momentum, is identified with the total spin of the M2 
brane on R 4 . 
Now, let 

A r(m L ,m R ) 

be the degeneracy of the 5- dimensional BPS states of M2 branes of charge (5 and 
spin the intrinsic (2j£,2jjj) = (mi,mjj) (where the spin refers to the spin of the 
highest state of the multiplet). To get an index, we will be tracing over the SU(2)r 
quantum numbers, so we get a net number 

N T L =J2(- 1 ) mR K mL ' mR} ( 5 - 5 ) 

m R 

of 5D BPS states, of the fixed SU(2)l spin m^. These integer invariants are the 
Gopakumar-Vafa invariants of Refs. [60l EU 

Each such 5D BPS particle can in addition have excitations on R 4 . Namely, for 
each 5D particle we get an analytic field 

Hz 1 ,z 2 ) 

on R 4 with z\ i as the complex coordinates. In the usual way, the modes of this field 

${z l ,z 2 )=Y,<xiuiA 1 4 2 (5-6) 

correspond to the ground-state wave functions of the particle with different mo- 
menta on R 4 . (We are suppressing a Gaussian factor that ensures the wave functions 
are normalizable) . Since f/(l) C SU(2)l acts on z\,Z2 with charge 1, the particle 
corresponding to 
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carries, in addition to the M2 branc charge /3 and intrinsic momentum to, a total 
angular momentum 0. 

n = ti + £ 2 + to. (5.7) 

Which of these 1-particle states are mutually BPS? The answer depends on the 
background, and a priori, we need to consider particles in four dimensions coming 
from both the M2 branes and the anti-M2 branes in M-thcory. Along the slice in 
the moduli space we have been considering, the central charges of the particle with 
M2 branc charge (3 and total spin n is 

Z((3, n) = C[fi) + n/R = (B(f3) + n)/R, (5.8) 

where C(/3),B((3) arc the C-field (B-field) flux through the 2-cycle (3. The states 
with 

Z(/3,n)>Q (5.9) 

all preserve the same supersymmetry and bind to the D6 brane (we could have 
picked the opposite sign, and than the particles would bind to anti-D6 branes). See 
Fig.[2ni for illustration. Note that although all the D0/D2 particles are mutually 
BPS, they preserve different supersymmetry from that of the D6-brane, which is 
why we have a stable bound state. 




Fig. 20. The central charges of the D6-brane and the D0/D2-brane particles, expressed in the 
complex plane. The D0/D2 particles, whose central charge lie on the real axis, is stable (unstable) 
in the light (dark) gray region. 



bc In the present case, we are restricting to Calabi-Yau manifolds with no compact 4-cycles. When 
the Calabi-Yau is furthermore toric, as in the cases discussed in Scction r5.4.1l and l5.4.2l the genus 
of the target space curve wrapped by the M2 branes vanishes. This means that the intrinsic spin 
of all the M2 branes vanishes as well. 
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It is easy to analyze the condition (|5.9p . For example, for 

B > 0, R > 

along side M2 branes with (3 > 0, and for sufficiently large n also the anti-M2 brancs 
with (3 < have positive Z > and contribute to the BPS partition function. So 
in general we need to consider both signs of f3. It is important to note that the 
degeneracies N™ of the 5D particles are independent of the background. The choice 
of background only affects which half of the supcrsymmetry the states preserve. 

Now, we can put these all together and compute the BPS partition function in a 
given chamber. Simply, in each chamber, the BPS partition function is the character 
in the Fock space of single particle states preserving the same supersymmetry! In 
fact, a useful way to go about computing the partition function is in the following 
steps: 

Step 1. Start with the unrestricted partition function - the character 

Z Fock = Tr Fock q Qo Q Q * (5.10) 

in the full Fock space. The oscillators of charge (3 and intrinsic spin m and arbitrary 
4d momenta contribute a factor 

J] (1- q e ^+ m Q^ = (1 - q l+m QP) lN ™. (5.11) 

In addition, both the M2 branes, and the anti-M2 branes contribute, and the total 
character is 

oo 

^Fock = n n t 1 - q i+m Q p y N ™ (5-i2) 

0,m 1=1 

Step 2. The 5d degeneracies N^ 1 of M-theory on X x R 4 ' 1 are computed by the 
topological string partition function on X K>0 | 61 | . This allows us to write 

^Fock = Z top (q, Q)Z top (q, Q^ 1 ). (5.13) 

In particular, the knowledge of topological string amplitude allows us to compute 
the BPS degeneracies in any chamber. 

The topological string partition function has an expansion 

oo 

Z top (q,Q)=M( q yW/* [] Ilt 1 -^^)^- ( 5 - 14 ) 

where q and Q are determined by the string coupling constant g s and the Kahlcr 
moduli t by q = e~ 9s and Q = e~* 0. The MacMahon function M(q) was defined 

As compared with Il3.34| l. we are here absorbing the minus sign in l|3.34ll into the definition of 
q. We will use this notation throughout the rest of this paper. Of course, the topological strings 
partition function itself is the same in both notations. 
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previously in (|3.32p . and x(X) is the topological Eulcr characteristic of X. Note 
that topological string involves only the M2 states with positive /3 > 0. On the 
other hand, the full Fock space includes also anti-M2 branes. Since M2 branes and 
anti-M2 branes are CPT conjugates in 5d, this gives another factor of Z top with 

Note that we also have states with j3 = 0. These are the pure KK modes, 
the particles with no M2 brane charge. To count the number of BPS states of 
this type, we note that, for each M 4 momentum we get a classical particle 

whose moduli space is the Calabi-Yau X. Quantizing this, wc get a particle for 
each clement of the cohomology of X. On a (p, q) form SU(2)r acts with the 
Lefshetz action, and SU (2) 1 acts trivially. We get that tcir eigenvalue of a (p, q) 
form on X is mji = p + q — 3. Therefore, the pure KK modes contribute with 
Np=o = — XpO- This agrees with the power of the MacMahon function M(q) we 
get from Z top (q, Q)Z top (q, Q^ 1 ). 

Step 3. We identify the walls of marginal stability as places where, the cen- 
tral charge vanishes for one of the oscillators contributing to Z top (q,Q) or 
Ztopiq^Q^ 1 )- 

Step 4. In any chamber, the BPS partition function is a restriction of ZFock to the 
subspace of states that satisfy Z(f3,n) > in that chamber. 

-Zbps (chamber) 

^Fock I chamber ; . , 

(5.15) 

chamber 

There is a simple way to keep track of the chamber dependence. For the book- 
keeping purposes, it is useful to identify the central charge with the chemical po- 
tentials. Then, in a given chamber, the BPS states are those for which 

q n Q 13 < 1 (5.16) 

where n = l+m = £i+£2 + m is the total spin (|5.7j) . As we vary the background, and 
cross into a chamber where this is no longer satisfied for some (n,/3) in Z top (q, Q) 
or in Z top (q, Q _1 ), we drop the contribution of the corresponding oscillator. 

For example, consider some special cases. When 

R>0, B^oo, (5.17) 
for all Kahler classes, Z((3, n) = (j3B + n)/R > implies that 

/3 > 0. 

In this case, only M2 branes contribute to the partition function. This is the chamber 
discussed in Ref . [SHI By taking the limit (|5 . 1 7[) in (|5.15[) , we find 

Z BPS (R > 0, B 00) = Z BT (q, Q) = M(q)*W/ 2 Z top (q, Q). (5.18) 



bc This is a degenerate limit of l|2.16ll . 
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The partition function in this chamber computes DT invariants. In Rcf. I79[ it was 
shown that, for a toric Calabi-Yau, Zbps is equal to Z top up to a factor which 
depends only on q. Here we derived the relation between Zbps arL d Z top including 
the factor oiM{q)^ x ^ 2 . 

On the other hand, when < B <C 1, the BPS partition function is given by 

Zbps(<7, Q) = Zncdt(<7, Q) = Z top (q, Q)Z top (q, Q~ x ). (5.19) 

This gives the non-commutative DT invariants studied in Section[3] When X is 
toric, the partition function is com puted b y the crystal melting picture E 4 | 32 | g en _ 
eralizing the previous result of Refs.ESESfor C 3 . In ScctionH we have seen that the 
thermodynamic limit of the partition function of the crystal melting model gives the 
genus-0 topological string partition function. This result was mysterious since the 
relation between Z top and Zbps was supposed to hold in the DT chamber discussed 
in the previous paragraph, not in the NCDT chamber. We now understand why 
there is such a relation in the non-commutative DT chamber also as in (|5.19l) . 

5.4. Examples 

In this section we give some examples of geometries without compact 4-cycles. We 
first study toric cases, namely resolved conifold and generalized conifolds. We also 
give a simple example of a non-compact, non-toric Calabi-Yau manifold as well. 
In each of these cases, we will use our methods to lay out the chamber structure, 
identifying walls where BPS states jump, and the BPS partition function in each 
chamber. In some of the cases we study, the jumps were studied by other means. 
We will show that they agree with the M-theory results. 

5.4.1. Resolved Conifold 

The topological string partition function for the resolved conifold was given previ- 
ously in (|3.33j) . which means that the only non- vanishing Gopakumar-Vafa invari- 
ants are 

iV° =±1 = l, iV° =0 = -2, (5.20) 

and that all BPS states in 5 dimensions has no intrinsic spin I 60 | 6 1 |_ Q ur f ormu i a 
(|5.15p then implies that BPS states are counted by 

Zbps(<7, Q) = Z t o P {q,Q)Z top (q,Q~ )| 

chamber 

II {i-q n Q f) ) nN °?- ( 5 - 21) 

(P,n):Z(p,n)>0 

The product is over /3 = 0, ±1 and n = 1, 2, ... such that Z(/3, n) > 0. 

The chamber structure is easy to identify in this case since the Kahler moduli 
space is one-dimensional. When 

R > and m - 1< B < m (5.22) 
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with some m > 1, the formula (|5.21[) gives 

oo oo 
Z B P S (q,Q) = M(q) 2 J](l II C 1 - ^Q^T- 

71—1 n— 771 

In particular, the chamber at m = oo counts the DT invariants ^§1, 

oo 

Z B P S (q,Q) = z DT (q,Q) = M( q ) 2 Y[(i -q n Q) n , 

n=l 

while the chamber at m = 1 counts the non-commutative DT invariants I 

oc oo 

Z B p S (q,Q) = Z NCDT (q,Q) = M(q) 2 JJ(1 -q n Q) n JJ(1 - g n Q _1 ) n . 

n— 1 n— 1 

On the other hand, when 

i? < and - m - 1 < S < — m 

with m > 1, we have 

m 

^BP S (g ) Q) = n( 1 -« n< 9)"- 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



In p articular, the chamber at m = oo counts the Pandharipande- Thomas invariants 



in p 



ZBPs(q, Q) = Z PT (q, Q) = fl(l- q n QY 



(5.28) 



n=l 



Summarizing, wall crossing of the resolved conifold is essentially a one- 
dimensional problem (see Fig. l21[) M . For later reference, we defined by chamber 



NCDT chamber 



DT chamber PT chamber 



Co Ci . C- 



C c 



C c 



Core Region 

C 2 Ci Co 



D2+D0 D2+2D0 D2+3D0 D2+nD0 D2+nD0 D2+3D0 D2+2D0 D2+D0 

R>0 R<0 

Fig. 21. Chamber structure for the resolved conifold. The charge 72 for the decay pattern 7 — > 
71 + 72 is shown for each wall. 



Strictly speaking, only half of the chambers are shown in Fig. 1211 There are also chambers for 
the flopped geometry, which correspond to R > 0, —n — 1 < B < — n and R < 0, n — 1 < B < n. 
This means that the there are two different ways to go from the core region to the NCDT chamber, 
and the chambers actually line up on the S 1 , rather than on R as depicted in Fig. 1211 
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C n and C n as a region where R < 0, — n — 1 < B < —n and R> 0, n — 1 < B < n, 
respec tive ly. The partition function is trivial in the chamber Co (the so-called core 
region S3) ; and each time we cross the wall we obtain one factor. On the left most 
side, we have Zncdt(9, Q) = Z top (q,Q)Z top (q,Q^ 1 ). 

As we will see in the next section, these results agree with the predictions from 
the wall crossing formula in all chambers. 

5.4.2. Toric Calabi-Yau without Compact ^-cycles 

We can also test our formula (|5.15|) for a more general toric Calabi-Yau without 
compact 4-cycles. A toric Calabi-Yau is characterized by a convex polygon on a 
square lattice, and the abse nce of compact 4-cycles means that there is no internal 
lattice point in the polygon I 30 | 31 [ By SX(2,Z) transformations of the lattice, one 
can move one of the edges of the polygon along the positive a;-axis, and one of the 
vertices to (x, y) with —y < x < 0. If we require that there is no internal lattice 
point, there are essentially two possibilities: (x,y) = (0, 1) and (0,2). In the former 
case, the polygon is a trapezoid of height 1, and the corresponding Calabi-Yau is 
the so-called generalized conifold, which has L — 1 P l! s where L is the area of the 
trapezoid. We will describe the resolved geometry in more detail below. In the latter 
case, we have a isosceles right triangle with two legs of length 2, which corresponds 
to C 3 /Z 2 x Z 2 . 

For the generalized conifold, the topological st ring partition function has been 
computed in Ref . 11611 using the topological vertex u^SI The counting of BPS states 
has been carried out in all chambers in Ref. [82] Thus, we will use this case to test our 
formula (|5.15[) . For C 3 /Z 2 x Z 2 , the counting in the non-commutative DT chamber 
has been done in Ref. I8T1 

Homology 2-cyclcs of the generalized conifold correspond to the simple roots 
cti, • • ■ , of the Al-i algebra. To identify them in the toric diagram, we divide 
the trapezoid into L triangles of area 1 and label the internal lines dividing the 
triangles as i = 1, • • • , L — 1. Each line i corresponds to the blowing up P 1 at c^. 
We will denote the D2 charge by 



In general, there are several ways to divide the trapezoid, and they correspond to 
different crepant resolutions of the singularity. If the two triangles across the line 
i form a rhombus, we have a resolution by 0(— 1, —1). On the other hand, if the 
two triangles form a triangle of area 2, the resolution is by 0(— 2,0). Both the 
topological string partition function and the BPS counting depend on the choice of 
the resolution. 

The topological string partition function for this geometry is given by 




(5.29) 



OO 



Z top (q, Q) = M(q) L / 2 I J(l - q n QiQ 



i+i ■•■Qj) 



71 N, 



(5.30) 



n — 1 i<j 
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where 

N .. = (_i)i+#{fce/|i<fc<i} ; (531) 

and / is the set of internal lines of the toric diagram corresponding to the resolution 
by 0(— 1, —1). Thus, the only non-vanishing Gopakumar-Vafa invariants are 

^ ro =±%,..^) = (5.32) 

and 

Nf=° = -L. (5.33) 

All BPS states in 5 dimensions carry no intrinsic spin. 
The central charge Z(/3,n) is given by 

Z{l3,n) = R- 1 (n + y^n l B(a = i)\ . (5.34) 



The formula (|5 . 1 5[) predicts that BPS states in the chamber characterized by B^s 
are counted as 

Z BPS (q,Q) = M(q) L [] (l-q n Q^. (5-35) 

(p,n):Z(P,n)>0 

Here the product is over all roots /3 of Am-i and n = 1,2, ... such that 2/(/3, n) > 0. 
This agrees with the result in Ref. |82] 



5.4.3. .A Non-toric Example 

Our discussion in Scction l5~2l and 15.31 are not limited to toric Calabi-Yau's, and 
applies to any Calabi-Yau without compact 4-cycles. In order to illustrate this point 
in a concrete setting, let us describe the geometry shown in Fig. 1221 This geometry 
arises by identifying two of the four external legs of the (p, q)-web of the resolved 
conifold. This is one of the simplest the non-toric (albeit formally toric) geometries 
studied in Ref. 11621 and it is straightforward to repeat the following analysis to 
other non-toric geometries discussed in Ref. 11621 




Fig. 22. The non-toric Calabi-Yau which arises by identifying two external legs of the (p, q)-web 
of the resolved conifold. 
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In addition to the P 1 of the resolved conifold, the geometry of Fig.[22]has another 
compact P 1 which arises from identification. Let us denote their homology classes 
by /3 riginai and /3 n ew, respectively. As a basis of the homology class, we choose 

Pi = /^original an d fa = /^original + Pnew 

The topological string partition function is given by LLs£al 

w..ft.«- w fna-^)-) ft ' 

(5.36) 



k.n—1 



where Q\ and Qi are the variables corresponding to fa and fa. The Gopakumar- 
Vafa invariants are therefore given by 

A$ = o = -2, ^=1, Ni 0i+kp2 =l, Nj&=-1 (keZ\{0}). (5.37) 

Notice that genus 1 Gopakumar-Vafa invariants are nonvanishing in this non-toric 
example. This arises because our geometry has an extra cycle compared with the 
resolved conifold. 

Again, the general formula gives (notice that m ^ in this case) 

Zbps(<7, Q) = Z top (q, Q)Z top (q, Q~ )| 

chamber 

II (l-q l+m Q^) lN P . (5.38) 

(l3,l,m):Z(P,n=m+l)>0 

The formula for the central charge is 

Z L8 = J2 mPi, n \ = R- 1 (^2 n.Bifa) + n^j , (5.39) 
and the equation 

Z{fi, n = I + m) = 0, (5.40) 

with corresponding Gopakumar-Vafa invariants nonvanishing, determines the posi- 
tion of walls of marginal stability. This is a new result which has not been discussed 
in the literature to the best of our knowledge. 



5.5. Wall Crossing and Crystal Melting 

We finally explain how the wall crossing phenomena described in this section is 
recast in the language of crystal melting described in Scction[3l Let us describe the 
example of the resolved conifold for simplicity; more general discussions, applicable 
to any chamber for generalized conifolds, is included in Section[7] as a special case 
of more general open BPS invariants. 

The crystal melting model described in Scction|3] is for the chamber Co (see 
Fig.rTT]). In the chamber C n , we have to change the ground state of the crystal. The 
new ground state (see Fig. [23]) has n + 1 atoms on top. 
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Fig. 23. The conifold crystal in the chamber C*2. We have 3 atoms on the top. Compare with 
Fig.QT] 

Using the same melting rule, we have the function Crystal, c„ • I n Ref- HOI it was 
shown combinatorially (see Section [S~51 for derivation from vertex operators) that 
the partition function takes a simple form 

A:rystal,C„ (.<2 I ?1 J ~ M WO ?1 ) 

n (> + (#w>)-r n (i + ( ? S"')'( 9 !"»)*«^ m " , ""'° , (5 - 41) 



X 

fe>l fe>l 



where <7q { are weights in the n-th chamber C„ . If we take 

g^=^ + V, ^ n) =%"gr +1 , (5.42) 

and substitute back into ([5.41)1 . we have 

M(qm) 2 U(^ + qU' 1 ) k l[0- + ioQi +1 ) k , (5-43) 

k>n fe>0 

which coincides with the result of (|5.23p under the parameter identification ([3.211) . 
This means that the wall crossing is represented as (1) a change of the ground state 
and (2) change of the weight in the crystal melting model. This is consistent with 
the discussion at the end of Section l373l where different choice of the fractional brane 
leads to Seibcrg duality of the quiver. The fact we have 3 atoms on the top of the 
crystal corresponds to 3 arrows starting from the D6-brane node in Fig.[8l 

When we send n — > oo, the ridges of the crystal approaches to the (p,q) web 
of the resolved conifold with larger and larger P 1 inside, and we begin to see two 
different corners of the crystal, each corresponding to the one for C 3 . This is how 
the combinatorial expression ([5.41)1 has been derived, and clarifies the connection 
between the crystal melting model of Section[3] and those in the topological vertex 
model. In other words, the limit n — > oo of our crystal gives the crystal melting 
model appearing in topological vertex constructions. 
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We can also construct crystal melting models for the chambers C n , where the 
partition function is a finite product. See Refs. [Ml [551 for discussion. 

5.6. Wall Crossing and Free Fermions 

There is yet another representation of the wall crossing phenomena, from the view- 
point of free fermions^. This is well-known for the C 3 crystal; the topological vertex 
can be computed in the transfer matrix approach^. I n this section we arc going to 
see that this can be generalized to our more general toric geometries. We moreover 
see that this gives anot her perspe ctive on the wall crossing phenomena . This section 
is based on the papers I 81 | 84 | 85 [ p or simplicity we specialize the conifold example 
in this section. More complete discussion, including the open BPS invariants, will 
be given in Section [7T4l See also Ref. 11671 for combinatorial description. 

We begin with two-dimensional free fermions 0, which have the mode expansion 

Hz) =$>fc + i/ 2 z- fe -\ r{z) = (5.44) 

k k 

satisfying the commutation relations 

H'k+1/2, ^-1-1/2} = (5-45) 
Starting from the vacuum |0) annihilated by all positive modes 

^+1/2 10> - r k +i/2 10> =0, k > 0, (5.46) 
we can construct a state 

L") = ]^[ V 1 * ai -i/2^-f»i-i/2|0), ai=m-i, bi = /4-i, (5.47) 

i 

corresponding to a two-dimensional partition /i = /i 2 . . .) o The mode expan- 
sion of the bosonized field d<f)(z) =: ip(z)ip*(z) : 

d4>(z) =^a„z-"- 1 (5.48) 

n 

satisfy the Heisenberg commutation relation 

[a n , a m ] = n6 n +m,o- (5.49) 

bg See Refs. 163 161, 165, 166]for connections between topological string theory and free fermions. 

bh Sce Refs. 168 169 for the formalism used here. 

bl See [Appendix A.5| for basic definitions about partitions. 
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Define vertex operators 



and 



q Tn 



r±(<7)=exp^ a Tn , (5.50) 



rl(g)=cxp^ V ' n a Tn , (5.51) 

n 

They have the property ^JJ 

r_(g)|AX> = ]T 5 lMh|A| l A )' r + (g)| A .) = ^ ^H A l|A), (5.52) 

and 

r'_(g)[/i>= E 9 HHA| |A), r+(g)|^> = £ ^ HA| |A>, (5.53) 

AM/i* A'Vjtx' 

where A >- fi denotes the interlacing condition 

Ai > /ii > A 2 > ^2 > • • ■ • (5.54) 
From the commutation relation (|5.49|) . we have 

F_(p),r+(?)] = i+p-\, [r'_(p),r + (?)] = 1 + ^^, 

1 1 (5 55) 

[r_(p),r + (g)] = - — — , [r'_(p),r' + ( q )} = - — — 

1 — p L q 1 — p L q 

where the commutator means the multiplicative commutator [x,y] = xyx~ 1 y~ 1 
here. This can be concisely summarized as 

n(px),TZ(p 2 )} = (1 - e^tfpZ)-^ 6 ^ (5.56) 

for ei.2 = ±l,ii.2 = ±1, if we write r±(p) = T±(p) and T±(p) = r ± (p). 

Now the claim is that the BPS partition function in non-commutative 
Donaldson-Thomas chamber can be expressed in the following form (as before, 
we define q = <7o<7i) : 

Zbps (90, <?l) = <0| • • ■ T_ (q- 2 )T'_ (q- 2 qi )T_ (q- 1 )^ {q- 1 q 1 ) 

xr + (i)r;( gi )r + ( g )r;( ggi )...|0), 

i.e. 

^bps(<Zo,<Zi) = (0| • • .A-{q- 2 )A-{q- x )A+{l)A+{q) . . . |0), (5.58) 



b 1 We follow the notation of Ref . 1841 which is different from those of Ref . 1811 1851 The notation here 
makes it easier to generalize our results here to those in Section r7.4.2l Compared with Ref. 1811 we 
have 

r here^ = rf eIe (q), r\ eIe (q) = r^ CIC (q' 1 ), r'i lcrc (q) = r^ harc (q), r^ crc (<j) = r'* crc (g- 1 ). 
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when wc introduce 

A+(x) = r + (x)T' + (x qi ), A_(y) = T_(y)T'_(y qi ) (5.59) 
This can also be expressed in a suggestive form 

Z BPS = (fi_|fi+) (5.60) 

when we introduce 

(0_| = (0| . . . A^q-^A^q-^A^q- 1 ), |fi+) = A + (l)A + (q)A + (q 2 ) . . . |0). 

(5.61) 

It is easy to check f|5 -58[) directly by explicit commutation. Indeed, from (|5.55[) 
we have 

\A (t\ A , li,W = 

(l-x-^y) 



[A_(x),A+(y)] = - L—(l + qiX -iy)(l + q -i x -i y ) (5.62) 



and therefore 

oo oo 

l[A_(q- m ),~[[A + (f 

.m=l n — 



oo oo 



= n n (1-^2 (1 + ™ m+n ){i + €\ m+n ) 

m=l n=0 V y ' 

oo 

n=l 

This is precisely the BPS partition function of the conifold. 

Of course, the next natural question is where the expression (|5.57p comes from. 
The answer is that it arises from a sequence of Young diagrams obtained as the 
slice the crystal. As shown in Fig. [M] we can slice the conifold crystal by an infinite 
number of parallel planes, and we have a set of Young diagrams {A(n + l/2)} ng z, 
where we choose the label the planes (and therefore Young diagrams) by half- 
integers. By looking at the crystal structure carefully the melting rule (|3.16[) is 
translated into the condition that Young diagrams satisfy the interlacing conditions 

. . . A(-5/2) ^ A(-3/2) < A(-l/2) ^ A(l/2) J A(3/2) £ A(5/2) y A(7/2) 

(5.63) 

where we used the notation A -< pL (A -{ fj) for A ~< ^ (A* -< /U*) , respectively. The 
vertex operators T± and r' ± in (|5.5T[) implement (in the language of free fermions) 
the interlacing conditions (|5.63[) . 

It is also easy to generalize this analysis to other chambers. For example, in the 
chamber C2 where the crystal takes the form shown in Fig.[23j the analysis of the 
melting rule for the crystal suggests that we consider 

z B p S (9o,9i) = (0|...r_( g - 1 )r'_( g - 3 gi )r + (i)r'_(< z - 2 qi ) 

, , (5.64) 

xT + (q)T_(q- 1 q 1 )T + (q 2 )T + (q 1 )r + (q 3 )r + (qq 1 ) . . . |0), 

We can again easily verify explicitly that this coincides with (|5.23j) . Note that the 
vertex operators appearing in (|5.64j) are the same as in (|5.57j) . and we are just 
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Fig. 24. Slice of the conifold crystal by parallel planes. For a molten crystal, each of the slice 
gives a (complement of a) Young diagram. 

changing their orders! This means that wall crossings are described by commuting 
vertex operators! We will see in Section l7.4.2l that they hold more generally. 

6. Wall Crossing Formulas 

In the previous section we gave a derivation of the wall crossing phenomena by 
using the M-theory viewpoint and reformulating the BPS state counting problem 
as a counting problem of free spinning M2-branes. We have seen that when we cross 
a wall we obtain/lose one factor 

Historically, the M-theory argument was not the first method to derive these 
results. There is an independent derivation of the results, using the wall crossing 
formulas. The aim of this section is to explain the wall crossing formula, and to 
apply it to the examples we studied in the previous section. 

6.1. Overview: What Wall Crossing Formula Tells Us 

Let us first explain the meaning of the wall crossing formula before discussing the 
explicit form of it. Suppose that we already know the positions of walls of marginal 
stability (for example, by the central charges (|2.16j) ). Consider two chambers sepa- 
rated by a wall. Due to the wall crossing phenomena, the BPS index 0(7) in general 
jumps along the wall. The wall crossing formula relates the BPS indices on both 
sides of the wall. In other words, given an index on one side of the wall, wall crossing 
formula tells us the index on the other side of the wall. 

This means that given the chamber structures and the BPS index in one chamber 
of the moduli space, wall crossing formula reproduces the BPS index in all the 
other chambers of the moduli space. Of course, we still need to know the BPS 
index in at least one chamber. In the example of the resolved conifold discussed 
in Section l5.4.1l as well as in the examples of generalized conifolds in Section l5.4.21 
there is a special chamber (core region) where the BPS partition function becomes 
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trivial, and starting from this chamber we can reproduce the partition function in 
all other chambers. 

In the literature there arc several different versions of the wall crossing formula. 
The first one is the primitive/semi-primitivc wall crossing formula, which is proposed 
by Denef and Moore in the context of supcrgravity and is discussed in detail in the 
next section. The second is the Kontsevich-Soibelm an wall c rossing formula SSI. 
There is a lso yet another formula by Joyce and Song | 170 | 171 | s ee a j so the recent 
paper E21. In this section wc only discuss the semi-primitive wall crossing formula, 
and comment on the more general Kontsevich-Soibelman formula in |Appcndix A.3| 



6.2. Semi- Primitive Wall Crossing Formula 

Let us begin with the primitive wall crossing formula. Suppose that at the walls of 
marginal stability the decay pattern is 7 —> 71 + 72, where 71 and 72 are primitive 
charges. Here primitive means that particle with that charge is stable at any value 
of the Kahler moduli. In this case, the primitive wall crossing formula says that the 
jump of 0(7) is given by 

Afi( 7 ) = (-1)^^|( 7i ,7 2 }|q(7i)Q( 72 ), (6.1) 

where ( , ) refers to the symplectic paring of the charge lattice T. We do not need 
to specify the value of Kahler moduli for f2( 7 i) and ^(72) since we assumed that 
7i , 72 are primitive. 

Semi-primitive wall crossing formula deals with the decay pattern 7 — > 71+^272- 
The formula says that 



fi( 7 i) + ]T AO( 7l + N 2 j 2 )p N2 = n( 7 i) II f 1 (-i) fc<71i7 V 



N 2 >0 k>0 



fc|<7l,72>|0(fc 72 ) 



(6.2) 



where p is a formal variable. Note that the term constant in p gives a trivial relation, 
and the linear term gives the primitive wall crossing formula. 



In general, the decay pattern is 7 — ¥ N\ji + -/V272 (see the comment below 
(|2.16j) in Section l2~2j) . and we need a more general wall crossing formula, the one 
proposed by Kontscvich and Soibclman. However, the story simplifies considerably 
for a noncompact Calabi-Yau manifold with a single D6-branc. Since we are counting 
D6/D2/D0 bound states, 7 has one D6-brane charge. The fact that the Calabi-Yau 
manifold is non-compact with a single D6-brane means that it is impossible for a 
D6-brane to decay into an anti D6-brane and two D6-branes, for example; that will 
cost infinite energy. This means that either 71 or 72 has D6-brane charge one, while 
the other has charge zero. Without losing generality we can assume that 71 contains 
a D6-brane, while 72 contains none. This means that the decay pattern is always 
7 — >• 71 + N 2 j2- Note that this is not the case if we have multiple D6-branes. 
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Two further simplifications occur in our setup. First, the wall of marginal sta- 
bility is determined by 72 alone. As we have seen in (|2.f 6[) . the wall of marginal 
stability is determined by 

ArgZ( 7l ) = ArgZ( 72 ). (6.3) 

However, since 71 contains D6-branc having an infinite mass, we have 

ArgZ( 7 i) ~ ArgZ(D6), (6.4) 

and the 71 dependence drops out from the story. 

Second, the values of the Gopakumar-Vafa invariants given in Section f5 .4.11 and 
15.4.21 ( (|5.20[) and (|5.32[) ) says that Gopakumar-Vafa invariants vanish when we 
have more than one (anti) D2-branes. This means that when 72 contains (anti) D2, 
semi-primitive wall crossing formulas simplifies to 



fi(7i)+ J2 A^( 7 i + A^272)F 

N 2 >0 



N 2 



«(7i) II V 1 (-l) fe<7l ' 7 V) (6.5) 



fc>0 

(71,72)10(72) 



0( 7l ) (l - (-!)<™>p) 



Let us now use apply these formulas to the resolved conifold, whose chamber 
structure is summarized in Fig. 1211 We start with the simplest example, the wall 
between C and C%. At C (core region), the degeneracies are simply 

0(D6) = 1, all others zero, (6.6) 

and the decay pattern at the wall of marginal stability is 7 — > 71 + -/V272 with 

7i=D6, 72 = D2 + D0, (71,72) = 1. (6.7) 

This means, when combined with the expressions for the Gopakumar-Vafa invariants 
(15 . 20)) and the wall crossing formula (|6.5I) . 

fi(7i + 72! Ci) = 1, fi(7i + N 2l2 ;C\) = for all N 2 > 1. (6.8) 
and we have 

Z B ps(q,Q;C 1 ) = l + qQ = l-(-q)Q. (6.9) 

We can repeat the same argument for the next chamber C 2 - At C\ , the degen- 
eracies are 

fi(D6) = 1, fi(D6 + D2 + DO) = 1, all others zero, (6.10) 



In this section we are concentrating on toric examples, and do not discuss the non-toric example 
of Section l5X3l 
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and the decay pattern at the walls of marginal stability is 7 — » 71 + N2J2 with 

71 = D6 or D6+(D2+D0), 72 = D2 + 2D0, (71, 72} = 2. (6.11) 
Therefore, (|6.5|) and Gopakumar-Vafa invariants (|5.20[) tell us that 

0(7i+7 2 ;C 2 ) = -2, fi( 7l +27 2 ;<5 2 ) = l, 

0(71 + ^272; Ci) = for all N 2 > 2, 

and we have 

fi(D6) = fi(D6 + D2 + DO) = 0(D6 + 2D2 + 4D0) = f2(D6 + 3D2 + 5D0) = 1 
fl(D6 + D2 + 2D0) = 0(D6 + 2D2 + 3D0) = -2. 

(6.13) 

The partition function is then given by 

^bps(<7, Q] C 2 ) = H (1 - i-qfQ) 1 . (6.14) 

i=l,2 

In general, on the wall 7 — > 71 + ^272 with 72 = D2 + n DO, we have from (|6.5[) 
£fi( 7 i)g 7l +A ]T fi(7 1 +7V 2 72)9 7l p Ar2 =E9 7l ^(7i) (!-(-!)»", 

71 JV2>0,7i 71 

(6.15) 

where q = (q, Q) collectively refers to the DO-brane chemical potential q and the 
D2-brane chemical potential Q. If we set p = q 1 ' 1 = q n Q we have 

^bps(<Z, Q; C n ) = Z BPS (q, Q; C n _i) (1 - (-l)'VQ)" . (6.16) 
This means in general 

n 

Z B ps(q,Q;C n ) = IX (1 - (-i)V'Q) 3 ' ■ (6.17) 

and in particular 

OO 

^bps(9, Q; <?«,) = n C 1 (-i)V'Q) j • (6.18) 
J'=l 

The chamber between Coo and Coo is special, since this is the wall where the D6- 
brane makes a bound state with pure DO-branes without D2-branes, and therefore 
we should use the original formula (|6.2j) instead of (|6.5[) . The wall crossing formula 
(ECU) gives 

OO 

^ B ps(g, Q; Coo) = M(-q) 2 J! (1 - (-i)V'Q)' ■ (6-19) 

Repeating similar analysis, we have 

OO CXD 

^bps(<Z, Q; C n ) = A/(-<z) 2 J] (1 •: L)V'Q)' [] ( X ("WO* ■ ( 6 - 20 ) 
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and in particular 

oc oo 

Z B p S (q, Q; Co) = M(-q) 2 [] (1 (-l)V'Q) j II ( X - (-ly^'Q" 1 )'' ■ (6-21) 

3=1 i=i 

These results coincide with those of Section[5](see Section [5.4.ip . The same argument 
applies to the generalized conifolds discussed in Section l5.4.2l 

6.3. Derivation of the Wall Crossing Formulas 

In the previous section we have seen that the semi-primitive wall crossing formula 
correctly reproduces the results obtained in Section[SJ However, we still do not 
see the physics behind the formula. In this section we give a derivation of the 
semi-primitive wall crossing formula following Ref. 1471 This will uncover tantalizing 
connections with the multi-centered black hole solutions described in |Appcndix A.2| 
In supergravity, a particle with primitive charge is described by a single-centered 
black hole. The BPS state with charge given by the sum of primitive charges, 7 = 
71 + -/V272, correspond to a multi-centered black hole. Since 71 contains the D6- 
brane and has an infinite mass, 71 is called the "core" of the black hole. The D0/D2 
particles, which are sometimes called "halos" , surround the core of the black hole 
(see Fig.ESl). 




Fig. 25. The D6-brane (the core) sits in the center of the black hole, and is surrounded by 
D2/D0-brane charges (halos). 

Let us begin with the primitive wall crossing formula. Suppose that on one side 
of the wall we have a stable bound state with charge 7 = 71 + 72 , where 71 , 72 
are primitive. Suppose moreover that this particle become unstable on the other 
side of the wall. In supergravity, the particle corresponds to a 2-centered black 
hole, whose distance between the two centers given in (|A.7j) . Since (|A.7j) contains 
the central charge, the distance depends on the value of the Kahler moduli. The 
distance between the two centers become infinite on the wall of marginal stability; 
this means that there is no longer a stable bound state, and the index jumps. 

Since the two centers are far away from each other on the wall, the Hilbert space 
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splits into direct products: 

AH(7 71 + 72) = [J12] ® ^(71) ® W(7a), (6.22) 

where we have included the rotational degrees of freedom coming from the angular 
momentum of the 2-centered black hole (|A.12[) . This means 

Afi( 7 ) = (-l) 2Jl2+1 (2J 12 + l)0( 7l )0(7 2 ) (6.23) 

which is the primitive wall crossing formula, when we use the expression for J\i in 
(lA~T2l 

We can also derive the semi-primitive wall crossing formula (|6.2p from super- 
gravity viewpoint. The formula (|6.2[) is written in the form of a generating function, 
and contains several terms when expanded with respect to p. For example, the 
quadratic term for p contains two different contributions, one with two factors of 
^(72), and one with a single 0(272). What this means is simple: we can distribute 
the charge 272 into two particles each with charge 72, or to a single particle with 
272. Since the distance between the centers in (|A.7[) depends on the charges, the 
corresponding halos have different radius. 

Consider the more general situation 7 — > 71 + 7V 2 72- We can again distribute 
the charge -/V272 to halos with different radius, each having £^72. Since each halo 
has mutually zero symplectic pairing, we can ignore the interactions and use the 
2-centered black hole solution for each halo 0. This situation is concisely expressed 
if we introduce the generating function: 

@P N * An at =^(7i)(g)^(/[^7 2 ](gm72)), (6.24) 

^ 7->7i+-/V 2 72 k 

where T (p k [Jj 1 ,kj 2 ] ^(72)) represents a Fock space spanned by halos with 
charge kj2, and (J), (^) represents a formal sum and a product, respectively. This is 
a generalization of (|6.22l) . 

Now the each halo [</ 7l .fc 72 ] ® ^(72) contributes a factor 

I _ ( — 1) fc (7l,72)pfc I 



The sign comes from the angular momentum fc(2Ji 2 + 1) = k 1(71,72)!, and the 
power is determined by the degeneracy fc(2Ji2 + 1)^(72). This immediately gives 
the formula 



o( 7 i) + y, + ^72)^ = 0(71) n ( x (-i) fc<7i ' 72> /) 

7V 2 >0 fc>0 

(6.25) 

which is the semi-primitive wall crossing formula (|6.2j) . 



'The two assumptions in Scction[5]2] are the underlying assumptions behind this derivation. 
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As this derivation shows, it is not a coincidence that the methods of Scction[5] 
and those of this section give the same answer. The derivation of the semi-primitive 
wall crossing formula uses the Fock space spanned by D2/D0 halos, whereas when 
lifted to M-theory we have a Fock space spanned by spinning M2-branes and the 
graviton. We should keep in mind, however, that the argument of Scction[5]does not 
require the validity of supcrgravity approximation. 

7. Open BPS Wall Crossing 

In this section, we generalize the analysis in previous sections to include "open 
BPS invariants". The BPS invariants in previous sections, which should be called 
"closed BPS invariants" more precisely, count D0/D2/D6-branes wrapping holomor- 
phic 0/2/6-cycles. In this section, we also consider D2-branes wrapping disks whose 
boundaries are on D4-branes wrapping non-compact Lagrangian 3-cycles. The open 
BPS invariants counts the degeneracy of such D-brane configurations. Since D4- 
branes break the supersymmetry to half and fill K 1,1 in R 3 ' 1 , from gauge theory 
viewpoint we are discussing the wall crossing in the two-dimensional Af — (2, 2) 
theory. 

The notion of open counting exists already in the topological string theory; of 
course, this is the open topological string theory. One of the most important results 
in the open topological string theory is that the partition function can be computed 
by gluing the topological vertex ^1, assigned to each trivalcnt vertex of the (p, q)- 
web. 

Motivated by this result in the topological string theory, we propose a general- 
ization of the topological vertex, which we call the "non-commutative topological 
vertex" and which capture the open BPS invariants in the sense mentioned above. 
The vertex is defined combinatorially using the crystal melting model proposed 
recently, and depends on the value of closed string moduli at infinity. The vertex 
in one special chamber gives the same answer as that computed by the ordinary 
topological vertex. 

We prove an identify expressing the non-commutative topological vertex of a 
toric Calabi-Yau manifold X as a specialization of the closed BPS partition function 
of an orbifold of X, thus giving a closed expression for our vertex. We also clarify the 
action of the Weyl group of an affinc Al-i Lie algebra on chambers, and comment 
on the generalization of our results to the case of refined BPS invariants. 

We also derive the wall crossing of open BPS invariants by the generalization of 
the M-thcory argument in Scction[5] This reproduces the results from the crystal 
melting model, and moreover predicts the wall crossing of open BPS invariants 
with respect to open as well as closed string moduli. We also comment on the 
consistency of this result with the work of Ref. [36l which discuss the wall crossing 
in two-dimensional Af — (2, 2) theories. 

This section is based on the papers Refs.[62j[63] See the paper Ref. IMlfor geomet- 
rical meaning of the open BPS invariants and the vertex operator representations. 
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The vertex operator representation for closed BPS invariants is also discussed in 
Ref.EU 

7.1. Overview: Open BPS Invariants 

We begin this section by reformulating the closed BPS invariants for generalized 
conifolds (Scction [7.2|) . In this formulation, the closed BPS partition function l bm 

7 C 

^BPS^cr'^')' 

which is defined in Section [2~2l as the generating function of the degeneracy of D- 
brane BPS bound states, depends on maps a' , 9' specifying a chamber in the Kahlcr 
moduli space0. The use of these maps is natural from the viewpoint of affine Al-i 
algebra. As we have seen in Section lOl (see the discussion around (|5.17|l ). in one 
special chamber Ct op of the Kahler moduli space, the BPS partition function is 
equivalent the topological string partition function|^j (up to the change of variables, 
which we do not explicitly show here for simplicity): 

^bps - = ^top- (7-1) 

C-top 

It is natural to expect that similar story should exist for open BPS invariants 
as well. Namely, we expect to define open version of the BPS partition function EB 

7° 

depending on maps a, 8 specifying the chamber in the Kahlcr moduli space, such 
that the partition function reduces to the open topological string partition function 
in a special chamber Ct op : 



7° 

^BPS 



^tV (7-2) 



The question is how to define open BPS degeneracies such that the generating 
function follows the conditions above. 

As a guiding principle of our following argument, we use the crystal melting 
model developed in Section[3] As we have seen, this crystal melting model generalizes 
the result of Ref. [5S]for C 3 to an arbitrary toric Calabi-Yau manifold. In the case of 
C 3 , the crystal melting partition function with the boundary con ditions specified by 
three Young diagrams Ai, A2, A3 gives the topological vertex ^Sl C\ 1 \ 2 \ 3 . By using 
these vertices as a basic building block, we can compute open topological string 
partition function with non- com pact D-branes wrapping Lagrangian 3-cycles of the 
topology M 2 x S 1 included E31 j n this story, generalization from closed to open 

bm The upper index c stands for 'closed'. 
bn See I Appendix Al4] for details. 

bo Actually, the topological string partition function depends on the choice of the resolution of the 
singular Calabi-Yau manifold X. This is related to the choice of the chamber, as will be explained 
in the main text. 

bp Thc upper index o stands for 'open'. 
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topological string partition function corresponds to the change of the boundary 
condition for the crystal melting model for C 3 . 

Now the results presented in Section[3]show that the closed BPS partition func- 
tion discussed above can be written as a statistical mechanical partition function 
of the crystal model. This model applies to any toric Calabi-Yau manifold, and for 
C 3 the BPS partition function coincides with the topological string partition func- 
tion. Similarly to the case of the topological string story mentioned in the previous 
paragraph, we hope to define the open version of the BPS invariants by changing 
the boundary condition for the crystal melting model. The invariants defined in this 
way will be defined in any chamber in the Kahler moduli space, and reduces to the 
ordinary topological vertex in a special chamber. We call such a generalization of 
the topological vertex "the non-commutative topological vertex" following "the 
orbifold topological vertex" named in Ref. 11771 

We will see that this expect ation i s indeed true. We adopt the definition proposed 
in the mathematical literature I 83 I 84J although we hope to have clarified the presen- 
tation by using the language of the crystal melting model. Our non-commutative 
topological vertex is defined for a Calabi-Yau manifold X without compact 4-cycles, 
and a set of representations A assigned to external legs of the toric diagram. As in 
the case of topological vertex, A encodes the boundary condition of the D4-branes 
wrapping Lagrangian 3-cycles. We propose our vertex as the building block of open 
BPS invariants. 

We can provide several consistency checks of our proposal (see Section l7.3.4l for 
more details). First, our vertex by definition reduces to the closed BPS invariant 
when all the representations A are trivial. Second, our vertex shows a wall crossing 
phenomena as we change the closed string Kahler moduli, and the vertex coincides 
with the topological vertex computation in the chamber where the closed BPS 
partition function reduces to the closed topological string partition function. Third, 
the wall crossing factor is independent of the boundary conditions on D-branes, and 
is therefore the partition function factorizes into the closed string contribution and 
the open string contribution, as expected from Ref. 11781 and the generalization of 
the arguments of Section[5l 

Given a combinatorial definition of the new vertex, the next question is whether 
we can compute it, writing it in a closed expression. We show that the answer 
is affirmative, by proving the following statement. For a Calabi-Yau manifold X, 
the non-commutative topological vertex Cbps,(o-,0;A)(^O is equivalent to the closed 
BPS partition function Z^ ps ^, e ,s{X') for an orbifold X' of X , under a suitable 

b iThe word ' non-commutative' stem s from the mathematical terminologies such as "non- 
commutative crepant resolution" ESI and "non-commutative Donaldson-Thomas invariant" G2] 
The non-commutativity here refers to that of the path algebra of the quiver, and is different from 
the more familiar one (e.g. Refs. 11741 1175H in string theory contexts. The quiver (together with a 
superpotential) determines a quiver quantum mechanics, which is the low-energy effective theory 
on the D0/D2-branc worldvolumc (see Section[3]l- See also Ref. 1 1761 for "closed topological vertex", 
which correspond to C 3 /Z2 X Z2. 
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identification of variables explained in the following sections F*1 : 

CbPS,(ct,0;A)PQ = -ZbPS,((t',0')PO- (7-3) 

We will give an explicit algorithm to determine X' and a' ,9' , starting from the data 
on the open side. Since an infinite-product expression for Z^ ps g,JX') is already 
known 182153 this gives a closed infinite-product expression for our vertex. 

The rest of this section is organized as follows. We begin in Section lT^l bv slightly 
reformulating the results obtained in Section[5]the closed BPS invariants, their wall 
crossings, and their relation with the topological string theory. In Section [7T3l we 
define our new vertex using the crystal melting model. We also perform several 
consistency checks of our proposal. Section [7T4l contains our main result (|7.3p , which 
shows the equivalence of our new vertex with a closed BPS partition function under 
suitable parameter identifications. We also give an explicit algorithm for construct- 
ing the closed BPS partition function starting from our vertex. In Section[7AT] we 
treat several examples in order to illustrate our general results. We also include Ap- 
pendices |Appendix A.4| |Appcndix A. 5 and |Appcndix A. 6 for mathematical proofs 
and notations. 



7.2. Closed BPS Invariants Revisited 

Before discussing the open BPS invariants, we summarize in this section the defini- 
tion and the properties of the closed BPS invariants. This section hopefully serves 
as a self-explanatory introduction to the closed BPS invariants, and prior reading 
of Scction[5] is not necessary for the first reading of this section. 

Throughout this section, we concentrate on the case of the so-called generalized 
conifolds (sec Scction l5.4.2[) . The reason for this is that wall crossing phenomena is 
understood well only in cases without compact 4-cycles (see Section[5]), which means 
X is cither a generalized conifold or C 3 /(Z2 x Z2) 0- 

By a suitable SL(2, Z) transformation, we can assume that the toric diagram of 
a generalized conifold is a trapezoid with height 1, with length L + edge at the top 
and L- at the bottom (sec If we denote by L = L + + L_ the sum of the 

length of the edges on the top and the bottom of the trapezoid, this geometry has 
L — 1 independent compact P^s. We label them by cti, borrowing the language of 
the root lattice of algebra. 

The language of the root lattice will be used extensively throughout this section 
. We can also make more P^s by combining them. For example, combining all the 
P^s between i-th and j-th P 1 (assume i < j), we have another P 1 which we denote 

br More precisely, we need to specify the resolution of X and X' . We also need to impose the 
condition that two of the representations A are trivial. See the discussions in the main text. 
bs See Section r5.4.2l for the proof of this statement. 

bt The Calabi-Yau manifold is determined by L+ and L_ as xy = z L +w L - . 
bu The root lattice of Al-1 is exploited in Refs. l82l l57l l83l See also |Appendix A.4| 
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rTTT> 

Fig. 26. The toric diagram of a generalized conifold, with L_|_ = 3, L_ = 5. 

by 

ctjj := a, + . . . + otj. (7.4) 

This corresponds to a positive root of -Al_i. 

Suppose that we have a Calabi-Yau manifold X without compact 4-cycles. We 
also consider a single D6-brane filling the entire X and D0/D2-brancs wrapping 
compact holomorphic 0/2-cycles specified by n £ Hq(X;Z) and (3 € H2(X;Z), 
respectively. We can then define the BPS degeneracy f2(ro, j3) counting BPS degen- 
eracy of D-brancs 0. The closed BPS partition function is then defined by (see 
1231) 

Z c BPS (q 1 Q) = Y,n{n,P)q n Q fi . (7.5) 

The closed BPS partition function for generalized conifolds is studied in Refs.l82l 
l57l To describe the results, let us first specify the resolution (crcpant resolution) of 



X Q Each of the L - 1 P^s is either an 0(-l, -l)-curve or an 0{-2, 0)-curve. 



In the language of the toric diagram, this choice is specified by the triangulation of 
the toric diagram. We specify this choice by a map 

a : {1/2, 3/2,..., L- 1/2} -> {±1}. (7.6) 

In the following we sometimes write ± instead of ±1. When a(i — 1/2) = 1 (a(i — 
1/2) = —1), the i-th triangle from the left has one of its edges on the top (bottom) 
edge of the trapezoid. This means that the z-th P 1 is a 0(— 1, — l)-curve (0(—2, 0)- 
curve) when a(i - 1/2) = -a(i + 1/2) (a(i - 1/2) = a(i + 1/2)). By definition, we 
have lo-- 1 ^ 1 )! = L±. 

For example, in the case of Suspended Pinched Point (L + = 1,L_ = 2) whose 
toric diagram is shown in Fig.[27l L = 3 and there are 3 difference choice of resolu- 
tions. This is represented by 



C3 



{1/2, 3/2, 5/2}^ {-,-,+}, 

{1/2, 3/2, 5/2}^ {-,+,-}, (7.7) 
{1/2,3/2,5/2}^ {+,-,-}. 



More precisely, this BPS degeneracy is defined by the second helicity supertrace. 
bw This is not essential, since by varying the value of the Kahler moduli we can go to the geometry 
with other choices of resolutions. We just need to specify an arbitrary resolution in order to begin 
the discussion. See | Appendix A~4] for more on this. 
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Fig. 27. The choice of resolutions of a generalized conifold (L+ = 1, L- = 2). 



Given cr, the topological string partition function is given by 



15611611 



top,(T 



•,Q 



Ull(i-q n Qy 

n=l /3 



(7.8) 



where iVf is the genus Gopakumar-Vafa invariant 0. For the 2-cycle /3 
Q4 + . . . + oij , the explicit form of Ng depends on a and is given by (see (|5.32l) ) 



N: 



g=0 

' = Oti-\-...OLj 



(7.9) 



1 "p=a z + ...a j - 



(_]_\l+lt{fc|i<fc<j, <r(A:-l/2)#(T(fc+l/2)} 
/ ^ l-\-${k\i<k<j, is a 0{ — 1,-1)— curve} 

By CPT invariance in five dimensions we have N 9 a ~° , , > 
We also have iVS_ = x{X), where the Euler character x(X) for a toric Calabi-Yau 
manifold is the same as twice the area of the toric diagram. 

As shown in Section^ the closed BPS partition function is given by 

Zbps(<7,<2) = ^top('7,Q)Ztop('7,Q" 1 )| chambcr 

II (l - q n Q f) ) nN K (7.10) 
(P,n):Z(p,n)>0 

where the central charge Z(j3, n) is given in (|5.8j) . There 1/ R denotes (up to propor- 
tionality constants) the central charge of the DO brane, and following the discussion 
in Section l5T2l we choose the complexified Kahler moduli to be real. 

Now suppose that 1/R is positive 0. From (|7.10|) and (|5.8[) . it follows that the 
wall crossing occurs when the integer part of the value of the B-ficld through the 
cycle change. For the cycle on + . . . + oij, this is given by 

[S(ai) + ... + B( aj )], (7.11) 

Since there are L — 1 P^s in X, there are L(L — l)/2 such parameters. 

We can take a special limit B{a,i) —> oo for all i. Let us denote this special 
chamber by C'top- As discussed at the end of Section RTBI in this limit the BPS 
partition function reduces to the closed topological string partition function: 



C'top 



= Z. 



topi 



(7.12) 



Higher genus Gopakumar-Vafa invariants vanish for generalized conifolds. 
by Under this condition we are discussing only half of chambers of the Kahler moduli space, which lie 
between the Donaldson-Thomas chamber and the non-commutative Donaldson-Thomas chamber. 
The other half arises when 1/R is negative. 
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just as advertised in (|7.1|) . 

For concreteness, let us discuss the example of the Suspended Pinched Point 
(L = 3) using the triangulation a\ in (|7.7j) . In this example, the topological string 
partition function is 

oo oo oo 

Ztop^fo.Q) = u{ q fi 2 JJ(i - g B Qi)-» JJ(i - <f Q 2 )" IJC 1 - <TQM n > 

n—1 n—1 n—1 

(7.13) 

where M(g) is the MacMahon function defined previously in (|3.32[) .The BPS par- 
tition function is given by 

oo oo oo 

Z B P S (q, Q) = M{qf JJ (1 - q n Q 1 )- n JJ(1- <f Q 2 )" J] (1 q n QiQi) n 

n—1 n—1 n—1 

oo oo oo 

x n (i-^Qr'r™ n (i-w 1 )" II (i-™^)- 1 ) 

n>[B( ai )] n>[B(a 2 )] n>[B(ai+Q 2 )] 

(7.14) 



The parameters [B(cti + . . . + ctj)] specify the chamber, but as we can see from 
the definition they are not completely independent parameters. Since we only have 
L — 1 real parameters it is likely that this parametrization is redundant. Indeed, 
as explained in the | Appendix AT4| we can specify the chamber by a map 9, which is 
specified by L half-integers, 6(1/2), 6(3/2), . . . , 9(L — 1/2), satisfying one constraint 

x>H)=£H)- <" 5 » 

i=l v ' i=l v 7 

This means we can indeed parametrize the chamber by L — 1 independent (half- 
integers, which is what we expected. As discussed in |Appendix A. 4 9 is an element 
of the Weyl group of Al-\- 



7.3. The Noncommutative Topological Vertex 

In this section we give a general definition of the non-commutative topological vertex 
using the crystal melting model. This definition is equivalent to the one given in 
Ref. [83] using the dimer model 0- See Ref. [84] for more conceptual definition in 
terms of Bridgeland's stability conditions and moduli spaces. 
To define our vertex, we need the following set of data: 

• A map 

a:{l/2,...,L-l/2}^{±}. 



See | Appendix A.l| for the equivalence between crystal melting model and the dimer model. 
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As already explained in Section [7T2l this gives a triangulation of the toric 
diagram, or cquivalently the choice of the resolution of the Calabi-Yau 
manifold. 

A map 9 : Z/j — > Z/,, where Z^ is the set of half-integers. As explained 
in | Appendix A. 4 in the case of closed BPS invariants, 9 and a specify the 
chamber structure of the open BPS invariants. 

A set of Young diagrams A, assigned to external legs of the (p, q)-web. This 
specifies the boundary condition of the non-compact D-brancs ending on 
the (p, g)-web. We denote by Ai,...,Ai the Young diagrams for the top 
and the bottom edges of the trapezoid, and by A+,A_ the remaining two. 
We sometimes write A = (/i, v), where /i = (fix, . . . = (Ax, . . . , A^) and 
v = (A+, A_). In the example shown in Fig.[28l there are 5+2 external legs 
and therefore 5+2 representations. 



A2 A5 




Ai A3 A4 



Fig. 28. Representations assigned to external legs of the (p, g)-web. The dotted lines represent 
the (p, q)-web. 



For later purposes, we combine h\,...,hl into a single representation /i by 

fi(i - 1/2 + kL) :=Hi(k- 1/2). (7.16) 

In other words, we choose fi such that L-quotients of fi give fix, . . . , /if,. By 
abuse of notation, we use the same symbol /x for a set of representations 
fix, ■ ■ ■ , as well as a single representation define above. 

Given er, 9 and A, we define the non-commutative topological vertex 

CBPS,(<T,6(;A)(<7: Q). 

In the following we drop the subscript BPS for simplicity. 

Before going into the general definition, we first illustrate our idea using simple 
example of the resolved conifold. 



7.3.1. Example: Resolved Conifold 

In this example there is only one P 1 and the BPS partition function depends on a 
single positive integer TV := [-B(ai)]. In the language of 9, 

(9(1/2) = 1/2 -TV, 9(3/2) = 3/2 + TV. (7.17) 
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Wc fix a to be 

a(l/2) = +, a(3/2) = -. (7.18) 

Without losing generality we concentrate on N > 0, since N < corresponds to a 
flopped geometry, where a is replaced by — cr (see | Appendix A^4| . 

The ground state crystal for N = 2 is shown in Fig. [29] (this is the same figure as 
Fig. 1231 reproduced here for the convenience of the reader). This crystal, sometimes 
called a pyrami d, cons ists of infinite layers of atoms, the color alternating between 
black and white ESlSOl j n the N-th chamber there are A^ + 1 atoms on the top. 




Fig. 29. The ground state crystal for the resolved conifold for N = 2. The crystal consists of an 
infinite number of layers, and only a finite number is shown here. The ridges of the pyramid arc 
represented by four lines extending to infinity. This is the same figure as in Fig. 1231 



The closed BPS partition function is defined by removing a finite set of atoms 
f2 from the crystal. When we do this, we follow the melting rule (|3.16[) such that 
whenever an atom is removed from the crystal, we remove all the atom above it. 
We then define the partition function by summing over such f2: 

z = Y / ^o N) ) Wo{n) (q[ N) ) wim , (7-19) 

o 

where wo(fl) and u>i(f2) are the number of white and black atoms in f2, respectively 
The weights q$ N ^ (q[ N ^) assigned to white (black) atoms in the N-th chamber 

are determined as follows. We can slice the crystal by the plane, and each slice 

is specified by an integer i (see Fig. [3UJ) . Wc choose i so that the slice i = cuts 

through the center of the crystal. 

The weight q^ N ^ depends on the chamber and is given by 

<tf°=G N G N+ \ q[ N) =$ +1 1?, (7-20) 

when N is odd, and 

= a^V, i[ N) -% N ii N+ \ (7-21) 
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Fig. 30. We can slice of the conifold crystal by an infinite number of parallel planes. 

when N is even. This coincides with the expression in (|5.42[) . up to the exchange 
of two nodes of the quiver when N is odd. For example, = q$, q[ ^ = qi when 
N = 0, and q^ = q^ 1 , = qfiqi- The change of variables arises from the Seiberg 
duality on the quiver quantum mecha nics^H geometrically mutations in the derived 
category of coherent sheaves l &6 | 54 | ( se e Section l3.3.4j) . or in more combinatorial 
language the dimer shuffling ^01. The parameters qo , q\ defined here are related to 
the D0/D2 chemical potentials introduced in Section lT2l by (see (|3.2f |) ) F*1 

q = qoqi, Q = qi- (7.23) 

Now let us discuss the open case. When non-trivial representations arc assigned 
to each of the four external legs of the (p, </)-web, the only thing we need to do is 
to change the ground state of the crystal. 

The crystal has four ridges, corresponding to four external legs of the (p, q)-web. 
When we assign a representation, we remove the atoms with the shape of the Young 
diagram. More precisely, we remove the atoms with the shape of the Young diagram 
in the asymptotic direction of the (p, q)-web, as well as all the atoms above them, 
so that the melting rule is satisfied. See Fig.|3f]for an example. 

The partition function is defined in exactly the same way by (|7.19j) . and the 
result is denoted by Cr a ,e-,X)- 

Several comments are now in order. 

First, let us explain the origin of the name "the non-commutative topological 
verted' . Recall that, in commutative case, topological vertex is defined for C 3 . For 

ca The equation 117. 2111 is the same for N odd and even if we suitably exchange the two nodes of 
the quiver diagram. The relation l|7.23|l can also be written as 

q = q ™ q ™, Q = (gi N Y(*l N Y + \ (7.22) 

when N is even, and g^ N ^ and q± exchanged when JV odd. This coincides with the expression in 
Ref.lMl 
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Fig. 31. The pyramid for open BPS invariants. A non-trivial representation (1, 1) is placed on 
with one of the four external lines. As compared with the previous figure, atoms colored gray, 
corresponding to the Young diagram, are removed from the crystal. The red atoms have no atoms 
above them. 

a general affine toric Calabi-Yau manifold X, we divide the polygon into trian- 
gles and assign a topological vertex to each trivalent vertex of the dual graph. We 
can get the topological string partition function for the smooth toric Calabi-Yau 
manifold Y (resolution of X) by gluing the topological vertices with propagators. 
Similarly, assume that a polygon is divided into trapezoids. Then we can assign 
a non-commutative topological vertex to each vertex of the dual graph and glue 
them by propagators. The BPS partition function defined in this way is related to 
the topological string partition function via wall-crossing F^l. In Refs. 11771 11791 thev 
study the case when a polygon is divided into (not necessary minimal) triangles. 

Second, it is possible to give more geometric definition of the vertex (see Ref.l84|). 
For the closed BPS invariants, the crystal arises as a torus fixed point of the moduli 
space of the modules of the path algebra quiver (under suitable ^-stability con- 
ditions). The moduli space is the vacuum moduli space of the quiver quantum 
mechanics arising as the low-energy effective theory of D-branes (SectionE} . The 
similar story exists in our case. Namely, the crystal is in one-to-one correspondence 
with the fixed point of the moduli space arising from a quiver diagram. For example, 
for conifold with A =B, the quiver is given in Fig. [32] 

Third, in the case of C 3 , our vertex reproduces the topological vertex for C 3 by 
definition. 



Given a division of a polygon into trapezoids, we get a partial resolution of X and a non- 
commutative algebra A over the partial resolution, which is derived equivalent to Y . The BPS 
partition function given by gluing non-commutative topological vertices counts torus invariants 
A-modulcs. 
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Fig. 32. Quiver diagram for the open invariant with A = (1, 1). This is the Klebanov-Witten quiver 
1 10" I with an extra node and extra three arrows starting from it. The three arrows correspond to 
three red atoms in Fig. 1311 Note that this is different from the quiver in Fig. [8] 



7.3.2. General Definition from Crystal Melting 

We next give a general definition of the vertex. Readers not interested in the details 
of the definition of the non-commutative topological vertex can skip this section on 
first reading. 

Given a boundary condition specified by a, 6 and A = (/i, v), we would like to 
construct a ground state of the crystal, and determine the weights assigned to the 
atoms of the crystal. 

The basic idea is the same as in the conifold example. First, the closed string 
BPS partition function is equivalent to the statistical partition function of crystal 
melting. The ground state crystal can be sliced by an infinite number of parallel 
planes parametrized by integers n G Z, just as in Fig.[30l On each slice, there are 
infinitely many atoms, labeled by integers (x, y) G Z> . Therefore, the atoms in the 
crystal are label by (n,x,y) G Z x Z> . 

Let us show this in the example of the Suspended Pinched Point. The crystal in 
Fig. E31 clearly shows this structure. 

Another way of explaining this is to construct a crystal starting from a bipartite 
graph on K 2 , shown in Fig. [34] El In this example, the bipartite graph consists 
of hexagons and squares, and periodically changes its shape along the horizontal 
directions. 

The atoms of the crystal are located at the centers of the faces of the bipartite 
graph, and it thus follows we can slice the crystal along the horizontal axis. Each slice 
consist of an infinite number of atoms labeled by two integers (x, y) G (Z>o) 2 , since 
there are two directions, the horizontal direction and the perpendicular direction to 
the paper Q 

Now consider the open case. In this case, we construct a new ground state 
by removing atoms from the closed ground state. By the melting rule, the atoms 
removed from the n-th plane should be labeled by (x, y) G V(n), where V(n) is a 



cc This is a universal cover of the bipartite graph on T 2 , which appears in the study of four- 
dimensional Af = 1 quiver gauge theories. See Section|3]2] 

cd In general the bipartite graph is determined by a. A hexagon (a square) corresponds to 0(— 2, 0)- 
curve (0(— 1, — l)-curve). In other words, the i-th polygon is a hexagon (square) if o~(i) = a{i + 
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Fig. 33. The crystal for the Suspended Pinched Point. We can slice the crystal along planes 
represented by lines, which come with three different colors. This crystal appeared previously in 
Fig.[l2] 




__ + _- +__+__ 



Fig. 34. The bipartite graph for Suspended Pinched Point. We here take a = <j\ and 8 = id. The 
red undotted (the blue dotted) lines have half-integer (integer) values of the coordinate along the 
horizontal axis. This figure is the universal cover the bipartite graph on T 2 , discussed previously 
in Fig. [5] and in Fig|6] 

Young diagram. Depending on the representations on external legs, V(n) increases 
or decreases as we change n. Thus the ground state crystal for open BPS invariants 
are determined by such a sequence of Young diagrams {V(n)}, called transitions 
below. In the following we make this idea more rigorous. 

Let us begin with some notations. Let p, and \i! be two Young diagrams. We say 
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+ 



/ if the row lengths satisfy 



/i >- fi' 



Mi > Mi > M2 > M2 > ■ ■ ■ > 



(7.24) 



and jjt y fx' 



1 if the column lengths satisfy 



Vi > Vi > V2 > V2 > • • ■ • 



(7.25) 



We define a transition V of Young diagrams of type (cr, 9] n,v) as a map from the 
set of integers Z to the set of Young diagrams such that 

• V(n) = V- for n <C and V(n) = v + for n ^> 0, 

ao6(h) 

• V(h — fxo 9{h)/2) y V(h + fio 9(h)/ 2) for h: half-integer. 

Then as shown in Ref . |531 there is a minimal transition V m i n of Young diagrams of 
type (cr, 9] fi, v) such that for any transition V of Young diagrams of type (cr, 9; /1, v) 
and for any n £ Z we have V(n) D V m i n (n). 

For a transition V of Young diagram of type (a, 9] fi,is), the corresponding crystal 
configuration can be defined by 



where a(n,x,y) denotes the atom at position (n,x,y). In particular, ^4(V m i n ) gives 
the ground state crystal. 

Having defined the ground state crystal, the partition functions is defined again 
as the sum over a subset Q of A(V m m) satisfying the following two conditions: l cc l 

• £1 is finite set, and 

• ft satisfies the melting rule (|3.16j) . In other words, if a' £ f2 and a' = aa 



For a crystal ft C A(V m i n ), we define the weight w(fl)i by the number of atoms 
with the color i contained in fl: 



ce It is straightforward to show that a subset Q C A(V m i n ) satisfies the two conditions if and only 
if 7i(V m i n )\f2 = A(V) for a transition V of Young diagram of type (cr, 9; /i, v). 

We can also explicitly write down the melting rule using the coordinates (n,x,y). Let us write 
a □ a' when there is a path (a composition of arrows) a such that a' = aa. The partial order Zl 
is then generated by 



A(V) := {a(n,x,y) n £ Z, x, y £ Z> , (x,y) £ V(n)}. 



(7.26) 




w(Q)i :— j}{a(n, x, y) £ fl \ n = i (modL)}. 



(7.27) 



a(h - 1/2, x, y) Zl a(h + 1/2, x, y) 
a(h + 1/2, x, y) □ a(h - 1/2, x, y) 



H 06(h) = + 
pa 0(h) = - 



a(h + 1/2, x, y) □ a(h - 1/2, x + 1, y) 
a(h + 1/2, x, y) □ a(h - 1/2, x, y + 1) 
a(h - 1/2, x, y) □ a(h + 1/2, x + 1, 3/) 
a(h - 1/2, x, j/) Z) a(/i + 1/2, s, y + 1) 
a(n, x, y) Z) a(n, x + 1, y) 



ao8(n- 1/2) = a o 6(n + 1/2) = + 
ao9(n- 1/2) = a o 6(n + 1/2) = - 



ixo 6(h) = +,ao 0(h) = - 
/x o 6»(/i) = +,cro 0(h) = + 
(io 0(h) = — , cr o 0(h) = - 
(io 0(h) = (to 0(h) = + 



a(n, x, y) □ a(n, x, y + 1) 
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Also, for 9, we put 

Qi : = <?e- 1 ( l -i/2)+i/2 • ( 7e- 1 (i-i/2)+3/2 <76»- 1 (i+i/2)-i/2 (7.28) 
when O-^i - 1/2) < fl" 1 ^ + 1/2), and 

1i '■= ?0-i(i_l/2)-l/2 ' ^fl" 1 ^- 1/2) -3/2 9 e -l( i+ i/2) + l/2 (7.29) 

when — 1/2) > £> _1 (i + 1/2). Here we defined qi for i 6 Z periodically, 

qi+L = Qi- (7.30) 

We then define the vertex by 

C\te^){q^ ■ • • , 9L-i) := E(9o) lu(n) ° (d-i) 1 "^- 1 - (7-31) 

n 

The parameters qo, . . . , gz,-i defined here are related to the D0/D2 chemical poten- 
tials introduced in Section l7T2l by 

g = go---?i-i, Qi = qi (i = l,...,L- 1). (7.32) 

7.3.3. Refinement 

We can also generalize the definition to include the open refined BPS invariants EE 

Let us start by recalling the meaning of the refined BPS counting, first in the 
closed case ^H As we have seen in Section[5l when the type IIA brane co nfigura tion 
is lifted to M-theory anc i w hen we use the 4d/5d correspondence 1 5 8 1 5 9 1 ^ ^ e 
D0/D2-branes are lifted to spinning M2-branes in R 5 , which has spin under the 
little group in 5d, namely SO (A) = SU(2) L x SU(2) R . The ordinary BPS invariant 
is defined as an index (|5.5[) : it keeps only the SU{2)l spin, while taking an alternate 
sum over the SU(2)r spin. The refined closed BPS invariants is defined by taking 
both spins into account. 

The situation changes slightly when we consider open refined BPS invariants. 
The D4-branes wrapping Lagrangians, when included, are mapped to M5-branes on 
R 3 . This means that 5*0(4) is broken to 50(2), and we have only one spin. However, 
there is an 50(2) R-symmetry for M = 2 supersymmctry in three dimensions, 
and in the definition of the ordinary open BPS invariants we kee p on ly one linear 
combination of the two, while tracing out the other combination LL°^I The refined 
open BPS invariants studied here takes both of the two charges into account. 

In the language of crystal melting used in this paper, the open refined BPS 
invariants are defined simply by modifying the definition of the weights. Here, we 
explain how to modify the weights in the case of fi = 0. 

For an integer n define 

w(CL) n :=%{a(n,x,y) eft}. (7.33) 



c sSee Refs. HSEZHHO] for c i ose d re fined BPS invariants. 
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We also define the weights by 

In '■= <7e- 1 (n-l/2) + l/2 ' 90" 1 (n- l/2)+3/2 ^0- 1 (n+1/2) - 1/2 (7.34) 

when e-^n - 1/2) < fl-^n + 1/2), and 

3n := '?e- 1 i(n-l/2)-l/2 ' ( n _ 1/2) -3/2 ^- 1 (n+l/2)+l/2 (7.35) 

when 0- 1 (n - 1/2) > -1 (ra + 1/2). Here we defined 

<?„ := g„ (7.36) 

when n ^ mod L (recall g„ is understood to be periodic with respect to n (|7.30p ). 
and 

g+ n > 

(g+g.) 1 / 2 n = 
g_ n < 

when n = mod L. We then define the refined vertex by 

C i<Ji9 . Av) (q + ,q.,q 1 ...,q L - 1 ) := £ JJ (Cr (n) "- (7-38) 

By definition, the refined vertex reduces to the unrefined vertex by setting q + = 
g_ = go- The reader can refer to Ref. [S3] for the definition of weights in general 
cases 



7.3.4. Consistency Checks of Our Proposal 

In |Appendix A.6| we gave a purely combinatorial definition of the non-commutative 
topological vertex. We now claim that is captures open BPS invariants in the fol- 
lowing sense. 

Consider a generalized conifold (a toric Calabi-Yau manifold without compact 
4-cycles) with representations assigned to each leg of the (p, g)-web. Each repre- 
sentation specifies a boundary conditio n on the non-compact D4-brane wrapping 
Lagrangian 3-cycle of topology M 2 x S 1 

In the absence of D4-branes, we are counting particles of D0/D2-branes wrap- 
ping 0/2-cycles, which makes a bound state with a single D6-brane filling the entire 
Calabi-Yau manifold. When the D4-branes are included, D2-branes can wrap disks 
ending on the worldvolumc of D4-branes. The degeneracy of such D-branc configu- 
rations is what we mean by the open BPS degeneracies. Note that supersymmetry 
is broken by half due to the inclusion of D4-branes; our counting of BPS particles 
makes sense because we are counting BPS states in lower dimensions, where the 
minimal amount of supersymmetry is lower. 

We can provide several consistency checks of our proposal. First, our vertex by 
definition reduces to the closed BPS invariant when all the representations A are 
trivial: 



CbPS,(<t,0;A=0) - ZbPS,(o-,0)- 



(7.39) 
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The second consistency check comes from the wall crossing phenomena. As shown 
in Rcf. |83j the vertex goes through a series of wall crossings as we move around the 
closed string moduli space (Kahlcr moduli of the Calabi-Yau manifold), just as in 
the case of closed invariants. It was also shown in Ref. [83] that in the chamber Ct op 
where the closed BPS partition function reduces to the closed topological string 
partition function, our vertex gives the same answer as that computed from the 
topological vertex (in the standard framing): 



C 



BPS,(cr,0;A) 



^topological vertex, A ■ (7.40) 

Ctop 



The third consistency check comes from the fact that the wall crossing factor is 
independent of representations. In other words, 

^ / \ C(t7,0;/i,i/)(?O! • • • > 1l-x) ln ,-.<. 

L(o,6;w){qo,---,qL-i) j; j-g — e y (7.41) 

does not depend on 8 I83 | 84 | ^ gee a i g0 Section !?. 4. 21 for an explanation) 0. This 
means that the open BPS partition function, which is defined by the sum over 
representations, takes a factorized form 

ryO _ ^BPS ryp _ ^BPS 7 1) C7AO\ 

^BPS — ~~7~Z ^top — ^top- 
^top -^top 

Since Z^ PS /Z£ op takes an infinite product form as explained in Section !7.2l and Z° op 
also takes the infinite product form HlHl Zgpg itself should take an infinite product 
form, which is consistent with a suitable generalization of the argument of Section[5] 
We will see discuss this in more detail in Section [T31 

Using (|7.42p . we can compute our vertex using the ordinary topological vertex 
formalism. In the next section, we give a yet another way of computing the non- 
commutative topological vertex. The advantage of our approach is that the final 
expression manifestly takes a simple infinite product form, and we do not have to 
worry about the summation of Schur functions. 



7.4. The Closed Expression for the Vertex 

In this section we give a closed expression for our non-commutative topological 
vertex. We do this by proving a curious identity stating the equivalence of our 
vertex for a toric Calabi-Yau manifold X with a closed BPS partition function for 
an orbifold of X 0. We also comment on another method via vertex operators. 



ch To be exact, we have to normalize the generating function by a monomial. See Ref. 1841 Corollary 
3.21, for the precise statement. 

"This is reminiscent of story of Ref. ll83l where the 'bubbling geometry' X' is constructed for given 
a toric Calabi-Yau manifold X such that the open+closed topological string partition function on X 
is equivalent to closed topological string partition function on X' . However, our story is different 
in that the vertex computes only a part of the full open BPS partition function; the partition 
function itself is given by summation of our vertices over representations. 
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Start from a non-commutative topological vertex for a generalized conifold X, 
which has L — 1 compact P^s. As we discussed above, for the definition of the 
vertex we need (1) a for a choice of the crepant resolution of X, (2) a map 9 
specifying the chamber together with a (3) a set of representations A = (jj,, u), and 
the resulting vertex is denoted by C cr ^-\(q 7 Q)- In the following we consider the 
special case v = (A+, A_) = (0, 0). 

Choose an integer M and consider the Zm orbifold X' of X. We choose the 
orbifold action such that when the toric diagram of X is a trapezoid with a top and 
the bottom edge of length L + and L_ respectively, X' has length ML + and ML_. 
We also choose map 

a' : {1/2, 3/2, . ..,ML— 1/2} -> {±1}. (7.43) 

and 

0' : Z h -> Z h , 9'(h + ML) = 6'(h) + ML (7.44) 

such that 

ao9 = a'oQ', /io0 = 0o0 / . (7.45) 

Then 

C{a,e$,n){q,Q) = c '(ff',fl'i0,0)(9 / j < 9')l 9 «= g '«' =«'?+ i =-9'?+ CM _i )i ' ( 7 ' 46 ) 

where £ = 0, — 1. 

See |Appcndix A.6] for an explicit method for choosing such M,cr',8' satisfying 
(|7.45|) as well as generalization of (|7.46j) to the case of refined BPS invariants. 

Since an infinite-product expression for closed BPS partition function for a gen- 
eralized conifold is already known (Section l7.2[) , we have a closed expression of our 
vertex when v = 0. 

7.4.1. Examples 

Let us illustrate the above procedure by several examples. 
Example: C 3 

First, we begin with the non-commutative topological vertex for C 3 . Since there 
is no wall crossing phenomena involved in this case EH, the vertex should coincide 
with the ordinary topological vertex, thus providing a useful consistency check of 
our proposal. For C 3 , we have L = 1, er(l/2) = —1, 6 = id. 

Take C 3 with representation A = (fj,, v = 0) with \i = (N, N — 1, . . . , 1) at one 
leg. The method in |Appendix~Q1 gives M = 2 (therefore X' = C 2 /Z 2 x C), and 

9' (1/2) = 1/2 -N, 6' (3/2) = 3/2 + N, a' = -1, (7.47) 

CJ There is one wall between R > and R < 0, however we do not discuss such a wall since we are 
specializing to the case R > 0. 
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and thus 

[-B(ai)] = N. (7.48) 
The weight is given by (|7.20[) or (|7.21[) . By solving for 

q't = <?'?' = q, (7.49) 

we have 

q'o = q- 2N+ \ q'l=q 2N+1 (7.50) 
in the case of N odd. Substituting this into the closed BPS partition function 

- 9'oViT 2n ?'o</i +1 )-" II (! - (7-51) 

n>0 n>0 n>iV 

we have 

iV 

M(q)/Y[(l-q 2n ~ 1 ) (N+1) ~ n - (7.52) 

n=l 

This coincides with the know expression for the topological vertex ^10- The case 
of N even in similar. 

Example: Resolved Conifold 

Now let us discuss the next simplest example, the resolved conifold. 
Consider the representation /i = (EP, □), with 

6 = id, cr(l/2) = +, cr(3/2) = -. (7.53) 

In this case, the method in |Appcndix A. 6| gives 

6>'(l/2) = —7/2, 0'(3/2) = -1/2, 0' (5/2) = 11/2, 6' (7/2) = 13/2 (7.54) 

with 

(j'(l/2) = +, a'(3/2) = - </(5/2) = - ^(7/2) = +. (7.55) 
Then we have 

[B(ai)\ = [B(a 3 )} = 0, [B(a 2 )} = [B(a 2 + a 3 )] = 1, 



[B(ai + a 2 )\ = [B(ai + a 2 + a 3 )} = 2, 

and 

q'o = % 3( ii 3 = <1~ 3 ' 1i = <ioql = qQ, 

q 2 = qlql = q 3 -. q'z = qo = qQ~ 1 - 

t 

ck In the normalization of Ref. 1561 coincides with M(q) <?~" M "/ 2 C^^^ 



(7.56) 



(7.57) 
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— L\ — n 



The closed BPS partition function corresponding to the B-field (|7.56p is the follow 
ing: 

n (i - <?"t 4 " n c 1 ? m ^) n n o- - ^q^t n ^ <™r" n (* ^qv) 

n>0 n>0 n>0 n>0 n>l 

n a n a ^q's i ) n n (! - ^i^r n (* wr 

n>0 n>0 ra>0 ri>2 

n (i - v' n Q2Q3T n ( x - ^"Va" i ) n 

n>0 n>l 
n>0 n>2 

(7.58) 

where as in (|7.32j) we defined 

q' = qWiqWa, Q'i = Qi (» = 1,2,3). (7.59) 

Substituting (|7.57p for (|7.58p under the identification (|7. 32[) , we obtain the open 
BPS partition function: 

^ncdt(?, Q) ■ (1 - 9)" 3 (1 - ^T^l - Q) 2 (l - 9 2 Q)(1 - q 2 Q- V ), (7-60) 

where 

^ C dt(9, q) ■■ = n (! - « n r 2 " n c 1 - i n ^ n n (! - ^g -1 )" ( ? - 61 ) 

n>0 n>0 n>0 

= z t t op (g,g)z t c op (g,g- 1 ). (7.62) 

This coincides with the expression computed from the result of Ref. 11611 

M{qf J] (1 g"Q) fl - )3 fl _ 3) ( X - Q)^ 1 - ^ 2 )(1 - ^~ 2 ) (7.63) 

up to the wall crossing factor and the normalization by a monomial as remarked in 
footnote [ch] 

7.4.2. Expression via Vertex Operators 

As we have seen the computation of our vertex from the relation (|7.46[) is always 
possible but time-consuming in general. One efficient method is to use the vertex 
operator formalism, explained in Section ^. 61 By straightforward generalization of 
the argument there we can write down a rather general formula for the vertex, for all 
chambers and for all toric Calabi-Yau manifold without compact 4-cycles. Indeed, 
from the definition of the crystal it follows that E3 

C(„, 9 .,r,u)(qo, • • • , Ql-i) = (A+| . . . 0(-l/2)0(l/2)0(3/2) . . . |A_>, (7.64) 

cl To be precise, we have to divide this result by an overall monomial when either A_|_ or A_ is 
non-trivial. See Proposition 3.15 of Ref. 1841 
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where 

0(h) :=r^::gj( ge - 1(h) ). 

The fact that (|7.4ip is independent of the closed string moduli (i.e. the map 6) 
can be verified from this expression. This result shows that wall crossing can be 
understood as the result of commuting vertex operators. We do not discuss the 
details here and interested readers are referred to Ref. H3J 



7.5. Open BPS Wall Crossing and M-theory 

In previous sections we discussed wall crossings of open BPS invariants with respect 
to the closed string moduli. However, we expect that there should be wall crossing 
associated with the open string moduli as well. Therefore, the question arises: At 
which values of the open string moduli is our vertex defined? How does the vertex 
change as we change the o pen string moduli? The answers to these question are dis- 



cussed in the recent paper 1^21 by generalizing the argument of Ref. [57] discussed 
in Section[5lF"l. To explain this is the topic of this section. 

Our main result in this section is stated as follows. Let us denote by Zbps the 
BPS partition function including all BPS particles, open and closed. Then M-theory 
lift predict that the chamber dependence of this is captured by the open+closed 
partition function as 

Z£l{q,Q 1 v) = Z°^{q,Q 1 v)Z°+;-*(q,Q~\v- 1 ) , (7.65) 

chamber 

where Z°^'*(q, whose closed string part coincides with Z£ op (q, Q" 1 ), 

will l)e defined shortly. 

We first begin in Section [7.5.1l with t he review the open topological string parti- 
tion function, focusing on the integrality anc [ the connection with counting 
M2-branes ending on M5-branes. We also generalize the work of Ref. to open 
topological string and open BPS bound states. Based on these results we derive 
in section 17.5.21 the wall crossing of the open BPS partition function. In the next 
section fSection [7.6[) we discuss the relation to the work of Cecotti and Vafa. 

As before, consider adding M D4 brancs wrapping a special Lagrangian subman- 
ifold L of X. The D4 branes fill R 1 ' 1 subspacc of the flat space, and break half the 
supersymmctry of the Calabi-Yau. The theory on the branes thus has Af = (2, 2) 
supersymmctry in two dimensions. We now get new kinds of BPS pa rticle s, cor- 
responding to D2 branes wrapping disks and ending on the D4 brane [ESI These 
particles can form BPS bound states with the D6 brane and the closed D2 branes 

cm Rclatcd results on open BPS invariants have been obtained independently in Ref. 1 1841 We thank 
R. Dijkgraaf, P. Sulkowski and C. Vafa for sharing this information with us. 

cn Chamber dependence of open BPS states, in the context of surface operators of Af = 2, d = 4 
gauge theories, was discussed recently in Ref. 11851 See also Refs. 186, 90 for related discussions. 
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and the DO branes we had before, and pin them to the D4 branes. These are the 
open BPS invariants that we would like to count 0. 

A D2 brane ending on a D4 brane is magnetically charged under 5d gauge field 
A on the D4 brane, and electrically charged under the corresponding dual two-form 
B, dA = *dB. Suppose b±(L) = r. We get r 1— cycles on L that are contractiblc 
in X and fill in to holomorphic disks in X. We will assume, for simplicity that 
r = 1, tough the generalization to arbitrary r is manifest. Integrating the B-field on 
the S 1 we get a U(1) M magnetic gauge group generated by f s± B. The D2 branes 
wrapping the holomorphic disks and ending on the D4 brane are BPS particles 
charges under the magnetic U(1) M . The U(1) M gauge fields in two dimensions sit 
in the twisted chiral multiplets £j, i = 1, . . . M , whose lowest components Ui enter 
the BPS masses of the particles charged under thcrrFp|. Since we are not interested 
in the gauge dynamics on the D4 branes, we will view the gauge symmetry as a 
global symmetry, and as background multiplets. The charges of a D2 brane under 
U{1) M keep track of which D4 branes the D2 brane has boundaries on, and how 
many times it wraps the corresponding S . More precisely, since the M D4 branes 
are identical, we have Sm permutation invariance so the particles are representations 
R of U(1) M /Sm (which can also be viewed as representations of U(M)). 

When we lift this to M-thcory on X x Taub-NUT x R, we now get an M5 brane 
wrapping the Lagrangian L and filling R 2,1 , where open D2 branes become M2 
branes ending on the M5 brane. On the M5 branes, there is a U(1) M gauge theory, 
with N = 2 supersymmetry, where the open M2 branes are BPS particles. Under 
the same assumptions as in the previous section (we will explain below why these 
arc justified), the computation of BPS bound states in IIA reduces to computation 
of BPS degeneracies of a gas of free particles in M-theory, now living on the flat 
R 2 ' 1 world volume of the M5 brane. Below, we will show, following Ref. 11781 that 
the degeneracies of the spinning open M2 branes on R 2 ' 1 are computed by the 
open topological string, corresponding to A-model on X with a Lagrangian L. To 
extract the degeneracies of the BPS bound state of one D6 brane with the open and 
closed D2 branes and DO branes, we have to restrict the Fock space to those states 
preserving the same supersymmetry. 

7.5.1. Open M2 Branes and the Topological String 

The BPS particles arc now labeled by their U(M) representation R, bulk class 
/3, spin s and i?-chargc r. The presence of the M5 branes breaks the symmetries: 
SO(4) = SU(2) L x SU(2) R is now broken to SO(2) L x SO(2) R : these now corre- 
spond to the little group of a particle in three dimensions and the R-symmetry of 

co The counting of open BPS states was also studied in Ref. 11871 In that context, the instead of 
the D6 branes wrapping X the authors had N D4 branes on a divisor in X. 

cp In fact, as we will need later, u = f dis ^ k + i f s ± A, where A is the electric U(l) on the D4 brane. 
We can write this, equivalently as u = J disk (fc + iBjvs)i since on the D4 brane Sjvs — is the 
gauge invariant combination. 
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the M = 2 theory. More precisely sl — s + r, sr = s — r where s is the spin and r is 
the R-charge. sr permutes the fields within the same multiplet, and sl annihilates 
the state. 

We denote the number of M2-branc particles with intrinsic spin (sl,sr) and 
wrapping /3 and representation R by Ni s ^' SR ^ , and define the index 

Now, again, each such 3d particle gives rise to a field <£, and excitations of this field 
on R 2 are the particles we want to count 

$(z)=]Ta„z". 

n 

Here we have set z = z\, since M5 brane wraps the and z 2 = subspace of R 4,1 . 

The partition function of these particles can be written as follows. Introduce 
chemical potentials q, corresponding to the spin s, Q corresponding to the charge 
(3 and vt corresponding to the charges of the M2 branes under the U(1) M . For each 
field in representation R, /3, s, we get a contribution 

oo 

nn(i-9 s+ "QV s ) % > 

where the product is over the weight vectors ha are of the representation R, ni^ R 
are the corresponding multiplicities and 

M 

i=i 

The full Fock-space partition function is 

oo 

^Foc k = n n ri( i - ^-^r^*. (7.66) 

f},s,R j: R n=l 

Note that, as in the closed string case, this should include M2 branes with both 
orientations. Flipping the orientation corresponds to sending j3 to —j3, s to —s and 
R to R, simultaneously, where CPT ensures 

N%, R = NZ; >R . (7.67) 

As explained in Ref. I178[ the degeneracies of open M2 branes ending on an M5 
brane wrapping L are computed by the open topological string partition function in 
the presence of a D4 brane wrapping a Lagrangian L. It was shown in Ref. 1 1781 that 
the open topological string partition func tion h as the following simple expansion, 
similar to the Gopakumar-Vafa expansion EQEH m the closed case: 

/ 00 V d \ 
Z° ov = exp PT ]T fBiq d , Q d )^B— , (7.68) 

\d=l R J 
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where 

s,/3>0 y y 

and where Nn t p, s is an integer counting M2-brane particles with and a"l The 
symbol TtrV denotes the holonomy of the gauge field on the D4-brane. Here V 
captures the BPS masses of the open M2 branes, namely 

V = diag(ui,«2) ■ • -,vm)- 

We can write 

Tr R V d = J2m %R l[vt k ", (7.69) 

where kf" are the weights of the representation R and rn£ R their multiplicities. Then 
expression (|7.68[) can be rewritten as the Fock space trace: 

^o P =n n n(w i+ w r y^\ ^ 

n=0 R,s,f3>0 fc H 

T tt k R 

where v R := Mi'V ■ This naturally corresponds to the half of the Fock space 
corresponding to M2 branes wrapping holomorphic curves only. In other words, we 
have 

^Fock = ZZ op (q,Q > v)Z^;( q ,Q-\v- 1 ), (7.71) 
where ^top contains the contributions of the anti-M2 branes, 

z%;( q ,Q- i ,v-i)=u n n i 1 - i n - s Q-^ R ) ** r > p " , (7.72) 

n=0 R,s,P>0 j: R 

where we used the (|7.67[) and the fact that Tt^V = Tr^l/ -1 . Note that Z°^p(q, Q, v) 
is almost the same as Z° op {q, Q, v), the only difference being that the spin s in the 
power of q is replaced by — s. 



7.5.2. Open BPS Wall Crossing and M-theory 

Up to this point we discuss all M2-brane particles. Consider the central charges of 
the particles, which determine which ones are mutually BPS. The central charge of 
the open D2 brane wrapp ing a holomorphic disk in class fc^, bulk class ft and DO 
brane charge n is given bM cr l 

Z(D2) = (u(k R )+t(ft)+n)/R, (7.73) 



cq We have shifted the definition of spin by 1/2 relative to Ref. ll82l In that paper, spins were taken 
to be half integral, while here we take them to be integral. 

cr As in Ref. 1571 there is a complex proportionality constant between R in the formulas here and 
the radius of the Taub-NUT; this is irrelevant for the discussion of wall crossing. 
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where 



and 



ik + B, 



s 2 



ik + B. 

disk 

Here, B is the NSNS B-field, and k is the Kahler class and R is the radius of 
the Taub-NUT defined above. The central charge receives contributions from the 
background twisted chiral multiplets where gauge fields and their superpartners 
reside (see, for example Ref. I188p . The central charge is of the form a(q) = q ■ a 
where q is the charge of the state under the gauge symmetry, and a the lowest 
component of the twisted chiral multiplet. The first term in (|7.73|) comes from the 
U(1) M gauge symmetry on the D4 branes. The second two come from the bulk 
gauge fields, reduced to two dimension^. To satisfy our Assumption 1 (page [50]) . 
that all the central charges of D2 and DO branes align, we had set the area of all the 
disks to zero and 2-cycles to be zero. This is possible, since t and u are the moduli 
of our our theory in two dimensions, which we can dial at will. The central charge 
is then given by 

Z(D2) = (B(k R ) + B{(3) + n) /R, 

in IIA, or equivalently 

Z(M2) = C(k R ) + C((3) + n/R, 
in M-theory. The BPS particles are particles with 

Z(D2) > 0. 

This means that the complete open+closed BPS partition function is 

where the Fock-space partition function is computed by the open+closed topological 
string partition function, 



Z(D2)>0 



cs The 4d M = 2 bulk vector multiplet splits into two twisted chiral multiplets. These determines 
the 2d Af = (2, 2) central charges of the 4d electrically and magnetically charged states. Namely, 
the twisted chiral multiplet coupling to electrons contains the the d = 4 complex scalar, and 
the longitudinal components Aq^ of the d = 4 vector. The twisted chiral multiplet coupling to 
monopoles is slightly more complicated. Start with the chiral multiplet containing two transverse 
components of the gauge fields A2,3- The dual twisted chiral multiplet determines the d = 2 central 
charge of the magnetically charged particles, in particular of the D6 branc. We will need this in a 
later section. 
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7.5.3. Examples 

In Section l7T3l wc discussed a crystal melting model for the open BPS invariants, 
generalizing the previous results for the closed invariants explained in Section[3] 
There the jump of the open BPS invariants under the change of closed string moduli 
was studied. However, it was not clear there at which value of the open string moduli 
the crystal computation corresponded to, or what happens as open string moduli 
are varied. The partition function computed there had two properties: a) in the limit 
large values of closed Kahler moduli it reduced to a single copy of open topological 
string, and b) as the closed string moduli are varied, only the degeneracies of closed 
BPS states jumped. 

From our results, property a) implies that BPS invariants studied in Ref. 1621 are 
in the special chambers where 

R > 0, u l ->■ oo, 

where u l are as defined in section 3. This is because then, fixing R, and u % and 
taking Re(t) to infinity one indeed recovers the topological string amplitude. For 
example, when < Rc(i) < 1, we have 

Z BPS = Z?+ c (q,Q,v)Z c (q,Q- 1 ). 

and when Re(t) — > oo, we have 

-ZWs = Z°+ p c (q,Q,v). 

In this case, taking Re(t) to infinity in addition, only M2 branes ever contribute, and 
no anti-M2 branes. This means that Zbps has the same open part as the topological 
string theory. 

We can consider another extreme < Rc(V) < 1. In this case, open part has 
contributions both from M2-branes and anti M2-branes. For example if < Re(i) < 
1 

Z B P S (q,Q,v) = Z°+ c (q,Q 7 v)Z^*(q,Q~\v- 1 ). 

We can also consider R < 0. For example, there is a chamber R < 0, Rc(u l ) — > 
oo,Re(t) — > oo, where the partition function simply becomes one. 

7.6. Relation to the Work of Cecotti and Vafa 

The wall crossing of the closed BPS invariants discussed in Chaps. \E\ [6] is a special 
case of a more general problem, the wall crossing of BPS bound states in four 
dimensional theories with N = 2 supcrsymmctry. In this context, Kontsevich and 
Soibelman ESI recently conjectured that the degeneracies on two sides of the wall are 
related by commuting certain symplectomorphisms of complex tori. In the special 
case of = 2 gauge theories in four dim ensions, t his remarkable structure was 
explained from several different perspectives I 48 | 49 | 50 | j n p ar ticular, in Ref. IBTil both 
the appearance, and the particular choice of symplcctomorphism, was beautifully 
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illuminated. In some cases, the statements that follow from Ref. [46] are particularly 
simple, and were predicted in physics literature S3 a while back. As we saw in 
Section]!)] when the state of charge 7 decays into two primitive states 71 and 72 = 
7-7i, 

7 -> 7i + 72, 

the degeneracies of single particle states 0(7) iumrF*! as | 46 | 47 | 

n( 7 )_>n(7) + ( 7 i,7i) 0(71) nfo). 

More generally, when the central charges of 71 and 72 align, they do so for any 
multiple of 71 and 72 as well, so one can a-priori lose an infinite number of states. 
While there are no known such simple formulas for the general case, Refs. SSJ 27J 
predicted that in the semi-primitive case, where only one of the charges is primitive, 
the degeneracies fi(7i + n.72) jump as 

00 

Z(q) -> Z(q) g")-"<Ti>^M"^), 
n=l 

where we defined 

00 

Z(q) = Y / n(ji+n l2 )q n . (7.74) 

As noted in Section[5] M-theory provides a derivation of this formula when 72 carries 
DO and D2 brane charges, and 71 in addition carries one unit of D6 brane charge. 
Note that the jumps are governed by a Mac-Mahon function type formula, JI^li (1 — 

In the present context, we count open BPS states. We have seen that M-theory 
predicts that the BPS states jump as 

0(7) Q( 7 ) + 0( 7l ) fi( 72 ). (7.75) 

in the primitive case, and more generally as, 

00 

Z(q) Z(q) ]J(l-q n )- Q ^\ (7.76) 

n=l 

in the semi-primitive case, where Z{q) is as in (|7.74j) . In other words, M-theory 
predicts that the open BPS states jumps are governed by eta-function n^Li(l — 
<Z n ) _1 type formulas. 

In fact, this could have been anticipated. The open BPS states are BPS states 
in the Af = (2,2) abelian gauge theory in two dimensions on the world- volume of 
the D4 branes. The open D2 branes (bound to the closed DO, D2 branes and the 
D6 brane) , as discussed above, are charged particles in this two-dimensional theory. 

ct In what follows, wc will be cavalier about whether we gain or lose the states in the jump, i.e. 
about the sign of AH, since that depends on the direction of crossing. 
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The particles in two dimensions are solitons. This is so even for fundamental matter, 
simply because a particle is a codimension one object, so the vacuum of the theory 
can change from one side to the other. Soliton spectra of the massive J\f = (2,2) 
theories were studied in Refs. ll89ll36l The fact that the spectrum of BPS states can 
jump was in fact discovered in this context. We will now first review the essential 
results from Ref. 1361 and then explain how to apply them in the present context. 

7.6.1. Review of the Results by Cecotti and Vafa 

Consider a massive AT = (2, 2) theory. BPS particles are solitons A^ interpolating 
between the vacua i and k at spatial infinities. Let 

Mifc =Tr. lfc (-l) F F 

be the "number" of such solitons, or more precisely the index that weights the 
solitons with their fermion number charge F. Only the BPS solitons contribute to 
the index. These live in short, two dimensional representations of the N = (2, 2) 
supersymmetry, and the index is the same as counting the multiplets weighted by 
(— 1) F , where F is the fermion number of the lowest component in the multiplet^^. 

The central charge, Z(Aik) of the soliton depends only on the vacua i and k, so 
we can define Wi, so that 

Z(A ik ) = Wi-W k . 

If the theory at hand is a Landau-Ginsburg theory, then Wi is the value of the 
supcrpotcntial W in the i-th vacuum. If, as we cross the line of marginal stability, 
Wj crosses the line in the T^-plane interpolating between Wi and W k , the soliton 
A ik decays as 

and moreover the number of solitons jumps as 

IMk -> (Mk + VijUjk- (7.77) 

More generally, vacua j\ , ji , . . . j'at can cross the straight line connecting Wi and Wj , 
where the order is set by the order of interception points. The state can have 
different decay channels, corresponding to all the different ways of getting from i to 
k while passing through the intermediate vacua. The index in the ik sector jumps 
as 

51 /'■.; /', .., , •••/', . ./. • (7-78) 

l<n<N 

1 <si<s 2 <...<s„<A r 

This has a simple interpretation. The jump corresponds to counting all the chains of 
BPS solitons interpolating from vacuum i to vacuum k via the intermediate vacua. 
More precisely, we are counting the lowest components of these multiplets, weighted 
by (— 1) F . Since the net fermion number is simply the sum of the fermion numbers, 
the result is simply the product of the BPS degeneracies in the individual sectors 
of the chain. 
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7.6.2. Counting of Open BPS States 

In the case at hand, we are counting massive BPS particles, which carry charges 
under the ?7(1) gauge fields from the D4 brane and the bulk. The corresponding 
central charge that enters the M = (2, 2) supersymmetry algebra is determined by 
the lowest components a a of the twisted chiral multiplets I] J"1 For a particle 7 of 
charge q a , the central charge is 



Per definition then, the solitons of the theory are the fundamental particles. Consider 
a vacuum of the theory Adding a BPS particle 7 to it changes the vacuum of 
the theory from \i) on one side to \j) on the other where, again per definition, 



We can think of the vacua as labeled by the charges of the particles needed to create 
it, from an arbitrary but fixed vacuum. Note that we have extra structure here, that 
is not present in an arbitrary massive M = (2,2) theory Namely, given any vacuum 
\i'), the BPS particle 7 takes it to a vacuum \f) where Wy and Wv differ by the 
BPS mass of 7, W r - W v = Z(j). 

Among the particles in the theory is a D6 brane wrapping the whole non-compact 
Calabi-Yau. The focus of our paper has been understanding the wall crossing of 
particles carrying one unit of D6 brane charge. However, the corresponding central 
charge is strictly speaking infinite. If \k) is a vacuum with one unit of D6 brane 
charge, then Wk is strictly at infinity in the W plane. To deal with this infinity we 
need a regulator. The correct way to do this, it turns out, is to cut off the W-planc 
at some large but finite radius A. A^ is then a non-compact cycle with a boundary. 
The relevant data that affects the problem is only the angle 8 that Wk — Wi makes 
with the real axis in the W-pland cv l (see Fig. |35|) . Moreover, it is easy to see that 
9 is independent of any additional D2 or DO brane charges that the D6 brane may 
carry, since the central charges of these arc finite. So, any soliton Ajk where Wj 
is in \W\ < A, corresponds to a state with the same value of 9, but whose central 
charge Z(Ajk) differs from Z(A^.) by a finite amount Wj — Wi. 

Fix now 7 = Ajfe, corresponding to one unit of D6 brane charge, and some 
D2 brane and DO brane charges. The degeneracy mk of domain walls interpolating 
between the vacua i and k is simply the degeneracy of the 1-particle states with 
charge 7, 



cu In this section, for simplicity of the notation, we use a to denote both open and closed twisted 
masses of the particles on the D4 branes. These include the central charges of DO, D2 and D6 
branes, bound to open D2 branes. In addition, since we are not interested in the dynamics of the 
gauge theory on the D4 branes, but only in the BPS particle content, we have been viewing E Q 
as non-dynamical, so that the charges q a are simply the global symmetry charges. 
cv This parallels the treatment in Ref. 1531 where the wall crossing of closed BPS bound states 
corresponding to D6 branes on non-compact Calabi-Yau was studied. 




(7.79) 



Wj -Wi = Z{i). 



Hik = 0(7)- 
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Fig. 35. In the W-plane, all the Wj i 's are on a single line emanating from Wi. is infinitely 
large, and we cut off the W-plane at large radius A, shown here in green. 

As we vary the moduli, suppose vacuum j inside the cutoff W^-plane crosses the 
straight line between i and k. This then implies that corresponding state 72 is a 
single particle state, with some DO and D2 brane charges, and no D6 brane charge. 
So, Ay = 72, and Ajk =7 — 72 = 7i carries one unit of D6 brane charge. The 
soliton degeneracies are simply the degeneracies of the corresponding one-particle 
states, so 

(Mj =0(72), (ijk =^(7i), 
and the jump (|7.78[) in the degeneracies of 7 predicted from Ref. [SHI agrees with 

GZ5). 

Moreover, if the vacuum \j) corresponding to 72 crosses the ik line in the W 
plane, a large number of other vacua will cross as well. These are vacua which can 
be obtained from \i) by adding n 72 particles. They are all collinear, located at 
points Wj n = Wi + nZ (^2) ■ This holds for any n, as long as Wj n is inside the cut-off 
W-plane. These additional critical points lead to more complicated decays. Let's 
fix 71 = A.j n k, the decay product with one unit of D6 brane charge, and consider 
all the different ways that we can obtain it from 7 = A^. There is is a two-soliton 
decay, 

A lfe -> Ay„ + A jnk . 

where 7 = 71 + 7172 decays into 71 and 7172 = Ay n This, according to Ref. [36] 
corresponds to a change in the degeneracies 

or 

0(7) ^0(7) + ft( 7l Mn7 2 ). 

This is exactly as predicted from M-theory. 

There are also channels where 77,72 is split into more particles, each carrying 
some multiple of 72 charge. Here we need to be careful. Much of the discussion of 
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Ref. [36] assumes that no three vacua are collinear in the VK-plancQ. Here, we are 
in the opposite regime, of a large number of collinear vacua. In this regime, a very 
naive application of (|7.78[) overcounts the change in the soliton number. 

Recall that, from (|7.78p . the number of solitons gained or lost, is times 
the number of ways of assembling a chain of solitons of total charge 77.72 out of 
the available ^(^72) one particle BPS states particles of charge £72, where k = 
1,2.... In counting the number of possibilities, the order of in which we string the 
72-solitons in a chain would matter. When the vacua are not collinear, the order 
indeed matters. However here this manner of counting would lead to a contradiction. 
Consider example, in a decay channel where 7 splits into 71 and n copies of 72, 

A ik -> A. in + A jlh +... + A jnk , 

The n 72 particles are all mutually BPS, so in counting the number of such n particle 
states, we can simply use the free Fock-space, generated by 1 particle states. The 
dimension of the corresponding n-particle Hilbcrt space is not $7(72)", as the naive 
application of Ref. l3"B"lwould suggest, but f2(72)(fi(72) + l) • • • (f2(72)+JT.— l)/n\. The 
number of states we gain or lose is just this times fi(7i). One should in principle be 
able to verify this by a careful computation of n-soliton contributions to the index 
in Ref. I189[ l36l More generally, we can consider more complicated splits, where n 
units of 72 charge are split among mi particles of charge £72, with J^e 171 ?^ = n - 
The number of states is just the product of the dimensions of nig particle Hilbert 
spaces. Of course, this is precisely the counting predicted from M-thcory. 



8. Discussion and Future Directions 

In this paper we have seen that many problems in the BPS state counting and wall 
crossing phenomena have natural and beautiful solutions, which are obtained by 
rapid developments in this field over the last a few years. However, there are many 
interesting and important topics which require further exploration. Of course, many 
of the problems are interrelated, and the solution of one problem is likely to shed 
light on other problems as well. 



The formula we proved in Section[5] 



^BPS - ^top 



h b ' (8 ' 1} 

chamber 



is highly reminiscent of the Ooguri-Strominger-Vafa relation 

■^bps = |^top| 2 - (8.2) 

cw If collincarity is coincidental, one can remove it without loss of generalities, simply by displacing 
the vacua slightly. However in the present case it is physical, and there is no small deformation of 
the theory that can remove it. 
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It would be interesting to clarify the relation between them. Of course, one 
should keep in mind that there are important differences between the two. First, 
the former relation is verified only for the Calabi-Yau manifold without compact 
4-cycles, whereas the latter is believed to hold for any Calabi-Yau. Second, the 
first relation applies to any chamber, while the second holds only in a special 
chamber of the moduli space. In particular, as emphasized in Ref. |47l we need 
to specify the chamber of Zbps hi order to even formulate OSV conjecture. 
Third, the arguments of Zt op arc different in the two relations. Finally, we 
are typically considering different set of D-brane charges; OSV is mostly for 
D4/D2/D0, whereas (JHHJ is currently only for D6/D2/D0. 
It would be interesting to extend our analysis to a Calabi-Yau manifold with 
compact 4-cycles and with multiple D6-branes. As we have seen, it seems 
to be rather non-trivial to extend the arguments in Scction[5] and Section|6] 
to these situations. This is an important problem since Calabi-Yau geome- 
try with compact 4-cycles frequently appears in geometric engineering and in 
black hole physics, and comp arison with scmiclassical entropy formula of black 
holes (e.g. BMPV black hole H2Ql) often requires large D6-brane charges. Sec 
j^ e f s | 191 | 192 |193| £ or recen ^ discussion of multiple D6-brane charges. 

One should remember that in simple examples, the answer is given by using 
the non-commutative topological vertex in Section r7.3l For example, if we have 
a canonical bundle over P 1 x P 1 and if the B-field through one of the P l! s is 
very large, we can divide the toric diagram into two along the P 1 and glue 
the two non-commutative vertices by a propagator as in the usual topological 
vertex. However, this simply gives the topological string partition function of 
P 1 x P 1 times the wall crossing factor for C 3 /Z2. The interesting question is 
what happens when the B-fields through both P l! s become small. 

One should also remember that crystal melting model we defined in Sec- 
tion[3] can be defined for toric Calabi-Yau manifolds with compact 4-cycles. 
The only question is how that partition function is related to the topological 
string theory, where we have fewer parameters. 

Further generalization is to consider compact Calabi-Yau manifolds. This 
will lead us back to the original motivation in the introduction. This is a noto- 
riously difficult problem in general, but the recent results on the computation 
of top olog ical string partition function on a some compact example up to high 
genus I I will be helpful for this purpose. 

It is an interesting problem to generalize our crystal melting model to string 
theory compactifications other than Calabi-Yau 3-folds. In principle, when we 
compactify string theory on a non-compact manifold preserving AT > 2 su- 
persymmetry and we have BPS particle of D-branes, the computation of BPS 
index should reduce to the computation of the Witten index of the supersym- 
metric quantum mechanics arising as the low energy effective field theory of 
D-branes. The possible drawback is that in general the connection with topo- 
logical string theory will be lost. Moreover, identifying the precise form of the 
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quantum mechanics, especially the supcrpotcntial, is a difficult problem in gen- 
eral. However, a technique generalizing the branc tiling method in Scction [3~2l 
is available for M2-branes probing Calabi-Yau 4- fold singularities @, g iving rise 
to a 3d Chern-Simons- matter theories generalizing the ABJM theory < We 
should als o no te that there is a quiver description for compact examples, such 
as quintic ^^1. 

It would be an important problem to check our prediction for non-toric exam- 
ples, for example the one in Section l5.4.3l Recently the paper Ref. [89] studies 
examples C 2 /7 x C, where 7 is a discrete subgroup of SU{2). The subgroup 7 
has a ADE classification, and even for D and E type singularities they proved 
the relation 

Zncdt(<7, Q) = Z top (q, Q)Z top (q, Q~ x ), (8.3) 
which is consistent with the results of Scction[5J 

In our discussion of wall crossing from wall crossing formulas in Section[6j 
semi-primitive wall crossing formula is sufficient since the Calabi-Yau mani- 
fold is non-compact and we have only a single D6-brane. However, if we want 
to discuss multiple D6-branes on compact Calabi-Yau manifold, we need the 
general Kontsevich-Soibelman wall crossing formula. In these case the simple 
relation like (|8.1[) would presumably not hold, however it is interesting to see 
if there is still a connection of the BPS state counting problem with topolog- 
ical string theory. The work of Ref. [50] is interesting from this viewpoint. It 
is also i mporta nt to physically understand the Kontsevich-Soibelman formula, 
see Rets.™. 

It would be int erest ing to see if there is a generalization of Gopakumar-Vafa 
large N duality including the background dependence, and whether there 
is a Chern-Simons interpretation of the wall crossing phenomena (see the recent 
paper Ref. [501 for related ideas). This will potentially lead to new persepctives 
in knot theory. 

Similarities and differences between wall crossings in Af = 2 and Af — 4 theories 
are interesting problems to pursue. For example in Af = 4 dyon counting in type 
II comp actification on K 3 x T 2 , there is a universal formula for the degeneracy 
l 206 l 207 l, stating the BPS degeneracy tt(P, Q) is given by 

r Ki(P 2 p+Q 2 a-2P-Qv) / \ 

(-innp,(|./«! m , «=Q- (M 

Here P, Q are magnetic and electric charges, $(M) is unique weight 10 Igusa 
cusp form of Sp(2, Z). Moreover, the wall crossin g, phenom ena is simply inter- 
preted as the change of the choice of the contourT2QHI2Q9] is an interesting 



cx See Refs. [[95] [[96] [193 DM] ESI EE] EE] for some early papers. Ref. [202] also discuss G 2 
holonomy manifolds. 
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problem to see if we can write the Af = 2 BPS partition function in a simi- 
lar form, namely a universal expression with the moduli dependent choice of 
contour. The cusp form 1/$ can be written as a square of the character of the 
Borcherds-Kac-Moody algebra, which is highly reminiscent of our expression 



Z t 2 op in (EB). i—,— i 

A related question is to identify the BPS state algebra EEEH f or the N = 2 
theories discussed in this pap er. In the J\f = 4 case this is given by a Borcherds- 
Kac-Moody algebra EMEil anc i f or ALE spaces co rresp onding to ADE affinc 
Dynkin diagrams, by an affine Kac-Moody algebra ^. Mathematically, the 
BPS state algebra is believed to be captured by the the motivic Hall algebra 
studied in Rcf. |46] See Ref. 12151 for more along these lines. 
The D6/D2/D0 systems studied in th is paper is T -dual to the D5/D1/D(-1) 
system, discussed in a series of papers I 216 | 217 |218| E u dicl ean D-branc instan- 
tons (Dl/D(-l)-instantons in this example) produce non-perturbative terms, 
which arc of phcno menqlogical i nterest since they can generate perturbatively 
forb idden couplings | 219 | 220 |221| an( j hayg applications to moduli stabilization 
l ^22 | It was noticed in particular in Ref. 12181 that topological string theory lies 
behind the continuity of the non-perturbative term across walls of marginal 
stability, which is similar in spirit to the work of Ref . [57] explained in Section[5j 
Since there has recently been significant developments in the understanding 
of D-brane instantons, these works can provide a different perspective on the 
problem discussed in this paper, as well as phenomenological motivations. 
Crystals similar to the one discussed in this paper appear in the context of 
counting gauge invariant operators in four-dimensional Af = 1 supersymmetric 
quiver gauge theories (see for example Refs. 12251 1224|) ^|. The difference is that 
they only consider gauge invariant operators, whereas in our construction of the 
crystal we used a path algebra, an algebra consisting of all possible operators 
(regardless of whether gauge invariant or not). This means that our crystal 
contains the crystal of Ref. 12231 as a subset (sublayers). Also, in the definition 
of the partition function we are using a set of atoms (satisfying the melting 
rule) , whereas Ref. 12231 counts the atoms themselves. This is attributed to 
the difference between the Euler character and the Hilbert series of the moduli 
space. 

Despite these differences, there are some similarities and it would be in- 
teresting to understand the relation better. For example, the paper Ref. 12231 
observes that their crystal changes as they change the value of the B-field and 
go through the flop transition. This is part of the chamber structure of the 
Kahler moduli space. Also, the counting of gauge invariants operators in toric 
examples reduces to counting of the lattice points in the convex polytope (im- 
age of the moment map), which is in the same spirit of the discretization of 



cy T his counting problem is also related to the index of 4d Af = 1 supcrconformal field theories 
[2251 



space we saw in Section[4] 

As mentioned in Scction[3l the atomic structure we have discovered would have 
possible implications for quantum aspects of black hole physics. Such discrctizcd 
structures are com mon in geometric quantizati ons, and is reminiscent of the 
bubbling geometry G3Q] anc | Fuzzball conjecture H^U 

Fluctuation of the dimer model from the thermodynamic limit is an interesting 
problem. In SectionSl we have seen that the thermodynamic limit of the crystal 
partition function gives a genus topological string partition function, and the 
dependence on the wall crossing disappears in the thermodynamic limit. By 
studying the fluctuation of the crystal from the limit shape, we should be able 
to see the higher genus topological string partition function as well as chamber 
dependence of BPS invariants EE This should be a doable problem since the 
dimer m odel is a kind of an exactly solvable model whose correlation function 
is 

known [227lIl3] 

in the mathematics literature. The higher genus cont ributions 
are also captured by the matrix model ("remodeling" the B-model) I 228 | 229 |^ 
and we expect a more direct connection between the dimer fluctuations and 
the matrix model. See Refs. 1230} 12311 for recent discussion. 
In this paper we only discussed the classical version of the crystal melting 
model, or cquivalently the d imer model. In the literature, there is a quantum 
version of the dimer model which originally arises from the resonating 
valence bond models for superconductivity and is known to have an interesting 
phase structure. There is a spec ulation that quantum dimer model is related 
to the topological M-theory l ^3o [ g ee Ref . Iffifl for recent discussions. 
Many of the four-dimensional M = 2 theories arise from six dimensions, by 
wrapping M5-brane on a punctured Ricmann surface. These theories have re- 
cently attracted considerable attention due to the work of Ref. 12341 which 
generalizes the result of Ref. 12351 and proposes a new family of interacting 4d 
conformal field theories. Furthermore, a new connection between 4d ga uge t he- 



ories and 2d Liouville theory on the Riemann surface was proposed 1231] ^ 
would be interesting to explore the implications of these results to the results 
presented in this paper. See the works Refs. l49l 11851 for discussion along these 
lines. 

The list can go on, but let us stop here. We hope that future developments will 
solve some of these problems, and will ultimately lead to a better understanding of 
string theory as a theory of quantum gravity. 



Acknowledgments 

This research project would not have been possible without the help, support and 
inspiration I obtained from many people. 



cz Sec Ref. 12261 for a recent attempt. 



Crystal Melting and Wall Crossing Phenomena 107 

First and foremost, it is a pleasure to express my gratitude to Prof. Hirosi 
Ooguri. He was my adviser during the last two years of my Ph.D. course at IPMU 
and Caltech, and this paper owes much to the discussions and collaborations with 
him. I learned a lot from him, and he offered me invaluable assistance, guidance 
and support. 

I owe my deepest gratitude to my adviser Prof. Tsutomu Yanagida. During 
my stay at the University of Tokyo I often enjoyed the discussion with him about 
physics, which was always inspiring and suggestive to me. I also had a wonderful 
opportunity to collaborate with him. 

I am also deeply indebted to Prof. Tohru Eguchi, who is my mentor since my 
Master course study at the University of Tokyo. The research on brane tilings and 
quiver gauge theories, which I have undertaken under his supervision, provides the 
basic ground which lead me to the study described in this paper. 

Special thank go to my collaborators: Mina Aganagic, Kentaro Nagao, Cumrun 
Vafa. It was a wonderful and an exciting experience to collaborate with you. I 
also obtained valuable feedback from the discussion with Alexei Borodin, Tudor 
Dimofte, Lotte Hollands, Sergei Gukov, Daniel Krefl, Andrei Okounkov, Kazushi 
Ueda, Jaewon Song, Piotr Sulkowski, Yukinobu Toda. 

I would also like to express my gratitude to the Ph.D. paper committee, espe- 
cially Yoichi Kazama, for criticism and comments on the manuscript which is partly 
incorporated into this final version. 

I am indebted to my colleagues at the University of Tokyo, both at the Depart- 
ment of Physics and the Institute for the Physics and Mathematics of the Universe. 
I would also like to thank the high energy theory group at Caltech for hospitality, 
where I enjoyed my one year stay in total. I would also like to thank the secretaries, 
Mami Hara at Tokyo and Carol Silberstein at Caltech. 

The contents in this paper are presented at numerous universities and workshops 
I attended. I convey thanks to all the participants for their interest, questions and 
feedback which greatly help to improve this paper. 

During my Ph.D study I have been supported by the JSPS fellowships for Young- 
Scientists, by the Global COE Program for Physical Sciences Frontier at the Uni- 
versity of Tokyo funded by MEXT of Japan, by DOE grant DE-FG03-92-ER40701 
and by the World Premier International Research Center Initiative of MEXT of 
Japan. I would like thank these funding agencies for support. 

Last but not least, I wish to thank my parents for their understanding, love and 
support. 

Note added for the journal version: I would like to thank the anonymous referee 
for careful reading of the manuscript and for pointing out numerous typos and 
mistakes in the original version of this paper. 



108 Masahito Yamazaki 



Appendix A. Appendices 
Appendix A.l. Perfect Matchings 

In this appendix we are going to explain the one-to-one correspondence between a 
molten crystal discussed in Section|3]and a perfect matching of the bipartite graph. 
This means that the problem of counting BPS states can also be reformulated as a 
problem of counting perfect matchings of the bipartite gra ph0 Such a description 
is useful for the discussion of thermodynamic limit in Scction[4] The contents of this 
appendix is basically a recapitulation of Ref. [MJ 

In Scction r3.4.1l we considered a quiver Q = (Qq, Q\), which is a universal cover 
of the quiver Q on T 2 . The dual graph of Q, which we denote by 7, can be made 
bipartite using orientation of arrows of Q, and is a universal cover of the bipartite 
graph 7 on T 2 described in Section l3.2.2l What we are going to do is to give an 
explicit correspondence between a perfect matching of the bipartite graph 7 and a 
configuration of molten crystal. 

We first construct a perfect matching from a molten crystal. Given a molten 
crystal as shown in Fig.[36l choose all the arrows of Q which are along the surface 
of the crystal. In the example of Fig.[36l such arrows are colored green in Fig.[36l 
while the remaining arrows are colored red. Take the set of the dual of edges colored 
red. It is proven by Ref. [34] that such a subset of edges of 7 is a perfect matching. 
This is the perfect matching we wanted to construct. 




Fig. 36. Given a configuration of a molten crystal, we can construct a perfect matching of the 
bipartite graph. Each arrow is colored green if it is along the surface of the crystal, and red 
otherwise. The set of dual of arrows colored red gives a perfect matching of the bipartite graph. 



da A bipartite graph is a graph with vertices colored either black or white and edges connecting 
vertices of different colors. A perfect matching of a bipartite graph is a subset of edges of the 
bipartite graph such that each vertex is contained exactly once. See Ref. [237 for an introduction 
to the dimer model. 
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In the case when no atoms are removed from the crystal, the perfect matching 
obtained by this method is called the canonical perfect matching, which we denote 
by rriQ. Since only a finite number of atoms are removed from the crystal, the perfect 
matching obtained from a molten crystal by the above method coincides with mo 
when sufficiently away from the reference point iq. 

Conversely, given a perfect matching m which coincides with mo when suffi- 
ciently away from the reference point ig, we can reproduce a molten crystal. Let us 
superimpose m with mo, and we have a finite number of loops, as shown in Fig. [37] 
in the case of Suspended Pinched Point. Define a height function h m such that 

(1) h m {j) = when sufficiently away from iq. 

(2) h m increases by one whenever we cross the loop and go inside it. 

The example of h m for the case of Suspended Pinched Point is shown in Fig. [37] 




Fig. 37. By superimposing a perfect matching of Fig. 1361 with the canonical perfect matching 
shown later in Fig. 1381 we have a set of loops, which defines a height function h m - From this 
function we can recover a molten crystal. 

By removing h m (j) atoms from each j £ Q 0l we can construct a molten crys- 
tal. It was proven in Ref. [34] that the set of atoms removed from the crystal so 
defined satisfies the melting rule of Section l3.4.2l This establishes the one-to-one 
correspondence between a molten crystal and a perfect matching of the bipartite 
graph, meaning that BPS states can also be counted by perfect matchings of the 
bipartite graph 7. 

Finally, let us finish this appendix by pointing out an interesting connection of 
the canonical perfect matching mo with toric geometry. The example of canonical 
perfect matching mo for the Suspended Pinched Point is shown in Fig. [38] In this 
example, the asymptotic form of the bipartite graph has four different patterns. Each 
of four patterns is periodic and therefore an be thought of as a perfect matching 
of the bipartite graph on T 2 . In the brane tiling literature, a perfect matching on 
the bipartite graph on T 2 is known to correspond to one of the lattice points of the 
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toric diagram IS2HH!|db| rcc0 g mze that the four perfect matchings are identified 
with the four corners of the toric diagram in Fig.@J(a) and that the borders between 
different patterns arc identified with the blue lines in Fig.0J-(b), which makes the 
(p, g)-web of the diagram. 




v — • 



Fig. 38. The canonical perfect matching of the bipartite graph for the Suspended Pinched Point 
singularity. Asymptotically, the perfect matching corresponds to one of the four perfect matching 
of the bipartite graph corresponding to vertices of the toric diagram. The blue borders between 
different choices of perfect matchings represent the (p, g)-web. 



In general, for an arbitrary toric Calabi-Yau manifold, we can use the same pat- 
tern to construct a perfect matching. Divide the universal covering of the bipartite 
graph into segments separated by the (p, g)-web of the toric diagram^. The perfect 
matching in each segment is periodic and is identified with one of the perfect match- 
ings of bipartite graphs on T 2 , which corresponds to one of the lattice points of the 
toric diagram and the lattice point in question is precisely the vertex surrounded 
by the two (p, g)-webs on T 2 [3 This determines a perfect matching. In particular, 



For this correspondence, we consider superimposition of perfect matchings and define a Z2- 
valucd height function, which is similar to the height function h m defined previously. 
dc We choose the diagram that corresponds to the most singular Calabi-Yau manifold. 
dd In a consistent quiver gauge theory, it is believed that the multiplicity of perfect matchings at 
the vertices of the toric diagram is one EEIl 
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this means that the ridges of the crystal line up with the (p, q)-web of the toric 
diagram. 

This pattern holds for all toric Calabi-Yau manifolds without compact 4-cycles, 
and we conjecture that the perfect matching constructed in this way is canonical. 
We would like to stress again that this conjecture is not needed to construct the 
crystal melting model. Here we arc simply pointing out that, in the examples we 
have studied, the crystalline structures fit beautifully with the corresponding toric 
geometries. 

Appendix A. 2. Multi- Centered Black Holes 

In this appendix we briefly summarize the multi-centered black hole solutions dis- 
covered in Rcfs. [2381 EHS EUl I24T1 

Why do we need to bother about these black holes? In many examples in the 
literature, including the Strominger-Vafa paper the D-brane con figura tions cor- 
respond to single centered black holes. This is not always the case ' \ however, 
and the BPS in dex f or the D-branes contain contribution from multi-centered black 
holes in general Es42J l n these cases, we need to discuss multi-centered black holes in 
order to correctly account for the BPS index of the D-branes. Indeed, most of the 
BPS states discussed in this paper correspond to multi-centered black holes. More- 
over, the supergravity viewpoint provides an intuitive description of BPS bound 
states and their wall crossing phenomena, and even yields the (semi-)primitive wall 
crossing formula as discussed in Section ROl We should also keep in mind, however, 
that there is a limitation to the supergravity approach, since some of the D-brane 
configuration discussed in this paper does not have large black holes and therefore 
validity of supergravity approximation is not guaranteed. 

Multi-centered black holes are specified by a set of charges 7^ at the "centers" 
Xi of the black hole . The Ansatz for the metric is of the form 

ds 2 = -e 2U (dt + u) 2 + e- 2U dx\ (A.l) 

where the function U(x) is a warp factor and u> is a one-form. They are deter- 
mined by BPS equations of motion. In this appendix we neglect higher derivative 
corrections to the equations of motion. 
The first BPS equation of motion states 

2e- u lm(e- ia n) = -H, (A.2) 

where a = Arg (J^i^ili)) and the ft <~ e~ lt is the period vector with suitable 
normalization S3 Another equation is 

* 3 dto = (dH,H), (A.3) 

where *3 is the Hodge star in R 3 , H is a function with values in the even cohomology 
(specifying D-brane charges) of the Calabi-Yau manifold, and ( , ) is the symplectic 
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paring of the even cohomology, as in the main text. The function H, which deter- 
mines tu through (IA.3I) . is a harmonic function with poles at the position x% of the 
charges. For A r -centered black holes it is given by 



(A.4) 



The equation (|A.3[) has a nonsingular solution if and only if the following integra- 
bility condition is satisfied: 

N 



E 

.3=1, JT^ 



\xi — X j I 



2Im(e- M Z( 7i )) 



where we acted d*s on both sides of (|A.3j) and used the relation 

z( 7l ) = <7i,n>- 

In particular, for 2-centered black holes this simplifies to 



ri2 ■= Fi - x 2 1 = 



(71,72) 



2Im(e- M Z(7i)) 



(7i,72> \Z{ji) + Z(j 2 )\ 



(A.5) 



(A.6) 



(A.7) 



2 Im(Z( 7l )Z( 72 )) ^ 

Since the distance between the center should clearly be positive, we have the con- 
dition 

(7i,72>Im (Z( 71 )Z( 72 )) > 0. (A.8) 

which is sometimes called the Denef stability condition 0. Wall crossing occurs 
when r\2 goes to infinity. This is the case when 

Z( 7 i) = cZ(<y 2 ), c € R_|_. 

This coincides with the result of Section [2~2l 

Finally, the black hole carries the angular momentum 

X q " X q 



(A.9) 



Ji 



In particular, for a two-centered configuration, 



I I 

(7i,72>|- 



(A.10) 



(A.11) 



This is a classical result, but there is a quantum correction 12421 w hich modifies the 
result to 



J12 = ^(|( 7 i,72)| - !)■ 



(A.12) 



dc In general multi-centered black holes, the integrability is not enough to ensure the existence of 
physically interesting solutions. For example, we have to check that the discriminant is positive, 
and that Kahlcr moduli remain in the physical moduli space. In general, it was conjectured in 
Ref. I47l that the components of the moduli space of multi-centered black hole solutions is classified 
by attractor flow trees. 
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This angular momentum plays crucial roles in the derivation of the primitive wall 
crossing formula in Section RJTSl 

Appendix A. 3. Kontsevich-Soibelman Formula 

In this appendix we briefly comment on the more general Kontsevich-Soibelman 
formula. Since the formula itself is not needed in the rest of this pape r, the d iscussion 
is brief and the interested readers are referred to the original papers 1 1 1 for more 
information. 

Let us start from the charge lattice F, which is the even coholomogy of the 
Calabi-Yau manifold. We have a natural symplectic pairing ( , } on T. Consider a 
Lie algebra with generators e 7 ,7 € F, satisfying the commutation relations 

[e 7l ,e 72 ] = (-l)<^( 7l , 72 >e 7l+72 . (A.13) 

We also define a multiplication by 

e 7l e 72 — e 7l + 72 . (A. 14) 

We define a group clement 

[/ 7 := cxpj^ — e„ 7 . (A. 15) 

n— 1 

which acts on e 7 > by an adjoint action. Using the commutation relation (|A.13[) . we 
have 

C/ 7 : e Y (-l) <7 ' 7 ' > e 7 ) <7 ' 7 ' ) e r . (A.16) 



As we have seen in Section l2T2l on the generic walls of marginal stability the 
decay pattern is 7 — > N171 + ^272 with Z (71) / Z (72) € R+. On both sides of the 
walls of marginal stability consider the product 

A:= (A.17) 

7=iVi7i+JV 2 72, A r i>0,A r 2 >0 

where the ordering of the product is given by the clockwise ordering of the central 
charge Z{^). When we cross the wall of marginal stability, the ordering of the 
central charges changes, and also the BPS invariants f2( 7 ;t) jumps. However, the 
wall crossing formula states that the whole product (|A.17|) stays invariant: 

A(t+)=A(t-). (A.18) 

where t± is the value of the Kahler moduli on different sides of the wall. 

Let us derive the primitive wall crossing formula from the Kontsevich-Soibelman 
formula. For a decay 7 — > 71 + 72 with 71,72 primitive, the wall crossing formula 
says 

jjQ(-y 1 ;t + )jj^hi+l2;t + )jjn( l2; t + ) _ jjn(7 2 ;t_)jj^(7i+72;t-)jjn( 71 ;t_) ^ 19) 
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Since what matters here is only the part generated by e 7l , e 72 and e 7l + 72 , we con- 
sider a quotient of the Lie algebra with all other generators. It is easy to see that this 
is consistent with the commutation relation (|A.13[) . We have f2(7,;t+) = il(ji;t^) 
for i = 1,2 because 712 are primitive. We simply use the notation f2(7i) in the 
following. We also note that 71 +72 is a central element of the quotient Lie algebra. 
This means that 

r/ AO( 7l + 72 ;t + ) = r/ 0(7 2 )r / 0(7 1 ) r/ -f2(7 2 ) r/ -f2( 7l ) 
71+72 72 71 72 71 

Using the commutation relation (|A.13|) again we have 

Afi( 7l + 72 ) = (-l)<^.-»> (71,72)0(71)^(72)- (A.20) 

This is (up to an overall sign, which depends on the definition of Afi) the primitive 
wall crossing formula (16.11) . 

It is also straightforward to derive the semi-primitive wall crossing formula from 
the Kontsevich-Soibelman formula 0. 

r7n(7i;t + ) 7-7^(71 +72 ;t+)rrO(7i+27 2 ;t + ) TT Q(-y 2 -t + ) Tr Q(2j 2 ;t + ) 

U ll U 7l+72 U 7l+272 ■•■ U 72 ^272 ■•■ (A 211 

_TTQ(r/a;t-) 770(73,*-) 770(71+272;*- ) 7-70(71+72;*- ) 770(71 ;t_ ) 

— U 72 ^72 ^71+272 U 7l+72 U 7l 

In this case, we consider the quotient algebra by e_/v l7l+ jv 2 72 

with Ni > 1. This 

means that e 7l +./v 272 's commute each other in the quotient algebra. Also, since 72 
is primitive 0(^72) does not have moduli dependence (at general value of Kahlcr 
moduli) . This means that 



770(71 ;i + )rrf2(7l +72;*+) 77^(71 +272;*+) 
U 1i U 71+72 ^71+272 



(A.22) 



= ^n( 72 )jjn(2 72 ) ^ [r n(7 1 ;t_) [r n(7 1 + 72 ;t_) /^n(7 2 ) [r n( 72 ) \ 1 
Since 

^0(71^)^(71+72^) ... =exp ^ n(7l + fc 72 )e 7l+fc72 I (A.23) 

\k>0 J 

in the quotient algebra, (|A.22[) with (|A.16[) gives the relation 

v -r-r / ,1 \ \ i (7l,72)0(i7 2 ) 

£fi( 7 i + k l2 ;t + )e 11+kl , = J] (l - (-l)%^> e , 72 J £ Q( 7l + k l2 ; t_)e 11+kl2 

k>0 l>0 k>Q 

(A.24) 

If we remove the common e 7l and replace e 72 by a formal variable p, we have 

x-^ 1 -r-r ( 1/ \ A '(7ii72)0(i72) _ — , 

XI 0(71 + ^72; t+) P k = [] (1 (-ly^V) ^ 0(71 + k l2 -t_) P k . 

k>0 l>0 k>0 

(A.25) 

df Thc following arguments are known to experts, but as of this writing there is no published 
account to the best of the author's knowledge. 
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Now recall that 0(71 + ^72) with k > is zero on one side of the wall. This means 
^(71 + ^72; t-) = for k > and we have 

n( 7l ) + X] A ^(7i + k l2 ;t + ) P k = H(i (-i)^V ) n( 7l ), 



fc>o ;>o 



(A.26) 



which coincides with the semi-primitive wall crossing formula (|6 . 2[) . This is what 
we wanted to show. 



Appendix A. 4. The Parametrization of Chambers by the Weyl 
Group 

In this section we explain the parametrization of chambers of closed BPS invariants 
by maps a and 9, as claimed in the main text. 

The map 9 is defined to be a map from the set of half-integers 1h to itself 

6 : Zh — > In-, 

satisfying the following two conditions. First, 

9(h + L) = 9(h) + L (A.27) 

for any h € Z^. In other words, 9 is periodic with period L. Second, 



Therefore is specified by L — 1 (half-) integers, namely L half-integers 

0(1/2),..., 1/2), 

subject to one constraint (|A.28j) . Let us assume for the moment that 9 satisfies the 
condition 

i)«(§) <...<» (i-1). (A.29) 

We will discuss other cases later. 

Given er, we have a specific choice of resolution having L — 1 P^s. Moreover, 
given a map we can determine the corresponding value of the B-ficld Bg by 

[B 6 (pn + ... + Oj)] = lH m e z I - V2) < »»L < 0(j + 1/2)}. (A.30) 

It is easy to see that this gives well-defined values of the integer parts of the B-field, 
which parametrize the chamber as explained in the main text. Conversely, it can 
also be proven that given any B-field, we can find a corresponding 9 uniquely 0. 

dg This should be replaced by 0(71 + £72;*+) = depending on the sign of (71,72)- The Denef 
stability condition l)A.8|l tells us on which side of the wall BPS bound states are stable. 
dh This comes from the fact that the action of the affinc Weyl group on the space of B-fields is 
faithful. 
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When the condition (|A.29[) is not satisfied, we need to change the choice of the 
crepant resolution. Choose a permutation S of {1/2, . . . ,L — 1/2} such that 

•Ki)) <»(=(§))< <»( s H)) 

holds. Then we replace 9 by 9' := 9 o £ and we choose cr' so that <j o 9 = a' o 9' 
F*l Note that the combination a o 9, which appears in the definition of the vertex 
in Section [7.3[ remain invariant under this process. This means that sometimes 
different a, 9 and a', 9' corresponds to the same chamber. We can either fix a and 
change 9 to parametrize chambers, or change both a and 9 for convenience. The 
latter parametrization is redundant, but sometimes useful. 

In the above discussion 9 appears somewhat artificially, but 9 is often used 
in the mathematical literature. The reason is that the maps 9 makes a group, 
which is the Weyl group of A^^\. As is well-known, the Weyl group of Al-i is 
the L-th symmetric group, which is a set of isomorphism {1/2,..., L — 1/2} — > 
{1/2, . . . ,L — 1/2}. The map 9 gives a generalization to the affine ca se. Th e affinc 
Lie algebra A^^i appears in the formula for the BPS partition function anc ] as 
we have seen specifies a chamber structure. This is reminiscent of the appearance 
of the W eyl group of the Borchcrds-Kac-Moody algebra in J\f = 4 wall crossing 
| 212 | 213 | would be interesting to explore this point further (see Section[8] for more 
on this). 

Finally, let us illustrate this parametrization with examples. Consider the re- 
solved conifold (L = 2) with the resolution cr(l/2) = +,cr(3/2) = — . Due to the 
condition ()A.28[) , we can write 

9(1/2) = 1/2 — N, 9(3/2) = 3/2 + N, (A.32) 

and this integer N parametrize the chambers. This integer iV is the same integer 
N appearing in (|7.17[) . When N > 0, condition (|A.29|) is satisfied and we are in one 
resolution a. When N < 0, (| A.29I) is not satisfied and by a flop transition we are in 
a different crepant resolution specified by a 1 (1/2) = — , a 1 (3/2) = +. If we define 9' 

by 

9' (1/2) = 1/2 + N, 6' (3/2) = 3/2 — N, 

then a, 9 for N > and a', 9' for N < parametrize the same chamber. In this 
sense, we can either choose a to be fixed and change 9, or change both a and 9, 
although the latter is a redundant parametrization. The chamber corresponding to 
topological string theory for one resolution a is given by N — > oo in (|A.32[) , and 
another resolution a' given by TV — > — oo. 

In fact, we can take cr 060Y;- 1 otf- 1 as cr, where 6 is the permutation induced by 9. 
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As a next example suppose L = 3. When 

9(1/2) = -5/2, 0(3/2) =3/2, 0(5/2) = 11/2, 

we have 

[B 9 ( ai )] = [5 e (a 2 )] = 1, [B 6 (ai + a 2 )] = 2. 

This is also given by 

0(1/2) = 11/2, 0(3/2) =3/2, 0(5/2) = -5/2. 

However, they parametrize different chambers in general 0. This is because the 
two 0's are related by a permutation S : {1/2.3/2,5/2} i-» {5/2,3/2,1/2}, and 
correspondingly we have to change the choice of crcpant resolution a by a o E as 
mentioned around (|A.31|) . 
More generally, if we have 

e (i) = i_ 2il _ i2 , 9 (5)^ + il _ l2 , e (£) = | + il + 2i , 

then 

[B e (ai)]=Li, [.80(0:2)] = £2, [-Bg(ai + a 2 )] = L x + L 2 , 
and if we have 

«(i)=-5" 2i '-^ «(£)=f^-^ »(f)-£ + i. + u* 

then 

[Bs(oO] = ii, [Stf(a a )] = £2, [B 9 (qi+q 2 )] =Li+L 2 + 1. 



Appendix A. 5. Young Diagrams and Maya Diagrams 



A partition (Young diagram) A is a set of non-increasing positive integers A = 
(Ai, A 2 , . . .), Ai > A 2 > . . . > 0. Given a partition A, its transpose A* is defined by 
A*j := #{j|Aj > i}, and its length |A| by |A| := Yli^i- For example, the transpose 
of (4, 2, 1) is given by (3, 2, 1, 1), and both have length 7. 

As is well-known, this is equivalently represented by a Maya diagram, namely a 
map 

A : Z h -> {±1} 

such that X(h) = ±1 for ±h 3> 1. We sometimes represent A by 

A=-- A(-5/2) A(-3/2) A(-l/2) A(l/2) A(3/2) A(5/2) 



where the symbol 



represents the position of the origin. For notational simplicity, 



we use the same symbol A for a Maya diagram. Our convention is shown in Fig.| 



^Sometimes they give the same chamber. This happens, for example for C 2 /Z2 X C, where there 
is a unique choice of crepant resolution. 
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Fig. 39. The convention of the Maya diagram in this paper. 



For example, 



rP = - 



+ H h + + • 

H h + + + ' 



(A.33) 
(A.34) 



For a Young diagram and a positive integer M, define the quotients A,(i = 
1,...,M) by 

Xi{h) = X(i - 1/2 + (h- 1/2)M) for ft G Z h . (A.35) 

As an example, let us consider A = shown in (|A.34[) . If you take M = 2, 
then (M = 2)~quotients are given by 



A x = ••• 
A 2 = • • • 



+ + + 



Suppose that M is chosen such that the representation A^ is trivial for all i 
1, . . . , M. This means that Xi can be written as 



X l {h) = %{h + MN(i)). 
The integers N(j) are called M-cores of A, and satisfies 



(A.36) 



M 



For example, if we take M = 3 for (|A.34I) . 

Ai = h + + + + + 

A 2 = h + + + 

A 3 = 

and 

N(l) = -1, N(2) = l, AT(3)=0. 
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Appendix A. 6. Proof of the Closed Expression for the Vertex 

In this appendix we give a proof of (|7.46[) . First, the following is clear from the 
definition: 

Proposition: Let a, 8 to be maps specifying the chamber for a Calabi-Yau manifold 
X. Choose an integer M and a', 8' for X 1 = X/Zm such that the following condition 
holds: 

ao8 = a'o8', fio9 = fi'o9'. 

Then we have 

Qo-',e';u',i/')(<?'; Q')\ f> = a' ' =0'°' =■■■(]">' = C(cr,0;u,I/)(9) Q)' 

This is because the both sides of the equation are defined by the same crystal 
with the same weights. The crystal for an orbifold is the same as the original crystal, 
the only difference being the difference of the weights; the crystal for the orbifold 
has more colors (variables). However, in the above equation we are specializing the 
variables so the weights are also the same. 

Lemma: For any a, 8, A, we can take M, a' , 8' such that 

a 08 = a' 08', /xo0 = 0o0'. (A.37) 

Proof. We choose an integer M such that all the M -quotients EB of /ij's become 
trivial, i.e. ML-quotients of the combined representation /j, (see (|7.16p ) is trivial. 
For example, this is satisfied if we define 

h- := min{/i €E Zh | fi o 8(h) = +}, h+ := max{/i S | /i o = — }, 

and take M so that ML > h+ — h- 0. This means that for any half- integer 1/2 < 
/i < ML - 1/2 we can take AT(/i) £ Z such that 

(N < N (h)), 



H o 8(h + NAIL) = 



\+ (N>N(h)) 
$(h+(N -N(h))ML). 



In other words, N(j) is the ML-core|ffj of fx 



Therefore the second condition of (|A.37j) holds if we define 8' : Zh — > Zh by 

8'(h) = h- N(h)ML 
for h G Z/i 0- It is clear that the first condition of (|A.37| determines er' uniquely. □ 



dk See [Appendix A.5| for the definition of M-quotients. 

dl The choice of M is not unique. But the final result is independent of the choice of M. For 
practical computation it is useful to take the minimum M. 
dm Sec [Appendix A.5| for the definition of M-coress. 

dn For practical computations, it is useful to further perform a permutation to 9 such that J A.29l > 
holds. Sce |Appcndix A~4| 



120 Masahito Yamazaki 



Our theorem follows from the Proposition and the Lemma given above. 
Theorem: For er, 9, A, take M, a', 9' as above. Then we have 



c. 



It is straightforward to generalize this theorem to the case of refined BPS in- 
variants discussed in Section r7.3.3l 
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